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Abstract 



The Einstein equations for a perfect fluid spatially homogeneous space- 
time are studied in a unified manner by retaining the generality of certain 
■ parameters whose discrete values correspond to the various Bianchi types 

^ ' of spatial homogeneity. A parameter dependent decomposition of the 

, metric variables adapted to the symmetry breaking effects of the non- 

. abelian Bianchi types on the "free dynamics" leads to a reduction of the 

equations of motion for those variables to a 2-dimensional time dependent 
Hamiltonian system containing various time dependent potentials which 
are explicitly described and diagrammed. These potentials are extremely 
useful in deducing the gross features of the evolution of the metric vari- 
ables. 



o 

u '. 1 Introduction 

bp. 

Although interest in spatially homogeneous cosmological models peaked in the 
early seventies, it was not until the late seventies that a more unified picture of 
. , the dynamics of these models developed, following the proper recognition of the 

■ role played in the problem by the gauge freedom of general relativity [1-3] . Since 

a number of books [4,5] and review articles [6-9] exist which discuss spatially 
homogeneous cosmology at various levels and from various points of view, it 
seems appropriate here to emphasize certain aspects of the subject which are not 
well covered in the literature. Remarkably there is still something new left to say 
on this topic nearly a decade after its most active period of research. The present 
discussion will seek to generalize many concepts which have already appeared 
in the context of particular symmetry types or special initial data and fit them 
together into a single unified picture of spatially homogeneous dynamics. It 
should be emphasized that although spatially homogeneous cosmological models 
are usually studied for a (nearly) discrete set of parameter values corresponding 
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to the (nearly) discrete set of Bianchi types [10 14], both the metric and field 
equations depend analytically on a 4-dimensional space of essential parameters 
(of which at most three may be simultaneously nonzero). By varying these 
parameters continuously, one may deform each of the various symmetry types 
into each other and thus relate properties of one Bianchi type to those of another, 
the more specialized symmetry types occurring as singular limits of more general 
types. 

Lagrangian or Hamiltonian techniques enable one to associate a finite di- 
mensional classical mechanical system with the ordinary differential equations 
equivalent to the spatially homogeneous Einstein equations and thus offer a con- 
venient means of visualizing the dynamics and of understanding its qualitative 
features. These techniques, which provide the framework of this exposition, were 
pioneered by Misner [15 18] and followed through by Ryan [19 23], leading to an 
alternative but equivalent description [24] of the qualitative results obtained by 
Lifshitz, Khalatnikov and Belinsky [25-29] for the evolution of certain spatially 
homogeneous cosmological models near the initial singularity using piecewise 
analytic approximations. This latter work was later confirmed and extended 
by Bogoyavlensky, Novikov and Peresetsky using powerful techniques from the 
qualitative theory of differential equations [30 35]. Qualitative studies of vari- 
ous perfect fluid models including the regime away from the initial singularity 
should also be noted [36-40]. 

The present disciission has as its foundation previous papers of the author 
[1-3, 41-43]. No attempt will be made to review the large body of important 
work which preceded them. References [4-9] adequately serve this function. In 
particular the bibliographies of the Ryan-Shepley book [5] and the recent review 
by MacCallum [9] provide an exhaustive list of relevant research papers. 

It turns out that the key to our problem is a simple adage of mathemati- 
cal physics: whenever a symmetric matrix is encountered, diagonalize it. First 
one diagonalizes the symmetric tensor density associated with the structure 
constant tensor whose components completely determine locally the type of 
spatial homogeneity. The three diagonal components plus an additional pa- 
rameter characterizing the trace of the structure constant tensor are the four 
parameters referred to above. Next one diagonalizes the component matrix of 
the spatial metric, maintaining the values of these parameters. This leads to 
a decomposition of the gravitational configuration space variables (namely the 
component matrix of the spatial metric) into two sets of three variables, one 
set parametrizing the diagonal values of the metric component matrix which 
are readily interpreted in terms of the action of the 3-dimensional scale group 
(independent rescaling of the unit of length along orthogonal directions) and 
another set specifying the diagonalizing matrix. This decomposition incredibly 
simplifies the field equations since the diagonalizing variables correspond to pure 
gauge directions, reflecting the effect on the metric variables of spatial diffeo- 
morphisms which are compatible with the spatial homogeneity. Scalar functions 
such as the spatial scalar curvature which are gauge invariant depend at most on 
some of the diagonal variables, for example. Finally one diagonalizes the DeWitt 
metric [44] on the conflguration space A4 of spatial metric component matrices 
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by choosing an orthogonal basis of the Lie algebras of the scale group and of the 
3-dimensional group G used in the metric diagonalization; the DeWitt metric 
is important since the "free dynamics" is equivalent to geodesic motion for this 
metric. The advantages of diagonal matrices over general symmetric matrices 
are obvious; all matrix operations (multiplication, determinant, inverse) become 
trivial and functions of these matrices have a much simpler dependence on the 
individual components. Diagonalizing a quadratic form kinetic energy function 
also greatly simplifies the equations of motion. 

A familiar example from classical mechanics which proves useful as an anal- 
ogy is the problem of the motion of a rigid body [45,46]. Such an analogy was in 
fact first introduced for general spacetimes by Fischer and Marsden [47] in their 
original discussion of the role of the lapse and shift in the three-plus-one formu- 
lation of the Einstein equations. It is even more appropriate for the spatially 
homogeneous spacetimes where the correspondence is nearly complete. The 
usual synchronous gauge spatial frame is analogous to the space-fixed axes in 
the rigid body problem. The "diagonal gauge" spatial frame which diagonalizes 
both the spatial metric and the symmetric tensor density associated with the 
structure constant tensor corresponds to the body-fixed axes which diagonalize 
the moment of inertia tensor. The special orthogonal group S'0(3, R) gener- 
alizes to the relevant 3-dimensional diagonalizing matrix group G. Since the 
components of the structure constant tensor must remain fixed under its action, 
this group is a matrix representation of a subgroup of the special automorphism 
group of the given Bianchi type Lie algebra (the subgroup is unimodular since 
its Lie algebra is required to be offdiagonal). The concept of angular velocity 
also has an analogue which is closely related to the shift vector field whose as- 
sociated time dependent spatial diffeomorphism drags the synchronous spatial 
frame into the diagonal gauge spatial frame by inducing the time dependent 
frame transformation which orthogonalizcs the spatial frame. However, since 
the action of G on the metric configuration space represents an orbital motion, 
the situation is more involved. 

At this point it is helpful to keep in mind the problem of the nonrelativis- 
tic motion of a particle in a spherically symmetric potential (the central force 
problem). Here the symmetry group SO{3,R) is a subgroup of the group of 
motions of the Euclidean metric on the configuration space . By introducing 
spherical coordinates one separates the configuration space variables into angu- 
lar variables describing the orbits of S0{3,R), namely 2-spheres except for the 
fixed point at the origin where the orbit dimension degenerates, and radial vari- 
ables which describe the directions orthogonal to the orbits. The components of 
orbital angular momentum arise from evaluating the moment function [40] for 
the action of SO{3, R) on R^ in the standard basis of its Lie algebra and may be 
interpreted as the inner products of the standard basis of rotational Killing vec- 
tor fields with the velocity of the system. Since 50(3, R) is a symmetry group of 
the dynamics, angular momentum is conserved and the problem is then reduced 
to 1-dimensional radial motion in a new potential, the angular momentum con- 
tributing an effective potential (the centrifugal potential) to the original radial 
potential. When the latter potential is absent, the case of the motion of a free 
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particle, it is of course simplest to consider only radial orbits which lead to the 
simplest representation of the free (straight line) motion, namely geodesies of 
the Euclidean metric. This restriction to radial orbits is possible because of the 
additional translational symmetry which allows one to transform the angular 
momentum to zero. 

In spatially homogeneous dynamics the Euclidean metric on is replaced 
by the Lorcntzian DcWitt metric on the 6-dimensional space A4 of spatial met- 
ric component matrices. The decomposition of the Euclidean space variables 
goes over roughly into the 3-dimensional space of "diagonal" variables and the 
3-dimensional space of "offdiagonal" variables as described above. The offdi- 
agonal variables describe the orbits of the action on ^A of the matrix group 
G and the diagonal variables describe the orthogonal directions. The moment 
function for the action oi G on A4 is the analogue of the orbital angular mo- 
mentum. Its components in a certain basis of the matrix Lie algebra ^ of G will 
be seen below to correspond to the space-fixed components of the spin angular 
momentum in the rigid body analogy, thus neatly intertwining these two classi- 
cal analogies. The free motion (geodesies of the DeWitt metric) is most easily 
represented as purely diagonal, using the larger isometry group SL{3, R) of the 
DeWitt metric to transform away the angular momentum associated with any 
particular subgroup G. The overall scale of the metric matrix represented by 
its determinant (product of its diagonal values) corresponds to the single time- 
like direction and the free motion is subject to an additional energy constraint 
requiring the geodesic to be null. Apart from the freedom to rescale and trans- 
late the affine parameter of these diagonal null geodesies, there is a 1-parameter 
family of them, parametrized by the angle of revolution of the 2-dimcnsional 
null cone in the space of diagonal metric matrices; these are the well known 
Kasner solutions [48]. However, a geodesic which has zero angular momentum 
for a particular choice of the group G will have nonzero angular momentum for 
almost all other choices of this group. 

In the central force problem the separation of radial and angular variables 
is clean. In our problem the symmetry group of the free dynamics is SL{3, R) 
and only the corresponding division of variables into a conformal metric (unit 
determinant) and a scale variable (the metric determinant which parametrizes 
the orbits of SL{3,R)) is clean. The division of variables into two orthogo- 
nal sets of three diagonal variables and three offdiagonal variables is instead 
highly ambiguous. There is essentially a 2-parameter family of subgroups G of 
SL{3, R) with 3-dimensional offdiagonal matrix Lie algebras whose orbits are 
almost everywhere transversal to the diagonal submanifold of Ai. The orbits 
of any of these subgroups may be used to perform the decomposition, which 
is automatically orthogonal with respect to the DeWitt metric. The addition 
to the free system of the spatial scalar curvature as a potential breaks down 
the SL{3, R) symmetry uniquely to one of those subgroups in the general case, 
although some degeneracy remains in some of the more specialized symmetry 
types (Bianchi types II and V; Bianchi type I is the free system and the symme- 
try remains unbroken). This breakdown of SL{3,R) symmetry to a particular 
subgroup G depends continuously on the parameters which specify the Lie al- 
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gebra of the spatial homogeneity group, just as the symmetry breaking scalar 
curvature potential itself depends continuously on these parameters. When the 
subgroup G has a compact subgroup, the diagonal/offdiagonal decomposition 
develops a singularity where the orbit dimension degenerates from its generic 
value three, similar to the coordinate singularity of spherical coordinates at the 
origin of R^. These points of the configuration space turn out to be associated 
with additional continuous symmetry of the spatial metric; they are protected 
by angular momentum barriers in the same way as is the origin of B? in the 
central force problem. 

The concept of angular momentum links the rigid body and central force 
analogies. The time dependent diagonalizing matrix is a curve in the matrix 
group (j, representing the orbital motion of the system. (In fact the matrix 
group G directly parametrizes the points of each orbit.) The tangent vector 
of this curve, namely the velocity associated with the offdiagonal variables, 
represents the orbital velocity of the configuration space point. Given a basis 
of the matrix Lie algebra g of G, one determines a corresponding basis of the 
tangent space at the identity and two global frames on G, one left invariant 
and one right invariant, which reduce to this basis at the identity. The left 
and right invariant frame (contravariant) components of the velocity tangent 
vector correspond respectively to the space-fixed and body-fixed components of 
the angular velocity in the rigid body analogy. These components are related to 
each other not by G itself as in the previously defined space-fixed and body-fixed 
components but by the adjoint representation of G. The association of left and 
right with space and body assumes that the diagonal gauge frame is related to 
the synchronous gauge frame by a passive transformation, corresponding to a 
right action of G on and since G is identified with its orbits this becomes 
right translation of G into itself. 

By the local identification of G with its orbits in M (the "offdiagonal vari- 
ables"), one may use the DeWitt metric on the orbit to lower the indices of the 
velocity tangent vector, leading to what are analogous to the space-fixed (left 
invariant frame) components and body-fixed (right invariant frame) components 
of the angular momentum. Thus the frame components of the DeWitt metric 
along the orbit act as the components of the moment of inertia tensor in the 
rigid body analogy. The "space-fixed components of the angular momentum" 
are the inner products of the left invariant frame vectors with the velocity. These 
vector fields generate the right translations and hence the right action of G on 
M and are therefore Killing vector fields of the DeWitt metric. The space-fixed 
components of the angular momentum are thus the components of the moment 
function for the action of G on and therefore the components of the orbital 
angular momentum, which are conserved for the free motion. Note, however, 
that the body-fixed components are related to these constant components by 
the adjoint transformation and so are in general time dependent. Furthermore, 
the right invariant frame components of the DeWitt metric must be indepen- 
dent of the orbital variables since right translation is an isometry. By properly 
choosing the basis of g, these components may in fact be diagonalized. 

One advantage of the present problem over the analogous classical problems 
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is that one is free to reparametrize the time variable by introducing a nontrivial 
spatially homogeneous lapse function. Lapse functions which depend only on 
the metric component matrix correspond to conformally rescaling the DeWitt 
metric and the scalar curvature potential. However, in general the lapse function 
may also depend on time derivatives of the metric components leading to a 
much larger freedom in the Hamiltonian system (a freedom not permitted in the 
Lagrangian approach) . Often special choices other than the usual cosmic proper 
time are suggested by the dynamics which help to simplify its description. Two 
such choices are the Misner fi-time [15] related to the logarithm of the metric 
determinant and his supertime [17], where the lapse is simply related to the 
metric determinant. The latter choice of time is also crucial to the Belinsky- 
Lifshitz-Khalatnikov analysis of the dynamics near the initial singularity. For 
the free dynamics, they are both affinc parameters for the geodesic motion. 

Of course all of these remarks will become clearer once explicit notation and 
formulas are introduced. The result of this formal manipulation is a reduction 
of the Einstein equations to a 2-dimensional Hamiltonian system with time 
dependent potentials associated with the spatial curvature, with the centrifugal 
forces arising from the "motion" of the diagonalizing spatial frame, and with 
the energy-momentum of the soiircc of the gravitational field, here assiimed 
to be a perfect fluid. This system must be supplemented by the equations of 
motion of the source of course. Explicit diagrams of the various time dependent 
potentials are extremely useful in deducing the gross features of the evolution of 
the metric variables; near the initial singularity they may be used to construct 
"diagrammatic solutions" of the field equations, as done by Ryan for the Bianchi 
type IX case [20]. 

The main body of the paper is divided into three sections. In the first of 
these the parametrized spatially homogeneous spacetime and field equations are 
introduced. In the second the parametrized decomposition of the metric vari- 
ables is introduced and used to reduce the Einstein equations to a 2-dimensional 
time dependent Hamiltonian system. The potentials of this system are then de- 
scribed in detail. In the third section their qualitative effect on the dynamics is 
discussed. 

2 The CD-parametrized Spatially Homogeneous 
Spacetime 

Before even beginning to discuss spatially homogeneous spacetimes, it is worth- 
while introducing some useful facts concerning the smooth action of a Lie group 
G on a manifold M as a transformation group [49-52]. A left action is just a 
homomorphism f : G ^ V{M) from the Lie group into the group of diffeo- 
morphisms of M, i.e. for a e G the corresponding transformation fa satisfies 
under composition fa^ o /a^ = faia^- By defining = one obtains an 
antihomomorphism /"^ satisfying o = /j^^^ which is the defining re- 
lation for a right action. When / is an isomorphism (G ~ /g = {/a | d € G}) 
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so that only the identity ag G G acts as the identity transformation on M, / is 
called an effective action. For example, any Lie group acts effectively on itself 
on the left by left translation La{a\) = aai and on the right by right transla- 
tion Ra{ai) = aia with G ~ Lg - Rg- (Note that R'^ : G V{G) is an 
isomorphism and a left action.) Introducing the redundant but useful notation 
a ■ X = fa{x) for the transformation /„ acting on a; e M, denote the orbit of x 
by G • .T = {a • .X I a e G}, namely all points which can be reached from x under 
the action of the group. 

Not all of the transformations are effective in moving the point x. Let 
Gx = {a G G I a • .T = .t} be the isotropy subgroup at x of the action of G on M , 
namely the subgroup of G which leaves x fixed. Intuitively one expects that at 
least locally the orbit of x is in a one-to-one correspondence with the smallest 
subset of transformations which can move the point x arbitrarily on its orbit. 
This notion is described by introducing the space G/Gx = {aGx | a G G} of 
left cosets of the subgroup Gx in G, where each left coset aGx = {ab\ b G Gx} 
is an orbit of the right translation action of Gx on G. All elements of a given 
coset map x to the same point of M and elements of different cosets necessarily 
map X to different points as one may easily check. One can therefore extend the 
domain of the map / from G to G/Gx when acting on .x, i.e. Fx{aGx) = fa{x) 
defines a map Fx : G/Gx G ■ x which turns out to be a diffeomorphism 
of the left coset space onto the orbit [49]. In particular the dimension of the 
orbit is the difference in dimension of G and Gx- Note further that fixing the 
point in question to be xq, so a general point of the orbit may be represented 
by X = /a(xo) = Fxo{aGxo), then the left action of G on M x ^ fa^^x) = 
/aia(xo) = Fxa{a,iaGxo) corresponds to the left translation aG^o — > cuaG^o on 
the coset space. 

In mathematics any orbit of a transformation group (therefore diffeomorphic 

to G/ H for some subgroup H of G) is called a homogeneous space, since all 
points of the space are equivalent under the transformation group. Here in the 
context of general relativity, a narrower notion of homogeneous space is required 
which incorporates not only the equivalence of the points of the space h\it of the 
geometry as well. A (pseudo-) Riemannian space (M, g) is called homogeneous 
if it is the orbit of an isometry group (invariance group of the metric g), in 
which case the action is said to be transitive. A simply transitive action is one 
in which the isotropy group at every point of the single orbit is trivial (contains 
only the identity ao of G); in this case the orbit and the group are diffeomorphic 
and the left action of G on M corresponds to left translation on G. Using the 
diffeomorphism Fxq : G ^ M for an arbitrary point xq of M to pull back the 
metric from M to G, one therefore obtains a left invariant metric on G. Thus 
a homogeneous (pseudo-) Riemannian space with a simply transitive isometry 
group is equivalent to a left invariant (pseudo-) Riemannian manifold involving 
that group. (For a right action one simply replaces left by right everywhere in 
the above discussion; the choice of left or right actions is a matter of convention.) 

However, not all transformation groups act transitively and the interesting 
question about a given action is how the orbits fit together to fill up the entire 
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manifold. One may introduce an integer- valued fimction dcix) = dim(G • x) on 
M whose value gives the dimension of the orbit to which each point belongs; 
all those orbits of a given dimension form a subspace called a stratum and the 
partitioning of M into the various strata is called a stratification [53] . Note that 
it is easy to verify that if X2 = 021 • xi and a G Gx^, then a2iaa2i G Gx^, i-e. 
Gx2 = 021^(^X1021, so the isotropy subgroups at different points of a given orbit 
are all conjugate (and therefore isomorphic) subgroups of G. Since dcix) = 
dim G — dim Gx, a decrease in the orbit dimension corresponds to an increase 
in the isotropy subgroup dimension. 

Often the action of a Lie group G on a manifold M describes a symmetry, all 
points of a given orbit being equivalent in some sense which depends on the con- 
text, and one is interested in how things change in the directions "orthogonal" 
or "oblique" ("transversal") to the orbits. It is therefore natural to introduce 
the orbit space M/G = {G ■ x \ x G M}; however, due to the varying dimension 
of the orbits, this is not a manifold. For nice enough actions one can usually 
choose a subspace of M (a submanifold with or without boundary or a collection 
of such subspaces which intersects each orbit only once or finitely many times) 
such that its intersection with the "generic" stratum of maximum dimension 
orbits is a submanifold whose tangent space is complementary ("transversal") 
to the orbit tangent space at each intersection point and hence this submanifold 
is a local slice for the action on the generic stratum [53]. This "slice" is very 
helpful in studying objects which are invariant under the group and seem more 
complicated when studied on the entire space M. 

For example, consider the rotations about the 2:-axis of R^. The group is 
SO {2, R) ^ acting as an isometry subgroup of the Euclidean metric on R^, 
the orbits are circles centered on the z-axis and lying in the planes of constant 
z, and the half plane y = 0,x > directly parametrizes the orbit space which is 
a manifold with boundary. On the other hand the full plane y = is a manifold 
intersecting the generic orbits (circles of nonzero radius) twice but having the 
advantage that the projection of all geodesies of the Euclidean metric onto 
this manifold are smooth curves, while those which intersect the z-axis suffer 
reflection at the boundary when projected onto the half plane. It is convenient 
to use the term "slice" to refer to either the plane or the half plane. 

The spatially homogeneous spacetimes or "Bianchi cosmologies" which are 
studied here have a 3-dimensional isometry group G acting simply transitively 
on a 1-parameter family of spacelike hypersurfaces (the orbits) which provides 
a natural slicing of the spacetime. Each orbit is a homogeneous Riemannian 
manifold and therefore isometric to a copy of G equipped with a left invariant 
Riemannian metric, namely the pullback of the induced spatial metric on the 
orbit. Rather than maintaining the distinction between G and each orbit, it is 
simpler to identify the spacetime manifold M with the product manifold Rx G, 
where R is the real line with natural coordinate t which parametrizes the family 
of copies Gi = {(t,x) | x £ G} of G in Rx G, on each of which G acts by 
left translation. However, this still leaves open the question of how these left 
invariant Riemannian manifolds fit together into a spacetime and how the left 
translations on each copy of G in ii x G fit together into a global action of G 
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on M. 

One needs to describe a threading of this natural shcing by a congruence of 
curves in the spacetime (which is nowhere tangent to an orbit) which will be 
identified with the t-lincs in Rx G, the i-coordinatc on Rx G corresponding to 
a given time function for the slicing. Each such identification leads to a different 
global reference system based on the same natural slicing of the spacetime. In 
order for the induced metric on each copy of G to be a left invariant metric, 
the class of threading congruences must be compatible with the action of G on 
the spacetime. One threading congruence and slicing parametrization is picked 
out uniquely by the symmetry, namely the invariant congruence of geodesies 
orthogonal to the orbits, the proper time along these geodesies measured from 
some initial orbit serving to parametrize the family of orbits. Identifying this 
parametrized congruence with the t-lincs oi RxG establishes the so called "syn- 
chronous reference system" [48] adapted to the spatial homogeneity, with the 
action of G on M = R x G being f-independent left translation on each copy 
of G. Any other threading of the slicing may then be viewed as a t-dependent 
diffeomorphism of the family of copies Gt of G in M relative to the synchronous 
threading [47]. Dragging along the t-dependent left invariant spatial metric by 
this diffeomorphism will lead to the t-dcpendent spatial metric in the reference 
system adapted to the new congruence. This will again be a i-dependent left 
invariant metric only if one restricts the spatial diffeomorphism freedom, i.e. 
restricts the allowed class of thrcadings, to be compatible with the group struc- 
ture of G. The compatibility condition is that such diffeomorphisms map the 
space of left invariant tensor fields on G into themselves. These consist of the 
left and right translations and the automorphisms of G, the latter diffeomor- 
phisms being those which preserve the group multiplication and form a finite 
dimensional Lie group Aut{G) = {h G T^{G) | h{a{) h{a2) = h{aia2)} called 
the automorphism group. The translations and automorphisms together form 
a semidirect product Lie group [51] Lq Xg Aut{G) = Rq Aut{G) = 'D{g). 

Before discussing in detail the structure of a spatially homogeneous space- 
time, it is worth understanding first the homogeneous Riemannian manifolds 
from which they are constructed. These are left invariant Riemannian 3-mani- 
folds (G, g), where g is a left invariant Riemannian metric on the 3-dimensional 
Lie group G. On each Lie group there is a natural identification of the tensor 
algebra at any given point, say the identity ao G G, with the algebra of either 
left or right invariant tensor fields. Given a tangent tensor at the identity, one 
can left (right) translate that tensor all over the group using the differential 
of the unique left (right) translation which maps the identity to each point of 
the group, thus obtaining a left (right) invariant tensor field which coincides 
with the original tensor at the identity. In particular, given a basis e = {e-a} 
of the tangent space at the identity and its dual basis {li"} of covectors (satis- 
fying a)"(ef,) = 5°b), one obtains a global left (right) invariant frame e = {eg} 
(e = {ca}) and its dual frame {w"} ({w°}) of left (right) invariant 1-forms on 
the group. The components of a given left (right) invariant tensor field in this 
frame are just the components (namely constants) of the original tensor at the 
identity with respect to the given basis of the tangent space there. For example. 
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a left invariant Riemannian metric may be expressed in the form 

g = gab <Siu)^, gab = g{ea,eb) (2.1) 

where the constant matrix g = (gab) is symmetric and positive-definite. This 

relation in fact establishes a diffeomorphism (for each left invariant frame e) 
from the space of left invariant metrics on G onto the space M of symmetric 
positive-definite matrices of the given dimension. For dimension three, A4 is 
a 6-dimensional submanifold of the space gl{3, i?) of 3 x 3 real matrices whose 
natural basis will be designated by {e^a}, in terms of which a matrix may be 
represented as A = {A°'b) = ^"6e*a- 

Let g and g denote the spaces of respectively left and right invariant vector 
fields on G, each isomorphic as a vector space to the tangent space at the identity 
TGao ^^'^ having corresponding bases e and e arising from some basis e of TGao ■ 
These vector spaces turn out to be closed under the Lie bracket operation and 
are therefore Lie subalgebras of the infinite dimensional Lie algebra X(G) of 
smooth vector fields on G. As a Lie subalgebra of X{G), each generates a finite 
dimensional subgroup of the group 'D{G) of diffcomorphisms of G into itself; 
g {g) generates the action of G on itself by right (loft) translation, with image 
diffeomorphism subgroup Rq{Lq). The Lie algebra g of left invariant vector 
fields on G is referred to as the Lie algebra of the Lie group G. 

A Lie group G is completely determined locally by the structure of its Lie 
algebra g. Given a basis e of g, this structural information is contained in the 
collection of (constant) components of the structure constant tensor defined by 

[ca, fib] = C^^abCc or C ab = i^^^ea, eb]) . (2.2) 

Since e is also a global frame on G with dual frame {w"}, a standard formula 
gives the dual relation 

dw" = -iC^bc'^'' A (2.3) 

Similar formulas hold for e and {Cj""} except for a change in sign of the struc- 
ture constant tensor components, while [ea,eb] = since the diffeomorphism 
subgroups Rg and Lq they generate commute with each other due to the asso- 
ciativity of the group multiplication. The structure constant tensor components 
arc not arbitrary but must be antisymmetric in the lower indices and satisfy a 
quadratic identity imposed by the cyclic Jacobi identity. Let 

e = {c\o I = = qlc%]4 (2.4) 

be the space of possible real structure constant tensor components, a 6-dimen- 
sional space for 3-dimensional Lie algebras. 

Of course one may always choose another basis = A~^^^eb of g leading to 
new structure constant tensor components 

C^bc = A'^dCfgA-^A-^ = jj^{C\e) (2.5) 

which describes the same Lie algebra structure. In fact when the structure 
constant tensor components of two different Lie algebras of the same dimension 
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arc related in this way, the Lie algebras are called isomorphic and represent the 
same abstract Lie algebra. (A simple change of basis leads to bases of the two 
Lie algebras with identical structure constant components.) When C be = C^bc 
so that the components of the structure constant tensor are invariant under 
the linear transformation, then A is the matrix of an automorphism of the Lie 
algebra into itself. In other words the isotropy group of the above left action 
j of the general linear group on C at C"^bc is just the matrix representation of 
the automorphism group Aut{g) of the Lie algebra with respect to the basis e; 
denote this matrix group by Aute{g)- The orbits of the action of the general 
linear group on C correspond to the isomorphism classes of structure constant 
tensors. These isomorphism classes are designated by their Roman numeral 
Bianchi type following the original classification scheme of Bianchi [10] . 

In three dimensions the structure constant tensor is easily decomposed into 
its irreducible parts under the action of the general linear group GL(3, i?), 
greatly simplifying matters. One may dualize the antisymmetric pair of in- 
dices leading to an equivalent second rank contravariant tensor density whose 
antisymmetric part may be represented as the dual of a covector, leading to the 
following decomposition due to Behr [13,14] 

(jab ^ Ica^^gbcd ^ (j{ab) (j[ab] ^ ^ab _^ ^abc^^ 

C\c = C^'^cdbc = cbcdn'"' + afSl:, aj = ^C^f^ (2.6) 
= ufTi^"- = afC^"- = afC^ab ■ 

The Jacobi identity requires that the covector be annihilated by the symmetric 
tensor density. When this covector is nonzero, one may introduce a scalar h by 
the following formula [38] 

dadb = ^heacd^bfgn''^n'^^ ■ 2.7) 

These objects transform under the left action (2.5) of GL{3, R) on C in the 
following way 

ttb = acA~\ = jji^{ab), h = h = jp^{h). 

One may always diagonalize the symmetric component matrix n = (n"'') 
by an orthogonal transformation with matrix O e 0{3,R) . (The eigenvalues 
of n change sign if detO = —1.) The Jacobi identity guarantees that the 
covector may be chosen to lie along the dual of one of the eigenvectors of n, 
thus reducing the components of the structure constant tensor to the following 
"standard diagonal form" 

n = diag(n(^\n(2),n(^)) , af = aS^f (a > 0) ^^^^ 



an 



(3)=o, a2 = /in«n(2) 



Denote the corresponding subspace of C by Cd', this subspace turns out to 
contain all the interesting information. (It is in fact a "slice" for the action of 
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1 




II 
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IV 


1 
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III=VI_i 


1 


-1 
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-1 


VIo 


1 


-1 











VI/i5^0,-l 
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-1 





a 




VIIo 


1 


1 











Vllfc^o 


1 


1 





a 


a" 


VIII 


1 


1 


-1 




















IX 


1 


1 


1 





















Table 1: Canonical values of the standard diagonal form structure constant 
tensor components for each Bianchi type. 

the orthogonal group on C.) If e is the basis of a Lie algebra whose structure 
constant tensor components are in standard diagonal form, then the Lie brackets 
of the basis vectors are given by 

[62,63] = n^^^ei -062, [63,61] =71(^^62 + aei, [61, 62] = 71^^)63 . (2.10) 

Standard diagonal form is preserved by all diagonal matrix transformations 
(provided the third diagonal component is positive when a > 0) and certain 
permutations. Such transformations may be used to further reduce these com- 
ponents to canonical values for each orbit. Only the absolute value of the signa- 
ture of n, the constant h when well defined and the vanishing or nonvanishing 
of the covector with component row vector (a/) are invariant under the general 
linear group. By normalizing the nonzero diagonal values of n to absolute value 
unity, permuting these diagonal values if necessary and changing their overall 
sign using reflection matrices of negative determinant, while normalizing a to 
unity when nonzero and h is undefined, one may arrive at the particular choice 
of canonical values of the structure constant tensor components listed in Table 
1 for each isomorphism class or Bianchi type. The apparently odd choice for 
Bianchi type II reflects a prejudice which tries to associate the third basis vector 
with a preferred basis vector of the Lie algebra (2.10) when possible. (A conflict 
arises for Type IV which does not allow a choice corresponding to the type II 
choice.) Figure 1 represents the space as a 3-plane (class A submanifold) 
and an orthogonal half 3-plane (class B submanifold) in i?'* and indicates the 
canonical points of C corresponding to Table 1. Although is the union of 
two manifolds, it is clearly not a manifold itself. It is convenient to think of 
Cb as stratified by values of the integer pair (rank n, rf^,^ ) , where rf^,^ is the 
reduced orbit dimension, namely the dimension of the intersection of an orbit 
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0, a>q (cLass B) 



Figure 1: The parameter space Cd as a 3-plane (class A) and an orthogonal 
half 3-plane (class B) in with coordinates {n^^\ n^^\ n^^\ a) , showing the 
canonical representatives of each Bianchi type. Of the 8 open octants in the 
class A case, 2 and 6 respectively represent type IX and VIII, while half of 
the 12 open faces bounding these octants represent type VIIq and the other 
half type VIq; the 6 coordinate open half lines represent type II and the origin 
type 1. Similarly in the class B case, half of the 4 open octants are associated 
with each of the 1-parameter family of Bianchi types Vlh and Vllh, a single 
isomorphism class corresponding to a constant value surface of the function 
h = a'^{n^^^n^'^^)~^ . A typical such surface is illustrated in one octant, the angle 
9 given by tan^ = \h/2\^^^; those in the remaining octants are obtained by 
rotation through multiples of 7r/2, h alternating in sign for a given magnitude 
\h\. The 4 vertical open faces bounding these octants all represent type IV and 
the positive a-axis type V, with the a = plane giving the class A limit of each 
type. 

with Cd ; the strata are then labeled by the Roman numerals (excluding III and 
omitting subscripts on VI and VII) of the Bianchi types. Only types VI/i<o and 
VII/i>o represent strata consisting of a family of orbits; the remaining strata are 
themselves orbits. 

The space Cd is very useful in describing the notion of Lie algebra con- 
traction [50]. Consider the effect on Cd of an arbitrary positive-definite di- 
agonal matrix transformation, i.e. an element of the 3-dimensional abelian 
"scale group" Diag{3, R)^ (the identity component of the diagonal subgroup 
Diag{3, R) of GL{3, R) whose Lie algebra diag{3, R) consists of the diagonal el- 
ements of .9/(3, R) ) which represents independent scalings of the standard basis 
vectors of i?"^ or of the basis vectors of any 3-dimensional vector space. Such a 
matrix may be represented in the form 



= diag(e'^\e'^',e'^') €Diag{3,R) 
i8 = diag(/3i,/32,/33) e diag{i,R) , 



+ 



(2.11) 
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Diin(C» n Cd) 



Class A 



Class B Canonical SAut^ig) C 5L(3, R) 



3 



IX • • • VIII 



50(3, R) 



SL{2,1) 



2 







1 




x,5L(2,3) 



T2 X 



SL{3,R) 



S 



A4 



T2 x,SL{2,3) 



Table 2: The "reduced" specialization diagram describing the possible Lie 
algebra contractions and deformations of the Bianchi types and their orbit di- 
mensions restricted to the standard diagonal form subspace Cd of the space 
C of possible structure constant tensors. Not shown are direct paths between 
Bianchi types which may be connected through other Bianchi types by the in- 
direct paths shown. The final column indicates the identity component of the 
canonical special automorphism matrix group of each type, with notation ex- 
plained in appendix B. 

and its left action on Cd via (2.5) is 



The barred components represent another point in the same orbit as long as 
the scale transformation is nonsingular. However, if one takes a singular limit 
a point on the boundary of an orbit can be reached resulting in a change of 
Bianchi type. This is called Lie algebra contraction [50]. 

When a given stratum consists of a family of orbits as is the case for types 
VIh<o and VII;i>o, in order to arrive at a point of the boimdary of a given stra- 
tum not at the boundary of the starting orbit, one must allow motion transversal 
to the orbits; such a motion is called a Lie algebra deformation. For these two 
Bianchi types, changing the parameter h represents a Lie algebra deformation. 
For example, a type IV or V point of Cd can be reached from a type Ylh^o or 
type VII/,^0 point only by a deformation. 

Apart from trivial permutations, each of the canonical points of Cd may un- 
dergo such Lie algebra contractions and/or deformations to arrive at canonical 
points lying in the same stratum or in lower-dimensional strata at the boundary 
of the given stratum. The various possibilities are illustrated in Table 2 follow- 
ing MacCallum [39]. Each class B Bianchi type has a corresponding class A 
limit obtained by the contraction (V, IV) or deformation (Vlh, VII/i) a — > 
shown in Table 2. (Unfortunately the type II components arising from this 



j^(n(«)) = e^ 



a = 



jp{a)=e-l'a. (2.12) 
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contraction of the canonical type IV components differ from the canonical type 
II components by a permutation.) By extending the scale group to the com- 
plex domain, one may perform a rotation in the complex plane which directly 
connects the canonical components of certain Bianchi types. For example, the 
scaling diag(e*^ e^^ 1) with G [0, f ] is a path connecting the canonical com- 
ponents of types VIII and IX, which represent inequivalent real forms of the 
same complex Lie algebra. The types VII/i^o and Vllft^o also require analytic 
continuation of the parameter a as well. Such pairs are connected by horizontal 
dotted lines in Table 2. 

If the real scaling matrix is highly anisotropic, i.e. "nearly singular" , its 
action on a canonical point of Cd may simulate a Lie algebra contraction on 
the space of functions on Cd- the value of a function of the structure constant 
tensor components at the image point will approximately equal its value at the 
nearest point of the boundary. Thus under very anisotropic scalings, functions 
of a given Bianchi type structure constant tensor approach those of a contracted 
type. This is a very useful way of viewing the behavior of highly anisotropic 
Bianchi cosmologies, as will be described below. 

Equations (2.10) represent a Cc-parametrized Lie algebra g^^ . Using a trick 
involving the linear adjoint group, one may realize this Lie algebra as the Lie 
algebra of left invariant vector fields on a C£)-parametrized simply connected Lie 
group G^^ , where the parametrization arises by introducing canonical coordi- 
nates of the second kind with respect to the basis e_ of gi_ whose brackets are 
given by (2.10). These coordinates have range and are global for all Bianchi 
types but type IX where the simply connected group manifold is instead 5^ and 
these coordinates form a local patch centered at the identity. These results are 
summarized in appendix A and illustrate the elegant consequences of the first 
diagonalization referred to in the introduction. 

The result of this long digression is the -parametrized simply connected 
Lie group G^^ which enables one to simultaneously describe all 3-dimensional 
(simply connected) Lie groups. One may next introduce the Cu-parametrized 
left invariant Riemannian 3-manifold {G^^ ) or "homogeneous Riemannian 
3-space" with metric 

gc^ =9a,^:^®^l^ , (2.13) 
where e^^ and } are the explicit C^i -parametrized fields given by formulas 

(A. 9) and the component matrix g = {gab) = gab^^a with determinant g = detg 

lies in the 6-dimensional space A4 C GL{3, R) of component matrices of positive- 
definite inner products on R^. This space, through (2.13), parametrizes the 
space Ml{G^^) of left invariant metrics on the Lie group G^^. Since it is a 
bit awkward, the subscript C^i will usually be omitted in what follows. Before 
moving on to the C^-parametrized spatially homogeneous spacetime, it pays 
to examine the curvature of the Cj^-parametrized homogeneous Riemannian 3- 
space and the isometry classes of the space of such Riemannian manifolds. The 
latter question leads to the second diagonalization mentioned in the introduc- 
tion. 

The isometry classes of A4l{G) are its intersections with the orbits of the 
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difFeomorphism group 'D{G) on the space of all smooth Riemannian metrics on 
G. Consider instead the largest subgroup of 'D{G) which acts on the space 
M.l{G), i.e. which maps all left invariant metrics into left invariant metrics 
under the dragging along action. This is possible only if it maps the Lie algebra 
g into itself and hence the space of all left invariant tensor fields into itself under 
dragging along. The orbits of this group on M.l{G) should correspond to the 
isometry classes of left invariant metrics; it is assumed that they do. 

The "symmetry compatible subgroup" of 2?(G) having this property, already 
designated by T>{g) earlier in this section, is the semidirect product Lie group 
of translations and automorphisms of G 

V{g) = Lg Xs Aut{G) = Rq Aut{G) . (2.14) 

The equality of the two semidirect products is connected with the adjoint group 

ADg C Aut{G) of G, also called the group of inner automorphisms of G. In 
addition to the effective left action of any Lie group G on itself by left translation 
L and inverse right translation which are commuting actions due to the 
associativity of the group multiplication, G may act on itself on the left by 
inner automorphism: ADa = o = R~^ o La- The image group ADg is a 
homomorphic subgroup of the automorphism group of G but is not necessarily 
isomorphic to G. They are isomorphic and the adjoint action effective when 
ADaQ is the only transformation acting as the identity, i.e. when the center 
C{G) = {a G G \ ADa = Id} of G is trivial (contains only the identity ao). 
Returning to (2.14), the fact that AD~^ o a e Aut{G) if a G Aut{G) together 
with the identity La o a = Ra o [ADa^ o a) explains the equality of the two 
semidirect products. 

This discussion may be repeated at the Lie algebra level for the Lie algebra 
X{g) which generates I? (5), a semidirect sum Lie subalgebra of the Lie algebra 
A'(G) of smooth vector fields on G 

A" (.9) = g ®s aut{G) =g®s aut{G) . (2.15) 

Here aut{G) generates Aut{G) and ad(G) =g-g= {X"(ea - §«) | (X") e R^} 
generates the adjoint group. 

When the group 'D{g) acts on g by dragging along, Lg has no action by 
definition, while Rg and ADg have the same action, inducing inner automor- 
phisms of the Lie algebra: R'^X = ADaX = Ad{a)X G g for X G g. The 
image subgroup Ad{G) of the general linear group of g is called the linear ad- 
joint group and when G is connected coincides with the group IAut{g) C Aut{g) 
of inner automorphisms of g ; its matrix representation Ade{g) with respect to 
a basis e of g is exploited in appendix A. Similarly by dragging along, Aut{G) 
induces the action on g of the full group Aut{g) of automorphisms of g when G 
is simply connected as is assumed here. Left invariant tensor fields undergo the 
transformation associated with the corresponding tensor representation of this 
group. Similarly when X{g) acts on g by Lie derivation (define ad(^)X = £^X 
for ^ e X{g) and X G g and let ade(C) be the matrix of ad(^) with respect 
to the basis e oi g), g has no effect, while g and ad(G) have the same action. 
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inducing inner derivations of g (the image Lie subalgebra ad(g) C aut{g)), while 
aut{g) induces the action of the full Lie algebra of derivations der{g) = aut{g) 
of g, with matrix representation aute{g). The Lie algebra aut{g) generates the 
automorphisms of g . 

Thus when 'D{g) acts on the left invariant metric (2.13) by dragging along, 
the component matrix g undergoes the appropriate transformation law associ- 
ated with the matrix automorphism group. If A G Aute(<)) is an induced matrix 
automorphism of g , this transformation law is 

geM ^ /A(g) = A-i^gA-^ (2.16) 

The orbits of this action of Aut^(g) on A4 through the correspondence (2.13) 
represent the isometry classes of left invariant metrics. However, just as the 
space C could be reduced to its essential structure by diagonalization, here too 
the "offdiagonal" metric matrix variables are superfluous and all the essential 
information is carried by the diagonal submanifold Md of M, assuming that e 
is a basis of g whose structure constant tensor components belong to Cd- This 
submanifold A4 d, like the space Cd, is also a "slice" for the natural action of the 
orthogonal group on the full space. (The reduction ofA^xCtoA^uxC^isa 
consequence of the well known fact that one can always simultaneously diagonal- 
izc two real symmetric matrices by an orthogonal transformation.) In fact A4u 
is a "slice" for the action (2.16) of any 3-dimensional subgroup G C GL{2>,R) 
whose matrix Lie algebra g has a basis {/Ca} with the following property: for each 
cyclic permutation (a, 5, c) of (1, 2, 3), the matrix /£„ belongs to span {e^c, j,}, 
where {e"(,} is the natural basis of gl{3,R) already introduced above. In ap- 
pendix B, the matrix automorphism group is described for the Co-parametrized 
Lie algebra g . It always contains such a subgroup G, which may be used to 
map a general point of M to the diagonal submanifold Md- It therefore suf- 
fices to consider those automorphisms which map Md into itself to determine 
the isometry classes. Since Aljj is mapped into itself by all permutations and 
diagonal transformations, it suffices to consider elements of Aute{g) of this type. 

M D consists of all diagonal matrices with positive entries and clearly coin- 
cides with the scale group Diag{3, R)'^ as a submanifold of GL{3, R) . They are 
best identified, however, in terms of the simply transitive action of Diag{3, i?)"*" 
on Md using the identity matrix 1 e A^d as a reference point, as described at 
the beginning of this section. The abelian group Diag{3, R)'^ is most naturally 
parametrized by its Lie algebra diag{3,R) which in turn parametrizes Md 

/3 = diag(/3\/32,/?3) , = diag{e^\el'\e^') , g' = /;/(!) = e'^ 

e e e 

diag{3,R) D Diag{3,R)+ = Md ■ 

(2.17) 

The prime on g' serves as a reminder of its diagonality, while the right action 
is used to conform with convention. Through (2.17) the single matrix 

simultaneously represents three different diagonal matrices. The special scale 
group SDiag{3, R)^ = Diag{3, i?)"*" fl SL{3, R) with Lie algebra sdiag{3, R) = 
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{0 G diag{3, i?) | Tr j8 = 0} is related in a similar way to the unimodular 
submanifold M.d = A^d H SL{3,R); the following notation proves convenient 

^=ilng' = /30l+^, f3° = ^lng' = ^Tr0, 

/?« = i In^;^ = + ^« , = /?« - = /3« - = i HgUg'bb) ■ 

(2.18) 

Misner [15] introduced a basis oi diag{2>, R) and sdiag{3, R) which is orthonor- 
mal with respect to the inner product ((, )) = | (> on 5^"^' -^)> where the 
DeWitt inner product and trace inner product on gl{3,R) are defined by 

(A,B)^^ =Tr AB-Tr ATrB , (A,B)=TrAB, A,B Ggl{3,R) 
{{,)) = l{,) ■ 

(2.19) 

This basis and the corresponding parametrization of diag{3, R) and sdiag{3, R) 
are given by 

= l3^eA = /3°eo + (3+e+ + /J^e^ 
{eo,e+,e_} = {1, diag(l, 1, -2), V3diag(l, -1, 0)} 

(riAB) = i{{eA,eB))^J = diag(-l, 1, 1) = (ry^^) 

= /3+ + \/3/3- , P^=P+- V3p- , = -2/3+ 
/323 = 3^+_^^- ^ = -3/3+-^/3/3- , = 2\/3/3- . 

(2.20) 

These may be generalized by the definitions 

/3^ = -i/3", /3- = (4V3)-W% ^2 21) 

^ = /3a ea+ + 0~ea- (no sum on a) , 

with /3+ = pf and the others obtained by cyclic permutation of indices. For 
each cyclic permutation (a, 6, c) of (1,2,3), the Taub submanifold A4T{a) = 
{g' G Alz) I (7^(, = and its unimodular submanifold A^T(a) C A4d may 
be equivalently defined by 0^" = or /3~ =0. They intersect at the isotropic 
submanifold M.i and A4i = {1} respectively, for which /3+ = /3~ =0. Since 

g' = e^^, translation along represents a conformal rescaling of the metric 
(2.13) under which all curvatures scale by a factor e'^^ where q is an appropriate 
dimension. Thus the nontrivial information about curvature is associated with 
the conformal submanifold A4d (namely /3° = 0). 

The plane, i.e. sdiag{3,R) ~ SDiag{3, R) ~ A1_d is illustrated in 

Figure 2, indicating the Taub and isotropic submanifolds and each of the pairs 
of coordinate axes associated with the three coordinate systems {/3+, Each 
pair of coordinate vectors, namely {ea+, ea_}, is orthonormal with respect to the 
inner product ((, )). Interpreting the /3+/3^ plane as A4d, a cyclic permutation 
of the basis eoi g leads through (2.13) and (2.16) to a rotation of the coordinate 
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Figure 2: The /?+/?- plane: sdiag{3,R) ~ SDiag{S,R)+ Md- Shown are 
the three pairs of orthogonal axes related by rotations of angle 27r/3. The 
/3+-axis is the projection of the Taub submanifold A^T(a) to Md, while the 
origin represents the projection of the isotropic submanifold J^j. Lines parallel 
to the sides of the triangle are associated with constant values of the three 
diagonal components of the conformal metric matrix while lines parallel to the 
bisectors of the vertex angles are associated with constant values of their ratios. 
Reflections and permutations of the basis e act on M.d as the symmetry group 
of the equilateral triangle of this figure. 

axes by ±27r/3, while a transposition of two basis elements, say Ch and Cc, leads 
to a reflection about the Taub submanifold A^T(a)! where (a,b,c) is a cyclic 
permutation of (1,2,3) . The action of the discrete group of permutations on 
M D thus coincides with the symmetry group of the equilateral triangle shown 
in Figure 2, wlios(^ sides are parallel to the constant value lines of the 
coordinatcs and whose orthogonal bisectors are the Taub submanifolds. (Note 
that rotations about one of the frame vectors by ±7r/2 have the same effect 
on the /3^/3~ plane as transpositions.) The translations of the /3^/?~ plane 
correspond to the action on Md of the special scale group. 

The Lie algebra contractions of the space Cd arising from singular limits 
of the action of the scale group on this space are now easily described. Let 
{/? = sb I s G (— oo, 0]} be a ray from the origin of diag{3, R) parametrized by s 
and extend all of the subscript and superscript notation of (2.17)-(2.21) to the 
constant diagonal matrix b ; then (2.12) becomes 

j-^(n(«), a) = {e^^'°+''i^n^'^\e^^''"-''t)a) . (2.22) 
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If fe" 7^ 0, then one might as wcU set 6° = 1. In order that this have a finite hmit 
as s = /?° ^ — oo leading to a singular scale transformation (dete^ = e^^ 0) 
which therefore induces a Lie algebra contraction, the following inequalities 
must be satisfied when the corresponding structure constant tensor component 
is nonvanishing 

n(«)^0: b+>~; a ^ : 6+ < i . (2.23) 

These inequalities are illustrated in Figure 3. The label of a given dashed line 
indicates the structure constant tensor component which remains fixed under 
the scaling (2.23) associated with the points of the line, all points to the origin 
side of the line leading to a limit where that structure constant tensor component 
goes to zero and all points to the other side not leading to a finite limit if that 
component is nonzero. Thus all points in the interior of the triangle 123 lead to 
the abelian limit (0, 0, 0, 0), vertex 1: (0, n^^), n^^), 0), vertex 2: {n^^\0,n^^\0), 
vertex 3: (n(^\ 0, a), open side 23: (n^^), 0, 0, 0), open side 31: (0, n(=^), 0, 0) 
and open side 12: (0, 0, n^^\0). 

On the other hand when 6" = 0, at least one of the components of must go 
to infinity as s ^ — oo so a finite limit can result only if one of the components 
of n is zero. If (a, b, c) is a cyclic permutation of (1,2,3) and n'"-* = 0, then the 
limit s — > — oo of (2.22) will be finite only if 6^ > and fe+ > 0, which is the 
sector of the plane between the positive 6+ and 6+ axes, the limit being (0,0,0,0) 
between the axes, but with j p{n'^''\_rS'^\a) having the limit (n(''',0, 0) on the 

6^-axis and (0,n(^\0) on the fe+-axis. This class of Lie algebra contractions 
might be called "pure anisotropy" contractions. 

Returning now to the question of isomorphism classes within the diagonal 
submanifold M. d , namely the orbits of the action of the diagonal and permuta- 
tion automorphisms on A4d, one has four different cases corresponding to the 
four categories of Table 2. Modulo discrete automorphisms, the diagonal auto- 
morphisms of the Co-parametrized Lie algebra are described for each of these 
categories in Appendix B. For the first category only permutation and reflection 
automorphisms act on A^^i so A4d itself locally parametrizes the space of au- 
tomorphism group orbits on M . (For the canonical type IX case the six sectors 
into which the three ,5+ axes divide the (3~^(3~ plane are all isometric for a given 
value of but in the canonical type VIII case only reflection about A4t{3) 
connects isometric points.) For the second category of Table 2, there exists a 
diagonal automorphism subgroup generated by the matrix I^^^ = diag(l, 1,0) 
when n'^^^ = 0, leading to translations along /3+ = 0^ in the f3^(3~ plane, so 
P~ = (5^ together with locally parametrize the orbit space. The result for the 
other components (n^ ') / 0) of Cd belonging to this category may be obtained 
by cyclic permutation. For the third category automorphisms induce transla- 
tions along both /3+ and (3~ so all points of the /?"*"/?" plane are equivalent and 

alone parametrizes the orbit space, while in the abelian case M. consists of 
a single orbit. 

For the upper two categories of Table 2, the generic points oiM.D belong to 
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Figure 3: The parameter space for the family of Lie algebra contractions of the 
space Cd- 

an orbit on M. having three dimensions transversal to A^d and at most a dis- 
crete isotropy group. However, on certain submanifolds oi Ado and for certain 
Bianchi types the full orbit dimension decreases and only one or no directions 
remain transversal to A^^. At these points the isotropy group has dimension 
greater than zero corresponding to additional symmetries of the metric (2.13). 
This occurs at the Taub submanifold M.T{a) when n*^^^ = n'^'^^ ((a, 6, c) is a cyclic 
permutation of (1,2,3)) corresponding to local rotational symmetry (necessar- 
ily the index a = 3 in the class B case) and at the isotropic submanifold M.i 
when n'-^-' = n'-^-' = n^'^^ corresponding to isotropy. For the lower two categories 
the isotropy group has generic dimension greater than zero so there are always 
additional symmetries. However, the Taub submanifolds are still relevant to 
spacetime symmetries. The choice of canonical components for Bianchi type 11 
was made so that M.t(z) is associated with additional spacetime symmetry for 
all canonical points oiCo- This is discussed in greater detail elsewhere [43]. For 
the noncanonical points of Cd the submanifolds relevant to additional symmetry 
change as described below. 

The discussion of isometry classes and additional symmetries is not just an 
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interesting aside, but is important for appreciating the symmetries of the scalar 
curvature R of the metric (2.13). This function on At x is a scalar under 
a change of basis e of g and hence is invariant under the action of the matrix 
automorphism group on A4 alone, having a constant value on each orbit. Using 
standard formulas one may easily evaluate the components of the connection of 
the metric (2.13), raising and lowering all indices with the component matrices 
g and 

Ve^et = r^abCc , = \C\t, + . (2.24) 

Introducing the unit alternating (pseudo-)tensor rjabc = d^^abc and the two 
matrices m = (7~2ng and A = acrj^le^a, the Ricci tensor and scalar curvature 
of this connection are then found to have the following expressions 

R = R^be^ = - mTr m - l(Tr - iTr^m + 2aca^) + [m,A] 

i? = TrR=-(Trm2-iTrV)-6aca% ^^"^^^ 

where the conventions of Misner, Thorne and Wheeler [18] are followed for 
curvature tensor definitions. Assuming as always that C^bc G Cd, these are 
Cu-parametrized functions on M. Note that for g e Md the Ricci tensor 
component matrix is diagonal except for the last term which contributes a 12 
(and 21) component in the class B case. In the class A case e is then an 
orthogonal frame of Ricci eigenvectors, while linear combinations of ei and 62 
must be taken to obtain such a frame in the class B case, leading to structure 
constant tensor components not belonging to C/j. 

One may also introduce a potential function and several 1-forms on A4 which 
may be interpreted as force fields 

Ug = -g^R 
G = -g'^C'Hgab = -g*(i?"' - hRg'''')dgab 
Q = Q^'dgab = 2gi {a'C^\ - 2a''a'')dgab (2.26) 

G = -dUc + Q . 

The scalar curvature potential function Ug serves as a potential for the Einstein 
force field G in the class A case where Q vanishes, but in the class B case the 
Einstein force field has a nonpotential component Q which generically satisfies 
Q ^ ^ dQ. As a C_D-parametrized function on A4d, the scalar curvature 
potential is given explicitly by 

UG = e^°{V* +6a^e^f'^) 

V* = 5EL1 (n("))2e-8^'^ - [n(2)n(3)e4/3^ +„(3)„(i)e4/3^ +n(i)n(2) 6^/33+] 
= 2e4'3"[i(n(i)e2v^/'" - n^''^ e-^^^~ )]'' 

_2„(3)e-2/3+[l(^(l)g2v^/3- +^(2)g-2V^/3-)]+ l(^(3))2g-8/3+ 

(2.27) 
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The values of the curvature function V* (a scalar density of weight | which may 

be expressed in the form V* = g~^QabcdP'°'^''^'''^ using (2.40)) at the canonical 
points of Cd are 



IX/VIII 
VII 
II 

I,V 



2"= 



2e4^'+sinh22V3/?- T 2e-2/3+cosh 2^3/3- , 
2e^''''sinh22\/3/3- VI: 2e^l^^ cos\v^2V^p 

\e-^^* IV: 
. 



(2.28) 



Suppose y is a function on x C which is a scalar under a change of basis 
e of g and therefore satisfies y(g, C^c) = ^(/a(s)' Ja(^"''c)) equivalently 
^(/^^(g),C'"fec) = yiE^JAiC^bc))- Focussing now on Md x Cd and letting 
jj^ be one of the Lie algebra contractions (2.22)-(2.23), one sees that the effect 
on the function Y of such a contraction is equivalent to an infinite translation 
oi M.D diag{3,R) ^ R^. The value of the function for the original point 
C°6c therefore approaches its value for the contracted points as one approaches 
infinity in Md diag{3,R) ~ R^ in the negative /?° direction. e~^^ V* is 
such a scalar function and hence as one approaches infinity in the plane 
along the directions parametrized in Figure 3, the density V* approaches a 
rescaled version of the potential of the corresponding contracted points of Co- 
Furthermore the difference between V* and its contracted value at the same 
point of the P~^l3~ plane as one gets far from the origin becomes very small 
compared to the value itself. 

This can be seen in the diagrams of Figure 4 which show suggestive contours 
of the potential V* for canonical points of Cd ■ Contours of the same five function 
values are shown for each type, together with two additional closed contours for 
the type IX region where V* is negative. As one proceeds along any of the 
positive /?+ axes in the type IX case or the /?+ axis in the type VIII case, 
the potential quickly approaches that of a permutation of the canonical type 
VII potential, while along the positive and axes the type VIII potential 
approaches a permutation of the type VI potential. For directions in between 
these three positive axes the potential quickly approaches that of a permutation 
of the type II potential. Similarly the type VII and VI potentials approach 
permutations of the type II potential for directions between the positive 
and axes and the positive /?+ and axes, but type I in the remaining 
sector. Finally the type II potential approaches type I along all directions in 
the positive /?+ half plane. Note further that for all the nonsemisimple types 
the potential V* simply scales under translation along /3+, so the contours are 
simply translates of each other, a consequence of the existence of the additional 
diagonal automorphism generated by the matrix I^^-* = diag(l, 1, 0), except for 
type II where the matrix is instead diag(l, 1, 2). The reflection and permutation 
symmetries of all of the potentials reflect the existence of discrete symmetries. 
The reflection symmetry about the axis for all types but IV is connected 
with a discrete automorphism whose existence motivated the choice of canonical 
type II components. 

The dashed straight lines in Figure 4 indicate "channels" of width 1 outside 



23 




24 



of which the contours arc essentially the same as the corresponding asymptotic 
Bianchi type II potential (the difference becoming exponentially small with dis- 
tance from the origin). These channels themselves are either open or closed; 
outside of the dashed circles of diameter 1 in the type IX and type VIII figures, 
the open and closed channels are essentially the same of the corresponding type 
VII and VI channels respectively, where these latter channels are rotated by 
±27r/3 for comparison. (Corresponding contours are separated by a distance of 
less that .01 as one exits the circle and the difference decreases exponentially 
with distance from the origin.) Figure 5 shows the type VII and VI contours 
of the same function value together with the asymptotic type II contours of 
the same function value, showing that the deviation of the open and closed 
channel contours from the type II asymptotes only becomes important within 
the channel itself. The only region of the type IX and VIII potentials which 
is essentially different from those of the remaining Bianchi types (arising from 
them by contraction) is the interior of the dashed circle which occurs at the 
intersections of the three channels; similarly only the channels of the type VII 
and VI potentials are different from the potentials of the contracted types I and 
II. 

The potentials for the noncanonical points of C/j are obtained merely by 
translating the origin of coordinates in the plane and rescaling the po- 

tential. Let C^bciccm) {n{can),aa{can)) be the canonical point of Cd in the 
same orbit as C^tc ^ (n, Ua) € Co- One can then define the nonsingular matrix 

J = 7-^6^ G diag(3, R) by 

C\c = jr\C\,{can)) , F e Diag{2,,R) , = diag(|rM, jr^sj, jr^g]) . 

(2.29) 

Then one has the identity 

y*(n, /?±) = e2TV*(n(con), /J± - 7±) , (2.30) 

showing that V* is obtained from its canonical value by an active translation of 
the plane by (7~'~,7~) and a rescaling to its function values, leaving the 

shape of its contours unchanged. A Lie algebra contraction then corresponds 

to an infinite translation. All of the Lie algebra c;ontractions parametrized by 
Figure 3 act on the potentials of Figure 4 to reduce the more complicated ones 
to successively simpler ones. 

Figure 4: Suggestive contours of the potential function V* on the P~^/3~ plane 
are shown for the canonical points of Co (not shown is the type IV case which is 
related to type II by an active rotation of 47r/3) and values of V* are plotted for 
the Taub submanifold Mt(3), namely the /J^-axis. Unit distances are marked 
on the coordinate axes. Arrows indicate the direction of the associated force 
field, pointing toward directions along which the potential decreases. Note that 
V* is nonnegative except in the type IX case where the contours are closed for 
y* < and open for V* > 0, the minimum value -1.5 occurring at the origin or 
isotropic submanifold A4i. 
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Figure 5: Open and closed channel contours and their asymptotes. 



Having exhausted the essential points regarding homogeneous 3-spaces, the 
discussion may proceed to the spacctimc level, introducing the C^i-paramctrized 
spacetime {M^^ g^^ ) . The spacetime manifold is M = i? x G^^ , with the 
natural coordinate t on the real line R parametrizing the 1-parameter family 
of orbits of the natural left action of Gcd on Mq^ , namely ^-independent left 
translation of each copy of Gco in the product manifold. The copies of R in 
Mcj3 are the t-lines which are interpreted as the normal geodesies to the family 
of orbits, with t coinciding with the proper time along these geodesies. The 
spacetime metric may therefore be written in the following form referred to as 
synchronous gauge (zero shift and unit lapse) 

^g^^ = -dt ®dt+ gab{t) ® ojI^ . (2.31) 

The vector field = eo = d/dt is the unit normal to the slicing of M by spa- 
tially homogeneous hypersurfaces. Proper time derivatives will be denoted by a 
small circle ° . A reparametrization of the time t — > t{t) may be accomplished 
by introducing a nontrivial spatially homogeneous lapse function dt = N(t)dt; 
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barred time derivatives will be denoted by a dot ' , so one has the relation 
a = Na for the time derivatives of a function a only of time. Eq.(2.31) with 
nontrivial lapse but zero shift will be referred to as almost synchronous gauge. 

In synchronous gauge, or almost synchronous gauge as long as the lapse 
function is an explicit function of gab and gab or other known quantities, the 
spacetime metric is completely determined by the parametrized curve g(i) in A^. 
Acting on the curve by a parametrized curve A(t) in the matrix automorphism 
group 

E(t)^m = fA(J)(s{t)) (2.32) 

is equivalent to the introduction of a shift vector N(t) belonging to X{g^^ ) and 
satisfying the matrix equation 

ade(iV(t)) = (c^-(i:^(j)e,)) = A(t)A-i(7) , (2.33) 

and of a new spatial frame and off-hypersurface frame vector 

ea = A-\{t)et, , eo=N{t)e^+N{t) (2.34) 

leading to the expression for the metric in a general symmetry compatible gauge 

"g = -N(tfdi (E)di + gab{t){uj'' + N^cR) ® (w'' + N^di) . (2.35) 

Here IV tU°(iV) are the barred shift components, while eq.(2.33) follows 
from the comoving condition [eo,ea] = 0. In this gauge the spacetime metric is 
determined by the parametrized curve g{t) in M together with the parametrized 
curve A(i)A(i)~^ in the Lie algebra of the matrix automorphism group, again 
provided the lapse is known. 

In synchronous or almost synchronous gauge, the extrinsic curvature is sim- 
ple and its matrix of mixed components 

K = {K\) = -ig-ig = -(27V)-ig-ig (2.36) 

represents a matrix-valued function on the gravitational velocity phase space 
TAi , namely the tangent bundle of the gravitational configuration space Ai. 
Here {gab, gab} are the natural "coordinates" on M lifted from the "coordinates" 
{gab} on A4. The ADM gravitational Lagrangian density is a Langrangian 
function Lq on the velocity phase space 

Lg = N{T - Ug) 

(2 37) 

NT = NQ'^^'^'^KabK^d = m-'g''>"'''gabgcd = N g'^ (K,K)^^ . 

The definition of the momentum canonically conjugate to g is simply the as- 
sociated Legendre transformation between the velocity and momentum phase 
spaces 

tt"' = dlG/dgab , ir = (7r%) = -5^ (K - (Tr K)l) . (2.38) 
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The gravitational phase space is just the cotangent bundle T*A4 on which 
{dab,'^"'''} are the natural "coordinates" lifted from the "coordinates" {gab} 
on M. Note that different choices of lapse change the Legendre map. The 
Hamiltonian function on momentum phase space which is associated with Lq 
is defined in the usual way 

Hg = TT'^'gab -Lg = N{T + Ug) = NHg 

(2.39) 

= S-blX'^"" + UG = g^ (Tr - iTrV) + Ug = 2g2 ^ . 

The scalar density Hg is the gravitational super-Hamiltonian. Both Lg and 
Hg are Cij-parametrized functions due to the potential Ug ■ The kinetic energy 
has no parameter dependence and is just a rescaling of the square of the DeWitt 
norm of the velocity vector of the system, where the DeWitt metric ^ on is 
given by [44] 

g = g-'-'^dgab dg,d , , 

(2 40) 

Qabca ^ (^a(c^<i)6 _ ^ab^cd^ ^ ^ ^^^^^^^^^ _ 1^^^^^^) _ 

The kinetic energy T generates the "free dynamics" (the vacuum type I case: 
C°'bc = 0) whose solutions are just the geodesies of the DeWitt metric which 
are afiinely parametrized by the time t in synchronous gauge. A nontrivial 
lapse function N which is an explicit function on M. corresponds to conformally 
rescaling the DeWitt metric g N~^g so that t is an affine parameter with 
respect to the rescaled metric. However, only the null geodesies are relevant to 
the free dynamics due to the free super-Hamiltonian constraint T = which 
requires the tangent vector to the geodesic to be a null vector with respect to the 
DeWitt metric. A general lapse function leads to an arbitrary parametrization 
of these null geodesies. 

The general linear group GL{3, R) acting on A4 through (2.16) is a group of 
homothetic motions of {^4, Q) and the special linear group SL{'3, R) is the iden- 
tity component of its isometry subgroup. For each B G gl{'i, R), the correspond- 
ing homothetic Killing vector field (simply Killing vector field if B e sl{3, R)) 
is given by 

e(B) = -g,^aB^)d/dgab , (2.41) 

where the vector fields d/dgab on A4 are defined by dgcd{d / dgab) = '^'j^'^^d) • 
The Lie bracket of two such fields satisfies [^(A),^(B)] = -^([A,B]) . The 
corresponding generator of the lifted canonical action on the momentum phase 
space and the Poisson brackets of two such generators are given by 

P(B) = -2Tr Btt , {P(A), P(B)} = P([A,B]) , (2.42) 

where the only nonvanishing Poisson brackets of the "coordinates" (5^06) Tr"*") 
are defined by {gab, tt'^'^} = ^'^^a^'^b) ■ The corresponding function on the velocity 
phase space when = 1 (and dot becomes circle) 

P(B) = e-i"n(B)a65cd (2.43) 
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is just the inner product of the corresponding homothetic Kilhng vector field 
with the velocity of the system. However, only the natural action of SL{3,R) 
on the gravitational velocity and momentum phase spaces is a symmetry of the 
kinetic energy for arbitrary lapse and so in general only its canonical generators 
are conserved by the (unconstrained) free dynamics. Null geodesies on the other 
hand conserve P{1) as well, so the free super-Hamiltonian constraint restores 
GL{3, R) as a symmetry group of the allowed solutions. 

In Misner's supertime gauge [15-18] N = I2g2 , the restriction of the metric 
N~^Q to A^D is just the flat Lorentz metric induced on M.d by its identification 
(2.17) with the inner product space {diag{3,R),{{,)) in terms of which 
{f}^} are inertial coordinates 

N-^g I = TjABrf/?^ ® d^"" 

'^D (-2.44) 

NT = IvabP^P"" = Iv^'^PAPB , PA = VabP'' . 

Here {(3"^,$"^} and {i3'^,pa} are the natural lifted coordinates on TM.D and 
T*M.D- The null geodesies on M.d are just null lines in diag{3,R) 

P^{t) = rj^^pst + (3^(0) , PA = = V^^'PAPB 
^{t) = e'^C*\ ^(i) = r?^VeB = b. ^^""^^^ 

In fact since GL{3,R) maps null geodesies into null geodesies (and A^d is a 
totally geodesic submanifold of A^), a general null geodesic is of the form 

g(t) - /a (g'©) . A e GL{3,R) . (2.46) 

The Kasner exponents (pi,P2,P3) are defined to be the eigenvalues of the 
matrix k = (Tr K)~^K for a null geodesic ( (k,k)^ = (k,k) — 1 = 0), therefore 
satisfying pi + P2 + Ps = P\ + P^ + P^ = 1- For the diagonal null geodesies 
(2.45) the extrinsic curvature matrix is K(I)' = —^g~^(t)fi(t) so one has 

diag(pi,p2,P3) = (Tr b)"^b 

= (-3po)"^diag(-po +P+ + \/3p-, -po +P+ - V3p-,-po - 2p+) 
= |l ~ gsgnpo diag(p+ + ^/3p-,p+ - \/3p_, -2p+) , 

(2.47) 

where the unit vector (p+,p_) = (p^ + p^)~2 [p_^_,p_) coincides with 
\Po\~^ip+jP-) for these null geodesies. This gives the Kasner exponents as 
functions on the unit circle in the p+p- plane, representing the iS^ -parametrized 
family of null directions in the diagonal cotangent spaces. 

A very useful parametrization of this circle was given by Lifshitz and Kha- 
latnikov [54] 

(Pi(w),P2(m),P3(m)) = (u^ + M+l)"^(-u, 1 + u,m(1 + m)) 

{p+,p-) = -{u^ + u + ly^u"" + u - ^,V3[u + ^]) . 
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Figure 6: The Lifshitz-Khalatnikov parametrization of the Kasner exponents 
and of the unit circle in the p+P- plane. 

The Kasner exponents as functions of the variable u are shown in Figure 6, 
together with the correspondence with the unit circle in the p+p- plane estab- 
lished by the second equality, which may be inverted to give u± = — ^(l=F[3(l — 
p+)/(l + p+)]2 J where the upper (lower) sign applies in the upper (lower) half 
circle. Each of the six sectors into which this circle is divided by the pi+, P2+, 
and P3+ axes (where the coordinates {pi+,pi-} and {p2+,P2-} are obtained 
from by rotations by 27r/3 and — 27r/3 respectively and {p3+,P3-} = 

{p+,P-}) represents a different ordering of the same interval of eigenvalues. The 
intersections of these axes with the circle represent the three permutations of 
the two inequivalent Taublike case null geodesies whose associated type I space- 
time metrics (2.31) are locally rotationally symmetric. If (a,6, c) is a cyclic 
permutation of (1,2,3) then an interchange of the basis vectors et and Cc leads 
to reflection across the axis in A4d and hence reflection about the pa+ axis 
in each cotangent space, under which the unit circle is invariant. In terms of 
the variable u parametrizing this circle, the reflections across the pi+, P2+ and 
P3+ axes respectively correspond to the discrete transformations P23{u) — u^^, 
Psiiu) = —u{u + 1)^^ and Pi2{u) = —{u + 1) , representing transpositions of 
the basis vectors {ea}.*-^^-* The transformations corresponding to cyclic permu- 
tations of these basis vectors leading to rotations of the (3^ and p± planes by 
±27r/3 may be obtained by combining two transpositions. The index permuta- 
tion (1,2,3) (2,3,1) is represented by P23i(u) = P23oP3i{u) = -u~^{l+u), 
corresponding to a positive rotation by 27r/3 while (1,2,3) (3, 1,2) is repre- 
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sented by Pauiu) = —(1 +u) ^ , corresponding to a negative rotation by 2n/3 
[15]. 

The diagonal geodesies are characterized by the vanishing of the moment 

function (2.42) for any offdiagonal matrix B since both K and w arc diagonal 
for these geodesies. Such values for the moment function will be referred to as 
angular momentum. The constant transformation fj^ applied to the geodesic 
(2.46) "transforms away" its constant nonzero angular momentum. For the 
spacetime metric (2.31) corresponding to (2.46), the new spatial frame e' ob- 
tained from e by this transformation 

< = A-\e, (2.49) 

is an orthogonal frame of eigenvectors of the mixed extrinsic curvature tensor 
(both g' and K' are diagonal) . Such a spatial frame is called a Kasner frame and 
its elements arc; called Kasner axes [55]. Note, however, that unless A G Aute{g), 
the new structure constant tensor components will differ from the old ones and 
will in general not belong to Cd- In the analogy with the central force problem, 
changing the spatial frame by a constant linear transformation corresponds to 
changing the origin with respect to which the orbital angular momentum is 
defined. 

In the supertime gauge GL(3, R) is the isometry group of the rescaled De- 
Witt metric so all of its generators are conserved by the (unconstrained) free 
dynamics. However, invariance of the rescaled curvature potential NUq un- 
der the full automorphism group Aute{g) requires an incompatible time gauge 
N (X g~'^, so SAute{g) is the largest simultaneous symmetry group of both the 
kinetic and potential energies. 

For a nonvacuum spatially homogeneous spacetime with spatially homoge- 
neous energy-momentum components T"/3 in the synchronous gauge, matter 
variables may usually be chosen so that the matter super-Hamiltonian Hm = 
—2kg2T ^ acts as a potential function for the matter driving force that ap- 
pears in the evolution equations. If cIm is the exterior derivative on M., then 
this potential must satisfy 

T* = kg^T'^^'dgab = -dMUu ■ (2.50) 

If this is not possible, one can simply introduce a matter component Qm of the 
nonpotential force [42]. The matter supermomentum components are Ti^ = 
-2kgiT\ . 

The components of the gravitational supermomentum in almost synchronous 
gauge are defined by [2,42,43] 

= 2gi „ = P{Sa) , S^ = K- 2atS\e\ 
Tr 5„ = , [<5„ , h\ = {eabdn''^ + 3af6{^^)S, , 

where the adjoint matrices kc, = C^ac^'^b are introduced in appendix A and given 
explicitly by (A. 3). The matrices Sa generate a subgroup of SL{3,R) which 
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coincides with the hnear adjoint group in the class A case and a projective 
automorphism subgroup in the class B case [2]. These matrices are linearly 
dependent for Bianchi types I {Sa = 0), type II {Si = when n^^^ ^ 0, etc.) and 
type VI_i/9 (jn^"'^-' 1 2^]^ + |n(^)|252 — 0) where the supermomentum constraints 
are degenerate, imposing additional constraints on the matter supermomentum. 

The Einstein evolution equations in almost synchronous gauge are then the 
equations of motion of the total Lagrangian/Hamiltonian system with the non- 
potential force Q and total Lagrangian and Hamiltonian 

L = Lg- NHm , H = Hg + NHm = NH , (2.52) 

namely 

-6L/5gab = (dL/dgab) ' - dL/dgab = NQ^' 

gab - {gab, H} , n'^" = {tt^\ H} + NQ'^" . 
These are subject to the constraint equations 

n = nG + nM = 0, n^ + H^ = o. (2.54) 

For the present paper the source of the gravitational field will be assumed to be 
a perfect fluid. The appropriate choice of variables and expressions for the fluid 
super-Hamiltonian, supermomentum and equations of motion are discussed in 
appendix C. 



3 Diagonal Gauge as an Almost Synchronous 
Gauge Change of Variables 

The Einstein equations for a spatially homogeneous spacctime in almost syn- 
chronous gauge have been put in the form of a C^i-parametrized constrained 
classical mechanical system driven by the matter variables which themselves 
satisfy certain equations of motion. The combined equations of motion are 
invariant under the action of constant elements of the automorphism matrix 
group Aute{g) representing the freedom remaining in the choice of the spatially 
homogeneous spatial frame in almost synchronous gauge for each fixed point 
C"*6c € Cd- However, only the special automorphism matrix group SAute{g) 
acts naturally on the Lagrangian/Hamiltonian system as a symmetry group, 
corresponding to the additional restriction that the 3-form o;^ A A oj^ remain 
invariant under change of frame, so in this context it is SAute{g) rather than 
Aute{g) which plays an important role [3]. The supermomentum constraint 
functions are associated with this symmetry, although the direct connection is 
not so obvious when the nonpotential force is nonzero. 

Whenever a dynamical system has a symmetry group, one may simplify the 
system by choosing new variables adapted to the action of the symmetry group 
on this system, a very instructive example being the central force problem. Ex- 
actly how to adapt the variables depends on the particular way in which the 
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symmetry group acts on the system. For all points of Cd except the set of mea- 
sure zero occupied by the type I, II and V orbits, the group G = SAute{g) is 
3-dimensional, closely connected to the three linearly independent supermomen- 
tum constraint functions, and has an offdiagonal matrix Lie algebra g which is 
such that any point g€ Ai may be diagonalized by one or more elements of this 
group acting as in (2.16). For the remaining three Bianchi types, SAute{g) has 
a larger dimension and for types I and II there are fewer linearly independent su- 
permomentum constraint functions, but there do exist families of 3-dimensional 
subgroups G with offdiagonal Lie algebras which may be used to diagonalize an 
arbitrary point of A4. Recalling the advantages of diagonal metric component 
matrices suggests that Md should be used to parametrize the orbits of the 
action of these 3-dimensional groups G on A^, leading to the following class of 
parametrizations S) e diag{3, R) x G ^ A4 

g = /-i(e^/') = /-i(gO = S^e^^S, (3.1) 

which decompose the metric matrix variables into "diagonal" and "offdiagonal" 
variables as described in the introduction. The space diag{3, R) ^ Diag{3, i?,)^ 
^ Md representing the diagonal variables has already been parametrized in 
various ways, leaving to be discussed the parametrization of G as well as its 
choice for the type I, II and V cases. For these latter types and the remain- 
ing nonsemisimple types, some of the diagonal variables are associated with 
additional diagonal automorphisms. 

The property that the orbit oi A4d under the action of G be A1 (in order 
that any point g £ M may be represented in the form (3.1)) requires that the 
offdiagonal matrix Lie algebra g have an ordered basis {/Ca} with the property 
that for each cyclic permutation (a, 6, c) of (1,2,3), then Kq G span{e''c, 6^5}. 
Consider the following Lie algebra basis valued function {/Ca} on Cd, defined by 



e«" =2-^(k°,kr)i =2-i[(n(''))2 + (n(^))2]i 

and satisfying 

u = diag(n(i),n(2),n(3)) , n^") = n^'^^e^" 



(3.2) 



(3.3) 



where it is clear from the context when (a, 6, c) is to be interpreted as a cyclic 
permutation of (1,2,3). This is well defined everywhere on Cjj except for the 
type I, II, IV and V orbits (precisely those points of Cd where rankn < 2 
and the scale matrix e" = diag(e" , e" , e" ) is singular) where it has direction 
dependent limits. Consider those points of Cd for which n^^^ = 0, for example. 
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One then has 

e° = 2-5diag(|n(2)|, \n^^)\,\{n^^^f + {n^^^f\i) 
{«„} = {-\/2sgnn(2) e^,\/2sgnn(^) e^, \/2(- cos(?ie^ + sin^!)eS)} 
n = \/2sgnn^^^n^^^diag(sin(^,cos^0) 

(cos0,sin(/>) = |(n(i))2 + {n^^^f\-^n^^\n^^^) . 

For type IV only the sign of either Ki or K2 does not have a well defined limit 

for a given orbit, so {/Ca} has two values (but a unique G) for each orbit. In 
addition to this sign indeterminacy (which does not make G multivalued), an 
S^-paramctrized family of limits exists for each type II and V orbit component 
(associated with the indeterminancy of the matrix for the type V orbit and of 
the matrix Ka for the type II orbit component on which n(°) is the only nonva- 
nishing structure constant tensor component), while an 5^-parametrizcd family 
of limits exists for the single type I point containing all well-defined values of 
this function. Thus {Ka} is a multivalued Lie algebra basis valued function on 
Cd whose values at each point determine the offdiagonal diagonalizing automor- 
phism matrix subgroup Lie algebras g. G is then a multivalued matrix Lie group 
valued function on Cd- (The space of values of this function is diffeomorphic 
to S*^ modulo reflection about the origin, namely P^.) For the types I, II, IV 
and V wlu^rc the multivaluedness occurs, one may pick any value to describe 
the dynamics. 

At the canonical Bianchi type IX point of C^i, G has the value 5*0(3, R) and 
the parametrization (3.1) reduces to the one introduced by Ryan for all Bianchi 
types [5]. The latter parametrization arose from considerations completely un- 
related to symmetries but by a fortunate coincidence agrees with the correct 
choice for the most interesting case: Bianchi type IX with canonical structure 
constant components. Given any frame e with metric component matrix g, 
there is a natural orthonormal frame related to e by the symmetric square root 
of g (unique up to ordering when the eigenvalues of g are nondegenerate) 

g = B\ B\ = B\, e': = B-'\ea. (3.5) 

Since any symmetric matrix can be diagonalized by an orthogonal transforma- 
tion, (B = O^BzjO ,0 e 50(3,7?) ,Bd e Diag{3,R)+), using the property 
one obtains the result g = O B^O. By setting = e'^, one 
arrives at the Ryan parametrization and its associated orthonormal frame 

g = O^e'^O , e': = {O^e^Or^e, = 0-^'=Je^O)-",e, . (3.6) 

The incompatibility of this parametrization with the symmetry at all points of C 
except for the type I point (where it is unnecessary for perfect fluid spacetimes) 
and the type IX points where n is proportional to its canonical value made its 
application to other points of C ineffective. (Compare the expressions of the 
first paper of ref.(23) for the nonpotential force with (3.29).) 
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If needed one may use the following parametrization of G 

S = e^'^V^'^^e^'"^ . (3.7) 

Suitably restricting the domain of the parameters {9"-} when G has compact 
directions leads to local canonical coordinates of the second kind on G (which 
are always valid in an open neigborhood of the identity if not globally) and 
hence through (3.1) to local coordinates on each 3-dimensional orbit of G on 
Ad. These local coordinates on M bc;come singular on G-orbits of dimension less 
than three, which represent fixed points of the action of the compact subgroups 
of G. The basis has the property that when the single parameter O"" for 
some fixed value of a is nonzero in (3.7), then (3.1) parametrizes the symmetric 
case submanifold M.s{a) = {g € | Qab = Qac = 0}, where as usual (a, 6, c) is 
a cyclic permutation of (1,2,3). These three submanifolds are associated with 
discrete spacetime symmetries [43]. 

Through the parametrization (3.1), the action of G on Al 

g = /s'(g') - /a (g) = /SA-(S') , A e G , (3.8) 

becomes inverse right translation on G itself. Since the DeWitt metric is 

SL{3, i?) invariant and G C SL{3, R) by virtue of having a tracefrce Lie algebra, 
if a right invariant frame on G is employed, the components of the restriction of 
the DeWitt metric to the orbit (identified with G) will have components which 
can at most depend on the diagonal variables. 
The relations 

. . (3.9) 

define the left and right invariant frames {Ea} and on G with respective 

dual frames {W°-} and {W"} and structure functions G"hc = W"-{[Ea, Eb]) = 
—W°'{[Eb, Ec]) which are naturally associated with the basis {/Co} of its matrix 
Lie algebra. They satisfy [Ea,Ei)\ = 0. If S(t) is a parametrized curve in G, 
then 

S{t)-^S{t) = KaW"(i) , S(i)S(i)-l =K„l^"(i) (3.10) 

defines the component functions W°'{t) and W°'(t) of the curve's tangent vector 
with respect to these frames. Let the corresponding functions on the cotangent 
bundle T*G be denoted by Pa and Pa respectively. If {O""} are local coordinates 
on G, then one has the following coordinate expressions for the right invariant 
frame quantities 

Ea = E\{e)d/de'' , VF« = w\{6)de' , iw\) = {e\)-' , 

W'^ = W\d'> , Pa = E\pb , 

where {0°', ^)"} and {6,pa} arc the natural lifted coordinates on the tangent and 
cotangent bundles. Dropping the tildes leads to the corresponding left invariant 
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frame quantities which are related by the Hnear adjoint transformation on G 



W" = R\W^ , Pa = PbR-\ , etc. . 

One may also use the local canonical coordinates {9°-,Pa) to evaluate the follow- 
ing Poisson brackets 

{Pa, Pb} = C\bPc {S, Pa} = «aS 

{Pa,Pb} = -C\bPc {S,Pa} = SKa {Pa,Pb} = Q. 

Since the basis {/Ca} is closely related to the standard basis of one of the standard 
diagonal form adjoint matrix Lie algebras with rank n > 1, explicit expressions 
may easily be obtained from the formulas of appendix A for S, R and the 
component matrices of the invariant fields in terms of the paramctrization (3.7). 
The case n^^^ = is given as an example in that appendix. 

The paramctrization (3.1) has the following meaning. Given the spacctime 
metric (2.24) is some almost synchronous gauge, i.e. given the parametrized 
curve g(Z) in Ai and the lapse function N(t), then A(t) = S(i) is the matrix in 
(2.27) which affects the change to diagonal gauge, where the new metric matrix 
g(t) — g'(i) is diagonal and the new spatial frame 

e; = S-\{t)eb (3.14) 

is orthogonal, with the associated shift vector field satisfying 

ade'(iV(t)) = /c„M^« , (3.15) 

according to (2.33). The matrix e Diag{3, R)~^ then normalizes this orthogo- 
nal spatial frame, leading to the natural symmetry adapted orthonormal spatial 

e^' = {e^)-\e', = {e^S)-\eb (3.16) 

which may be used to introduce spinor fields on the spacetime in "time gauge" 

[42,56,57] or to introduce a natural Newman-Penrose null tetrad or to put the 
Einstein equations in a simple form without a Lagrangian/Hamiltonian formu- 
lation [58]. The offdiagonal velocities W^" are linearly related to the angular 
velocity of this natural spatial triad relative to one which is parallelly propa- 
gated along the normal congruence [42]. Note that in the canonical type IX 
case, this orthonormal spatial frame differs from the frame (3.6) which is asso- 
ciated with the Ryan paramctrization by an additional rotation, and like that 
frame, is unique modulo ordering of the frame vectors and barring degeneracies 
for each value of G. 

Let a prime indicate components with respect to the diagonal gauge spatial 
frame (3.14). As noted in the previous section, spatial curvatures have simpler 
expressions in diagonal gauge; in particular the scalar curvature potential func- 
tion Ug is independent of the offdiagonal variables and is simply given by the 
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formula (2.27). The extrinsic curvature (2.36) and the kinetic energy have the 
following expressions when evaluated with respect to the primed frame 

K' = SKS-i = -N-\0 + K* 'Wn , K*' = \{Ka + e-2^«„ ^e^^) 
NT = 6^/3° (K',K')^^ = N-\^^\Qr,AB0^^'' +gabW''W'') 

Gab = (k* 'jHf = {nf ',Kf ) , 

(3.17) 

indicating that the DeWitt metric itself is 

Ig = e^'f'" {6vABdl3^ O +gabW'' (S> W'') . (3.18) 

The components of the rescaled DeWitt metric Q = jg~iQ along the orbit 
directions are functions on the f3^P~ plane and are diagonal due to the choice 
of basis (3.2) of g 

Gaa = ie-2""(nWe^" - nWe-''")^ ^ (^ ^ ■ (3.19) 

When (e-"°n(''),e-"''n('=)) equals respectively (\/2,0), (1,-1) and (1,1), as 
occurs at canonical points of C^j, this expression has the values e"^^ " , 2cosh^/3'"^ 
and 2sinh^/3^'=. When is a compact generator, which means that e~" n^^^ 
and e~"°n(°) are nonzero and of the same sign, then Qaa vanishes for 

Pa = Pao = -(4\/3)-i In |e-«''n('')/e-«"n(^)| (3.20) 

which for canonical points of Cd represents the 2-dimensional orbit of the 

Taub submanifold A^T(a)- Such points of A4 represent singularities of the 
parametrization (3.1). This same condition on /?~ picks out the submanifold 
of Md for which dQaa = = rfa"^"" when e-^^n^*) and e-«"n(^) are both 
nonzero. 

The velocity-momentum relations following from the above kinetic energy 

are 

PA = d{NT)/d$^ = Ue^'^'"N-^f]AB$'' , - (12e3'3'')-i7Vr,^^pB 

Pa = a(ivr)/avK" = 2e^'^''N-^gabW'' , = ^e"^'^" Ng-^'^^Pt, . 

(3.21) 

Re-expressing the gravitational momentum (2.38) in terms of these momenta 
using (3.15) leads to 

tt' = SttS-i = ^(ry^^p^es + 3poeo) + ^^-'^'Pawf ' • (3.22) 

Defining 

P'(A) = -2Trir'A , A G gl{2,, R) , (3.23) 

one finds 

Pa = -P'iBA) , Pa = -P'(«a) , (3-24) 
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the minus sign arising from the inverse aetion used in the parametrization (3.1) 
to conform with other conventions. Note that the canonical generators of the 
action (3.8) of G on are instead (using (2.43), (3.12), (3.22) and (3.23)) 

P(/Ca) = P'(SKaS-l) = P'{Kb)R\ = -Pa , (3.25) 

namely the left invariant frame momenta. Finally, re-expressing the kinetic 
energy as a function on momentum phase space leads to 

NT = iiV(12e3'5'')-i(77^VpB + QQ-^'''' PaPb) (3.26) 

The action (3.8) of G is a symmetry of the free Hamiltonian dynamics generated 
by the kinetic energy function alone and hence the canonical generators Pa are 
conserved; this is clear from (3.26) since Pa commute with NT provided that 
A'' is G-invariant. The momenta Pa are related to the conserved momenta by 
the time dependent adjoint transformation. 

To summarize, {d/d0^,Ea} represents through the parametrization (3.1) 
an orthogonal frame on A4 adapted to the thrcc-plus-thrcc decomposition of 
M into orbits of the action of G (offdiagonal or "angular" variables) and their 
orthogonal submanifolds (diagonal or "radial" variables). In Misner's super- 
time time gauge N = 12e^^ , the diagonal part of the kinetic energy is simply 
the standard kinetic energy of the flat Lorentz geometry of the metric space 
{diag{3, R), ((, ))ow ) expressed in inertial coordinates {P^}, while the offdiago- 
nal part of the kinetic energy is just the one associated with the right invariant 
metric / = ^Gab^'^ ® W'' on G, the two pieces coupled together by the (3^ 
dependence of the components lab = \Qab- 

The offdiagonal dynamics is thus governed by the natural generalization 
of the rigid body dynamics at the canonical type IX point of where G = 
SO{3,B) to the (multivalued) matrix group G at each point of Cd- This is 
discussed at length for a general Lie group by Abraham and Marsden in a more 
abstract notation [46]. The spatial frames e and e' are respectively identified 
with the space-fixed and body-fixed axes, related by the passive transformation 
S which corresponds to Goldstein's matrix A of his equation (4.46) [45], the 
value assumed by S for the canonical type IX point of when expressed in his 
Euler angle parametrization of S0{3,R). Here (-VF°,-VF") and {-Pa, -Pa) 
play the roles of the space-fixed and body-fixed components of the angular 
velocity and spin angular momentum of the rigid body, and in a general time 
gauge, /(a) = 2N~^e^^ Qaa play the role of the principal moments of inertia, 
namely the eigenvalues of the moment of inertia tensor /. Notice that it is 
the space-fixed components of the angular momentum which are conserved by 
the free dynamics, exactly as in the force free rigid body case; however, it is 
the principal (body-fixed) axes of the moment of inertia tensor will allow one 
to solve the equations of motion for that case. It is this fact which motivates 
diagonal gauge in the present problem. 

The identification of S with the passive coordinate transformation from 
space-fixed to body-fixed coordinates leads to the minus sign which arises here 
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and in (3.24). The paramctrization (3.1) with / rather than would lead to 
the identification of S with the aetive transformation from the space-fixed basis 
vectors to the body- fixed basis vectors (Goldstein's matrix A~^), eliminating 
the mimis sign and interchanging left and right in the above discussion. For 
the type IX case, this rigid body analogy is simply the restriction to symmetry 
compatible diffeomorphisms and subsequent translation into frame language of 
the Fischcr-Marsdcn discussion for general spatially compact spacetimcs. The 
finite dimensional situation here allows the analogy to be carried much further. 

On the other hand, the left invariant frame momenta Pa generate the canon- 
ical action of G on the momentum phase space and thus act like orbital angular 
momenta in the central force problem, leading to the analogy in which the radial 
and angular variables correspond to the diagonal and offdiagonal variables of 
the present system. The fact that this action of G arises from the action of spa- 
tial diffeomorphisms on the spacetime metric (2.31) makes the decomposition 
into diagonal and offdiagonal variables relevant to the Einstein equations. 

A key feature of this decomposition is the resulting trivialization of the 
supermomentum constraints (which are in general intimately connected with 
the spatial diffeomorphism group) and of the nonpotential force. These involve 
the matrices {Sa} introduced in (2.51) and given explicitly by 

51 = — n(^)e^2 + n''^^e^3 — 3ae'''i= e"^Ki — 3a2^^ sgn{n'-^^)K2 

52 = -n(3)ei3 + n(i)e3i - ■Sae^2= e"'/C2 + 3a2-^ sgn(n(2))/ci (3.27) 
^3 = — n^^^e^i + n^^^e^2 + ae+ = e"' /C3 + ae+ , 

which leads to the definition of a matrix p by 

Sa = Kbp'^a + aae+ . (3.28) 
The nonpotential force is then evaluated as follows 

Q = 2e'^\^S,,g-'dg) = 2ae''°+^''USs,Ss-'dgS-') 
= 4ae^°+^^^ {6adp+ + e"%3W'') = Q+dl3+ + QaW^ . 



(3.29) 



Note that the /3+ component of the nonpotential force exactly cancels the /?+ 
component of the Einstein force arising from the exterior derivative of the 
term in the expression (2.27) for Ug 

Q+ - d/dl3+i6a''e^"+^^*) = , (3.30) 

i.e., only the first term e^°V* is relevant to the equations of motion for the 
variables. (Subtraction of this second term from Ug in the Lagrangian 
or Hamiltonian would require the introduction of a component of Q; if /S" 
is determined by integrating the super-Hamiltonian constraint rather than by 
equations of motion, this is irrelevant, leaving only the component Q3 as relevant 
to the remaining equations of motion and which vanishes for Bianchi type V.) 
The primed gravitational supermomentum components are given by 

' = P'iSa) = -Pbp\ - aaP+ . (3.31) 
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As already noted, these are linearly dependent for Bianchi types I, II and Vl^i/g. 
This is reflected by the vanishing of det p = [e" e" +9a^sgn(n(^'n(^')]e" which 
has the class A value e«'+«'+"' and the class B value |n(i)n(2)|(l + 9/i)e"'\ This 
also vanishes for Bianchi type V where e" = 0. 

For the remaining types one may solve the supermomentum constraints for 
the offdiagonal momenta or velocities 

Pa = -p-\{nf ' + a,p+) = p-'\{n^' ' - a,p+) 

and thus determine the diagonal gauge shift vector field (3.15) modulo time 
dependent elements of g (spatial Killing vector fields). This is essentially just 
an application of the thin sandwhich formulation of the initial value problem 
[59,60]; initial values of (g',g',7i^ ') are specified arbitrarily and the supermo- 
mentum constraints are used to determine the diagonal gauge shift variables. 
For Bianchi type V the diagonal/off diagonal decomposition is not compatible 
with the thin sandwhich interpretation since the diagonal momentum p+ rather 
than the offdiagonal momentum P3 is constrained, while for types I, II and 
VI_i/9 certain shift components (i.e. offdiagonal momenta) remain freely speci- 
fiable initial data due to the degeneracy of the constraints. The important 
distinction between the present case and the thin sandwhich interpretation is 
that here the shift variables have been made velocities associated with dynam- 
ical variables and the supermomentum constraints allow one to determine the 
shift variables directly at each moment of time rather than by integrating their 
(first order) equations of motion. 

The idea of incorporating the shift variables into the dynamics as the ve- 
locities associated with gauge variables was developed by Fischer and Marsden 
[46]. Their construction when adapted to the spatially homogeneous case ex- 
tends the configuration space to x i? x Aute{g) with natural variables 
(g, A) and velocities {g,i = N,AA~^ = ade(iV)) whose interpretation is the 
same as in (2.31)-(2.35). Imposition of the G-diagonalization gauge condition 
(g, ade(.^)) e A^D X ^ on this extended system then leads to the results of the 
present approach. 

Thus the diagonal/offdiagonal decomposition of the almost synchronous 
gauge gravitational variables leads automatically to the diagonal gauge equa- 
tions of motion, while the diagonal gauge supermomentum constraints then 
determine the diagonal gauge shift variables which appear in these equations 
of motion, except in the degenerate cases (Bianchi types I, II, V and VLi/g) 
where first order equations of motion for the undetermined shift variables must 
be used to evolve them. These first order equations may be obtained directly in 
Hamiltonian form by using the Poisson bracket relations (3.13). (For the vac- 
uum dynamics these are just generalized Euler equations with time dependent 
principal moments of inertia, driven by the additional nonpotential force com- 
ponent Q3 in the class B case.) To complete this scheme one needs to introduce 
the diagonal gauge matter variables by a time dependent change of variables 
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involving the transformation matrix S, leading to the appearance of shift terms 
in the new matter equations of motion through (3.15). These terms are analo- 
gous to the centrifugal force which appears in the rotating body-fixed frame of 
the rigid body problem. 

The combined Einstein and matter equations then involve only diagonal 
gauge variables and the diagonal gauge shift variables which are in general de- 
termined by the supcrmomcntmn constraints (or by first order equations of 
motion if not). The spacetime metric may then be represented directly in di- 
agonal gauge form, without integrating the equations which determine S, the 
transformation back to almost synchronous gauge. The offdiagonal part of the 
kinetic energy then acts like a time dependent effective potential for the diago- 
nal gauge gravitational variables, which from the central force analogy has been 
called the centrifugal potential by Ryan [20,21]. The matter super-Hamiltonian 
acts as another time dependent potential for these variables. 

For a perfect fluid source as discussed in appendix C, the appropriate vari- 
ables arc {n,l,Va); the diagonal gauge variables are defined by (C.3) with 
A = S, namely 

{n',l',v',) = {n,l,VbS-'\) . (3.33) 

The primed components of the matter supermomentum and the matter super- 
Hamiltonian then have the expressions 

2kl{^i^ + g'°^>;)5 - 2kpe'^'^'' (3.34) 

In the Bianchi type I dust (zero pressure) case, Ryan has called the first term in 
the matter super-Hamiltonian the rotation potential since it directly influences 
the gravitational variables only when the fluid is tilted {va 7^ 0) [61] and tilt 
is equivalent to rotation of the fluid in this case. However, for the remaining 
types tilt and nonzero rotation {C^^Vb 0) are no longer synonymous [43], so 
tilt potential is a better name. Let designate this potential. The second 
term or pressure potential only affects the equations of motion for 

To evaluate the equations of motion for the offdiagonal momenta, one needs 
the following Poisson brackets which are a consequence of (3.13) and (3.33) 

{Pa, V'^} = v',Ka\ {P„, V'^} = v'.Kf.R'a • (3.35) 

One then finds the result 

(P„)- = \Ne-^P°P,C%a'Q-^^'^Pd + N{Qa + Fa) 

= P,C%aW'' + N{Qa + Fa) 
{PaY = {PbR''ay = N{Qa + FbR\) 

Fa = {Pa,W^*} = 2kliv^)-W,Ka^v" 

For the free dynamics and class A vacuum dynamics {Qa = = Fa), these 
equations generalize the Euler equations for the body-fixed components of the 
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angular momentum of a rigid body, and like those equations conserve the space- 
fixed components of the angular momentum. For the class B vacuum dynamics 
the component Q3 appears as a driving term and the conserved quantities are 
instead the unprimed supermomentum components 

= p{Sa) = p'{Sb)s\ = p'isSaS-') 

= P'{SKbS-')p\ + aaP'{Se+S-') 

= P'{Kc)R-\p\ + aaP'{e+ - 3(0i/ci + ^^kj)) ^^'^ 
= -Pbp\ - aa{p+ - S{0'Pi + 0^P2)) . 

The momenta Pi and P2 arc conserved by the vacuum dynamics, being constant 
linear combinations oiHi and H2 , but P3 is not unless Q3 = 0. Because of the 
supermomentum constraints (3.32), the equations of motion for the offdiagonal 
momenta Pa are very closely related to the primed fluid equations of motion 
which may be obtained by transforming (C.5) 

i = Nliv^)-^2a'^v', {v'J =Nv',[{v^r'C'',av''-K,\W^] (3.38) 

The lapse A'' is still arbitrary and the super-Hamiltonian constraint remains 
to be dealt with. Several options are available. Introduce the function 7^ 
depending on an arbitrary real parameter h by 

7^ = -(sgnpo)[24e3^"(H + h) + Po"]^ , (3.39) 

enabling the super-Hamiltonian constraint to be expressed in the compact form 

n^^e-'''^"{-Po^+I,^) = , (3.40) 

which suggests that it be used to solve for po, the result being po = —Iq- {po 
is negative for expansion or increasing /3° with increasing proper time and is 
positive for contraction.) 

Suppose one initially assumes unit lapse function, so that the time of the 
classical mechanical system coincides with the proper time t on the spacetime. 
The super-Hamiltonian constraint may then be used to eliminate another degree 
of freedom from the system by a standard technique of classical mechanics [1] 
which replaces the time t as the integration variable by some function t of 
(accomplished by changing the independent variable in the action). This is 
called an "intrinsic time reduction" of the system [62]. The usual and obvious 
choice is t = f3^ = — f2 or negative f2-time. (il-time is useful for approaching the 
initial singularity.) This is clearly equivalent to specifying a new lapse function 
. The reduction technique leads to the reduced Hamiltonian Iq and equations 
of motion 

/3± = 7o} , p± = {p±,Io} + 5+±NQ+ , (3.41) 

where the new lapse function N is determined by the first order equation for 
the original time t 

N = i = dijdh 1^^,= 12e3'5°7o-^ (3.42) 
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and the dot refers to the new time derivative. 

This same result can also be obtained by simply setting TV = —12e^^ Pq~^ 
in the original Hamiltonian H = NTi 

/3° = {/30,i/} = i(l+Po-2V), ^ • ) 

the result = 1 following from the imposition of the constraint po = —Iq. 
For any variable y other than po, one has dH/dy = {—pQ~^Io)dIo/dy, showing 
that Iq acts as the Hamiltonian for all variables but /S" when the constraint 
Po = —Iq is imposed. Misner's supertime time gauge N = 12e^^ differs differs 
only by eliminating the factor —po in the /3°-time lapse function, leading to a 
Hamiltonian of familiar form rather than a square root Hamiltonian as above 

H = -W + W- (3-44) 

In this time gauge, ^Iq^ acts as the Hamiltonian for the variables other than 
[P'^iPo) and the super-Hamiltonian constraint /3° = —po = Iq may be used to 
evolve Note that 



(3.45) 



The term Hgff acts as an effective Hamiltonian for the variables other than 
(/3°)Po); the corresponding Hamiltonian equations of motion for the degrees 
of freedom contain no nonpotential force. 

Whatever the time gauge, one may introduce an effective Hamiltonian for the 
(3^^ degrees of freedom by subtracting out the p^^ term, the pressure potential 
and the term in the scalar curvature potential 

H^ff =H- N[-±e-^0"po^ + 6a^e^°+*^^ - 2kpe^0"] 

= NiUe^^^Y^i^ip^ + + 3a-'"'Pan + 12e*^V* + 24e'f'"n'^'] . 

(3.46) 

As discussed above, the use of H^ff rather than H eliminates the need to 
consider the nonpotential force component Q^, so H^jf alone generates the 
correct equations of motion for the degrees of freedom. The component Qs 
is still relevant to the equations of motion for the offdiagonal variables, when 
needed, for which either H or -ffe// may bo used. A nonvanishing value of Qt, 
is connected with the single nonvanishing offdiagonal component of the spatial 
Ricci tensor in diagonal gauge (see (2.25)). 

One advantage of the supertime time gauge is that the kinetic energy for the 
diagonal variables corresponds to a relativistic particle in the flat 3-dimensional 
spacetime {diag{3, R), r^Asd^'^ d(3^). Neglecting the remaining terms in the 
Hamiltonian (the free system with zero orbital angular momentum) , one obtains 
afflnely parametrized geodesies for the system trajectories and the Hamiltonian 
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constraint just requires that the particle velocity be null. Except for small 
anisotropy in the type IX case, the remaining terms are all positive (although the 
pressure potential is negative, the total fluid Hamiltonian is positive) and their 
effect is therefore to make the particle velocity timelike. It is future pointing in 
the case of expansion increasing) and past pointing in the case of contraction 
decreasing), reversing direction only in the type IX case where recollapse 
occurs. 

The /3°-time gauge, on the other hand, uses the natural time variable /3° 
of the 3-dimensional flat spacetime to parametrize the system trajectories; for 
the geodesies it is also an affine parameter. This has the advantage of making 
the reduced Hamiltonian for the degrees of freedom explicitly rather than 
implicitly time dependent. 

The curvature potential V* has been described and diagrammed in the pre- 
vious section and enters the Hamiltonian in the supertime time gauge as the 
time dependent potential Ug = 12e^^ V*, which apart from constants has been 
called the gravitational potential by Ryan [20,21]. The centrifugal potential 

.3 3 

Ue = ^g-^^^PaPb = 3 ^ = (3.47) 

n=l n=\ 

is the sum of three terms whose dependence is given by (3.19) and is repeated 
in a more suggestive form here. 



(i) e-2ar„(b)„(c) > Q . g-iaa ^ 1 1 e-2"° „W „M | "iginh-^ (2^/3(/3- - ^-q)) 
(a) e-2"''n('')n(^) < : G'^'"' = i|e-2"''n('')n('=)|-icosh-2(2V3(/3- - /J-g)) 

2\/3/3-o = -iln|e-«"nW/e-«°n('=)| . 

(3.48) 

These are illustrated for the canonical points of Co in Figure 7. The po- 
tentials for noncanonical points arise from these by a combined rescaling and 
translation, infinite translations corresponding to Lie algebra contractions. All 
of these potentials are exponentially cut off in one or both directions. Case 
(?) corresponds to a compact generator Ka in which the potential behaves like 
(Pa )^^ s-t the origin, exactly as in the central force problem where the cen- 
trifugal potential has the dependence on the radial coordinate r. Thus the 
singularities of the parametrization (3.1) are shielded by angular momentum 
barriers exactly like the spherical coordinate singularity at the origin in the 
central force problem. Case (ii) also acts like an angular momentum barrier 
but can be overcome if the energy of the system is sufficiently large. In case 
(iii) the barriers have been translated out to infinity. In the supertime time 
gauge, all of the time dependence of the centrifugal potentials arises from the 
angular momentum factors. Increasing increases C/i"^ and causes a given 
value of the potential to move in the direction of decreasing values. The barrier 
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Figure 7: Plots of Q for canonical points of Cd, at which the pair 
e~"°(n^''\ n^'^^) assumes the values shown. 

(i) may be crossed provided that the relevant angular momentum component 
vanishes at the moment of crossing. 

For the canonical points of Cd, the gravitational potential (in supertime 
time gauge) is closely approximated by the simple time dependent Bianchi type 
II exponential potential J7g"^ = 12n^°-'^e^'^ -8/3^ between the positive P'^- and 
/3+-axes and outside the channels. Because of the simple exponential depen- 
dence of this potential, a given value of this potential has a constant value of 

— 2/3+, i.e. its contours move with velocity d(3'^ /d(3^ = ^ in the flat 3- 
dimensional spacetime {diag{3, R), rjABdfi^ ® d[3^). Similarly the tilt potential 
(in supertime time gauge) in the limit of large anisotropy or extreme tilting 
along the diagonal axis is closely approximated by the exponential potential 
^tiit — '^■'^kle^'^^°^^i^\vj_\ between the negative (3^- and /3+-axes, and hence its 
contours move with velocity d[3'^ /d(3^ = — 1 plus whatever additional velocity 
is due to the time variation of v'^ (and / in the class B case). These exponential 
potentials are indicated by single equipotential lines in Figures 8 and 9, the di- 
rection of decreasing values (the direction of the associated force) indicated by 
arrowheads on these lines. Figure 8 provides a key which indicates the labeling 
of the various potentials of Figure 9. Although the tilt equipotential locations 
are shown as symmetrical, decreasing vj[ translates U^^^^ in the positive di- 
rection, so that when an individual (primed) component of the fluid current 
vanishes, the corresponding tilt potential moves out to infinity. A similar state- 
ment holds for the exponential gravitational potentials and their dependence on 
|n*^°^|. One must actively translate these canonical potentials by the translation 
(and rescaling) (2.30) to obtain the noncanonical potentials, with Lie algebra 
contraction causing the potentials to move out to infinity. Equipotential lines of 
these potentials move with /3°-velocity ^ in the direction of decreasing values, 
while the tilt potentials in the exponential regime move with /^''-velocity 1 in 
the direction of decreasing values neglecting the additional time dependence due 
to the fluid variables. The centrifugal potentials in supertime time gauge have 
no explicit dependence on (3^ and hence their motion is entirely due to the G- 
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Figure 8: The key to Figure 9. Together with Figure 2, this shows the 
coordinate associated with each potential represented in Figure 9 by wall con- 
tours. In general the wall contours associated with tilt and type (Hi) centrifugal 
potentials will be translated from their symmetrical positions. The dotted lines 
for the type (i) and (ii) centrifugal potential contours in Figure 9 have double 
arrows since they actually represent a pair of contours symmetrically located 
with respect to the point at which the potential becomes infinite or assumes its 
maximum value. 

angular momentum factor. When this latter factor is constant as occurs in the 

class A symmetric cases to be described shortly, these potentials are stationary. 
The type VI and VII potentials in each sector between the negative (3^- and 
/3+-axes, namely ?7g(VII,n(«) = 0) = 2|n('')n(^)|e*(^°+^a )sinh22\/3(/3- - /J-q) 
and C/g(VI,n(") = 0) = 2\n^^^n'''''^\e'^''f^''+f^ihosh^2y^{j3- - P'o) simply trans- 
late along the /3+ axis with velocity = — 1 in supertime time gauge, 
so the channels move with unit /3°-velocity in this time gauge (outside their 
common intersection at semisimplc points of Co)- 

In the /3°-time gauge, all of the potentials become multiplied by an addi- 
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Figure 9: The walls associated with the gravitational, tilt and centrifugal poten- 
tials for the canonical points of Cd ■ Arrows indicate the direction of decreasing 
values of each potential. Although shown as symmetrical about the origin for 
simplicity, the tilt and type (iii) centrifugal walls will in general assume unsym- 
metrical positions. Figure 8 explains the conventions for this figure. 
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tional factor of (—po)~^, leading to additional implicit time dependence and an 
additional contribution to the velocities of their contours. 

Before discussing the qualititive effect of the various potentials on the dy- 
namics, it is important to have in mind a classification of the specializations 
which can occur. Because of the existence of various additional continuous (lo- 
cal rotational symmetry and isotropy) and discrete spacetime symmetries, the 
dynamics admits various special cases where one may consider the phase sub- 
space associated with a submanifold of the configuration space such that initial 
data in this phase subspace remains in it when evolved by the equations of mo- 
tion. These special cases arc important since they arc governed by much simpler 
systems of differential equations which often admit exact solutions. These ex- 
act solutions are the basis of analytic approximation schemes or "diagrammatic 
solutions" . 

Since Aute{g) is a symmetry group of the equations of motion, any such 
phase subspaces which are related by the action of an automorphism are equiv- 
alent and it suffices to consider only one representative from each equivalence 
class. These have been described elsewhere for Dirac spinor [42] and perfect 
fluid [43] sources for the canonical points of Cd and are summarized for the 
perfect fluid case in Table 3. 

For types I and IX all symmetric and Taublike cases (and for type I all 
Taublike symmetric cases) are equivalent, while only two inequivalent symmet- 
ric cases exist for the remaining class A types. The existence of nontrivial 
continuous isotropy subgroups of the action of Aute{g) on ^A for types I, II and 
V leads to the following additional equivalences for perfect fluid sources. For 
type I the general case is equivalent to the diagonal case. For types II and V 
the symmetric case M s{3) is equivalent to the diagonal case, while the Taublike 
type V case is equivalent to the isotropic case when ^3 = 0. All Taublike cases 
except for Bianchi types VIo and VI^ are associated with local rotational sym- 
metry and all isotropic cases with spatial isotropy. There also exist equivalences 
between diflferent Bianchi types [42]. The type VIIq and type I Taublike cases 
are equivalent, as are the type Yllh and type V Taublike cases. The symmetric 
cases correspond to the situation in which the G-angular momentum is aligned 
with one of the body-fixed axes, i.e. Pa has only one nonvanishing component 
(a class A constant of the motion); this preferred body- fixed axis coincides with 
one of the space-fixed axes and is a Kasner axis. The remaining two Kasner 
axes move in the plane orthogonal to this axis. The diagonal case corresponds to 
vanishing angular momentum, in which case the frame e = e' is a Kasner frame. 
The general case is characterized by the fact that the G-angular momentum is 
not aligned with any of the body-fixed axes so that the Kasner axes are distinct 
from the body-fixed axes. 

For the class B case, certain cases are characterized by the vanishing of the 
nonpotential force, so that the system becomes an ordinary Hamiltonian sys- 
tem. The symmetric case submanifold A4 5(3) is described by the conditions 
6^ = 9^ = in terms of the combined parametrization (3.1) and (3.7). The 
restriction of W'^ to this submanifold is simply L/on= S^^dO^, using an 
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General case 



Symmetric cases: g G A^s(a) , f 6 = f c = 
class A: all types (a = 1, 2, 3) 
class B: all types (o = 3) 



T 



Generalized Taub-like cases: g G A^GT(a) ; "^6 
class A: I (a = 1, 2, 3), V (a = 3) 
class B: V, VI/^ (a = 3) 



Diagonal case: g e AId 
class A: all types; Va = 
class B: type V; t;i = W2 = 



Taublike symmetric cases: g G A^TS(a) > Vb = Vc 
class A: type VIq (a = 3) , type I (a = 1, 2, 3) 
class B: types VI/i,V (a = 3) 







T 



Taublike cases: g € A^T(a) 

class A- I^^P'''' IX, I (a =1,2,3), 
Class A. {^^^^^^ ^ I VI 

class B: P'^J^L't// ^ -1 = - " " } 
Hype in=VI_i {a = 2); vi = V2 = V3 = 0^ 



= V2 = V3 = 
V2=0 



Isotropic case: g G Mi, Va = 
class A: types IX, VIIq, I 
class B: types VII/j, V 



Table 3: Specialization diagram for the perfect fluid dynamics at the canon- 
ical points of Cd , letting (a,b,c) be a cyclic permutation of (1,2,3) when 
appropriate. For noncanonical points of Cd , the submanifolds analogous to 
■^T(a)j -^TS(a)j ^ud A4gt(o) characterized by the condition (3.20) rather 
than p~ = . For noncanonical class A points of , the type VIq Taublike 
symmetric cases occur for the index value coinciding with the index associated 
with the vanishing diagonal component of n and the Taublike cases occur for 
the index a such that n^**^ and n^'^^ are of the same sign or both vanish. The 
manifold analogous to the isotropic submanifold is obtained by applying to Mi 
the diagonal transformation which transforms the canonical point into the given 
noncanonical point of Cd- 
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obvious abbreviated notation for the restriction, so the restriction of the non- 
potential force is 

(3.49) 

For this symmetric case the fluid must satisfy t;i = t;2 = or equivalently = 

V2 = 0. For type V, one has = and Q [^^3^ is proportional to rf/?"*"; when 

V3 = 0, the third supermomentum constraint becomes /3+ = and the constant 
value of Z?"*" may be transformed to zero by the action of an automorphism, 
leaving an ordinary Hamiltonian system in the remaining degrees of freedom. 
(This case is equivalent to the diagonal case.) The Taublike type Vllh case, 
being equivalent to the Taublike type V case, is also Hamiltonian. (^33 = 0.) 
For type VI/i (or type V) the Taublike symmetric case is described by the 
additional condition dGss = (namely (3.18) which is equivalent to n"a = ^^^') 
or e~" n'^^^e'^^^ = — e~"^n^^^e~^^'^ = e~'^^ q, where q is defined by the 
two conditions q^ = — n^^^n^^) and sgng = sgnn^^\ so that ^33 = e~^" q^ 
and the restriction of the nonpotential force to the Taublike symmetric case 
configuration space is 

Q |j,g(3^= -4aV'+^'^''dlnu , Inu = -6/3+ - 2Age-"'6l3 , (3.50) 

while the third supermomentum constraint becomes 

n'i I ^5(3) = 2oe2(/5°+/5+) (In u) ' = 2klv3 . (3.51) 

When ^3 = 0, then this constraint requires one to restrict the Taublike symmet- 
ric case configuration space by the condition du = 0, leading to the vanishing 
of the nonpotential force and an ordinary Hamiltonian system. Setting q = 
gives the type V limit which is again equivalent to the diagonal case. 

The Taublike symmetric case M.st{3) is diagonalized by the following linear 
transformation AB, which transforms the structure constant tensor out of the 
space Cd except for type V [q — 0) 

K3 = e-"' (-n^i^e^i + n^^^e\) , n = diag(n(i) , n^^) , 0) 

K3 = AB/«3(AB)-i = qe.-"' {e^/VS) , n = ABn(AB)^ = qie\ + e\) 
g = (AB)-i^g(AB)-i = e2(/3"eo+/3+e++re_) ^ ^- = ^^--'0^/^ . 

(3.52) 

The generalized Taublike case Mct(3) is mapped onto the symmetric case 
A4s(i) by this linear transformation. The spatial frame e is a Kasner frame 
consisting of eigenvectors of the spatial Ricci curvature. 

When (73 = 0, then u = — 3(/3+ + \V3^~) has a constant value which may 
be transformed to zero by a translation in 6^. For Bianchi type HI = VI_i, 
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A = sgnnS^^ = ±1 holds and u = 5/?^ (A = 1) or u = —6/3^ (A = — 1) 
and one has a Taubhke case A^t(2) (i-e- P2 = 0) ■^t(i) (i-e- /^r = 0) 
respectively when referred to the barred frame. This Taublike case for Bianchi 
type III is exactly the case where the fourth linearly independent Killing vector 
field which exists for all spatial metrics of this type [10] becomes a spacetime 
Killing vector field, leading to local rotational symmetry. Since the (positive 
curvature) Kantowski-Sachs metrics arc related to these by a simple analytic 
continuation in a particular coordinate system, as discussed in appendix D, the 
Kantowski-Sachs case is also Hamiltonian. 

4 Qualitative Considerations 

With the machinery that has been developed in the preceding sections, one can 

qualitatively or numerically study the time evolution of the general spatially 
homogeneous perfect fluid spacetime. Three regimes in this evolution are of 
particular interest, namely the limit — > — 00 towards the initial singularity in 
the class of initially expanding models (po < as /?" ^ —00), the limit /?" ^ 00 
away from this singularity in the same class of models [38] (omitting the Bianchi 
type IX subclass, which reach a maximum of and recollapse) and the small 
anisotropy limit in the subclass of models which arc compatible with spatial 
isotropy [63] The first regime is especially simple to treat qualitatively, the key 
ideas for which will now be sketched. 

Lifshitz and Khalatnikov. later joined by Belinsky, were led to the qualita- 
tive study of spatially homogeneous models by their investigation of the nature 
of a generic initial cosmological singularity. It seems that the approach to the 
singularity in adapted comoving coordinates induces a sort of contraction of 
the Einstein equations which essentially decouples the evolution of the three- 
plus-one variables at different spatial points while mimicking the behavior of a 
particular spatially homogeneous model at any given point. This work is con- 
cisely summarized and updated in a recent article by these authors [64]. At 
about the same time as the main thrust of the BLK work began, Misner ap- 
proached the problem of the dynamics of spatially homogeneous models from 
an entirely different point of view, paving the way for Ryan to study the qual- 
itative behavior of the approach to the initial singularity in general Bianchi 
type IX models. This is the same class of models to which BLK confined most 
of their attention, the relation between the two approaches being described in 
ref.(24). The formalism of the present article allows one to extend this work to 
the general spatially homogeneous model. 

In studying the regime (3° —00, it is convenient to use the variable 
f2 = — which increases towards the singularity. As one approaches the sin- 
gularity, the "anisotropy" increases and the scale over which the motion of the 

variables occurs becomes much larger than the one over which the curvature 
of the valleys and corners associated with the channels of Figure 4 is apparent. 
The various potentials reduce to time dependent exponential or exponentially 
cutoff potentials in the positive or negative /3+-sectors (defined with respect to 
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a shifted origin) which arc relevant to each potential (as indicated in Figures 
8 and 9). These potentials move with either constant or time varying velocity 
on the /?"'"/?" plane in supertime time gauge. Because of the sharp rise of an 
exponential, a given potential has little effect on the motion of the universe 
point until a "collision" or "bounce" occurs, during which the universe point is 
significantly affected by the potential and then returns to a state in which that 
potential has little effect on the motion. Again due to the exponential nature of 
the potentials, essentially only one potential usually affects the universe point 
at any given time. When no potential exerts a noticable effect on the universe 
point, the evolution reduces to the free dynamics whose exact solutions arc the 
Kasner solutions, characterized by straight line motion in the P~^P~ plane with 
unit /3°-velocity (a null curve in the flat 3-dimensional spacetime diag{3,R) ). 

Exact solutions also exist for each of the cases in which only one potential 
is present. These may be viewed as scattering problems and used to relate the 
asymptotic Kasner solutions before and after the collision with that potential. 
When the velocity of a given potential is constant and timelike (/3°-velocity 
less than one), one can always transform to its rest frame where the problem 
reduces to an ordinary 1-dimensional scattering problem in a fixed potential at 
constant energy [17]. If the velocity is constant but null, one can transform 
to null coordinates and solve the problem. Because of the super-Hamiltonian 
constraint, the system trajectory is always a timelike or null curve in diag{3, R) 
except near the point of maximum expansion in Bianchi type IX models where 
the spatial curvature is positive and the trajectory becomes spacelike as it makes 
the transition from expansion to contraction (see Figure 3 of ref.(17)). 

One must also consider the possibility that the system trajectory is affected 
by one of the channels, although this becomes increasingly less likely as the scale 
of the motion increases (simply because the width of the channel is fixed). 
In this case, referred to as the case of small oscillations by BLK and as a mixing 
bounce by Ryan and others, only a qualitative solution exists for the exit of the 
universe point from the channel [15,29], apart from a numerical calculation of 
a particular open channel mixing boimce [5,22]. This case can be complicated 
by the presence of a centrifugal potential in the channel, as discussed by Chitre 
and Matzner [65] and others [40]. 

Af isncr introduced the extremely useful idea of associating a moving "wall" 
with each potential by selecting a particular contour line which marks the point 
at which the potential has a large enough value to significantly affect the motion 
of the universe point. For a given potential and a given moment of time in super- 
time time gauge, one can consider the simplified super-Hamiltonian constraint 
in which only the expansion energy and that potential appear and solve it for 
pwaii^ the value of the particular /3+ or f3^ coordinate on which the potential 
depends at which the constraint is satisfied. This locates the "wall contour" or 
simply "wall" associated with this potential 

This contour is the one at which a turning point of the motion would occur if 
U were time independent, po were constant, all other potentials were zero and 
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the motion were orthogonal to the contour. In fact since the terms omitted in 
the super-Hamiltonian constraint are positive (away from the point of maximum 
expansion in the type IX case), the "turning point" of the orthogonal component 
of the motion must occur at a smaller value of the potential, namely before the 
system point can make contact with the wall contour. (Because of the motion 
of the walls, the "turning point" may only be an actual turning point in the 
rest frame of the potential.) However, this gap between the wall and "turning 
point" is on a scale which is usually much smaller than the scale over which 
the motion is occurring and so provides a useful marker of where a collision 
occurs between the potential and the universe point. 

Since po is a constant far from the wall contour during a Kasner phase, the 
value of the potential at the wall contour remains fixed and the wall contour 
coincides with a given contour of the potential, thus moving with the same ve- 
locity as the potential. However, po may change during the course of a collision, 
a decrease in |po| requiring a decrease in the potential value of the wall contour 
and causing an additional inward motion of the wall contour towards the uni- 
verse point. In supertime time gauge, the equation of motion for |po| = —Po in 
the expanding phase is given by 

\po\- = -{po,H} = dH/d(3° 

= 12[4e*^\v* + 6a'e'0^) + e^^^HuttiS - {v^r'g"''v'A) - Ukpc^^^"] . 

(4.2) 

Using the definitions of Appendix C, the fluid term inside the brackets of eq.(4.2) 
may be expressed in the obviously nonnegative form 

67(//w^)''e3^"(2-7) [2 - 7 + Ij^-^v^Va] • (4.3) 

Thus all of the terms inside the brackets of eq.(4.2) are nonnegative (apart 
from the positive curvature regime at Bianchi type IX points), indicating that 
\po\ increases with supertime and therefore decreases as il — > oo, except in the 
vacuum type I case where \pq\ is a constant. Note that the equation of motion 
for f2 in this time gauge, namely f2 = dfl/dt = po, determines f2-time as a 
function of supertime; these are affinely related only when po is a constant. 

The walls associated with the various gravitational, centrifugal and tilt po- 
tentials and their fi-velocities are as follows [20,21] 

Jj{a) ^ l(^(a)^2g-8(f;+-i/3") 

(/3+)7" = -if7-ilnbo/n('^)| , 
d(/3+)7"/dl7 = -i - id ln\po\/dn 
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Vcosn/ 

Vcosn ^/ 
C=bo/Pa||e-'"''nWn('^)/3|^ 

(d/3-)r"/rff^ = ±(2V3)-i r^]' + J d In \Pa/po\/dn , 

' ^' (4.4) 

.^|po/^'„||(2V3)-ie-"° (^J'J)l 

(/3-)r" = ±(2%/3)-iln.-i , 
d{p-)r"/d^ = ±{2V3)-^d In \Pa/po\/dfl 

24fc/e^(/^°+'^J)|<| 

(/?a+)tT"=J^ + 5liibo'(48fc/|<|)-i], 

dmmi'ld^ = l + \d \nW/\v'^]/d^ 

Note that there arc two wahs symmctricany located about p'^ = /3^q for the 
centrifugal potentials of type («) and {ii). (For the latter case these walls coincide 
when C = 1 and then disappear for ^ > 1.) These pairs are represented in 
Figures 8 and 9 by single dotted lines at /3~ = (/3~q = at the canonical 
points when well defined). More detailed figures indicating these pairs would 
generalize Figure 1 of ref.(24). Although these walls move away from each other 
as Ipo I decreases for fixed Pa , their separation decreases on a scale which expands 
with n to keep up with the other walls. All of the centrifugal walls move in the 
direction of decreasing values of the associated potential as |po I decreases. 

Between collisions the universe point has unit fi-vclocity while the gravi- 
tational walls have f2-velocity ^ and the component of the fi-velocity of the 
centrifugal walls associated with the explicit dependence on is 0, so the 
universe point will overtake and collide with these walls, the length of Kasner 
segments of the /J+/3~ motion increasing with f2 as Q ^ oo. The tilt wall on 
the other hand also moves with unit O- velocity between collisions neglecting the 
additional component of the velocity due to the time dependence of |f^| (and I 
in the class B case) , so the universe point can never catch this wall unless signif- 
icant changes in its motion occur due to the time dependence of Ipol, \v'a\ and I. 
This means that the tilt potential becomes less and less important as f7 ^ oo, 
as discussed in detail by Ryan [20,21]. All of the walls for each canonical point 
of Cd are shown in Figures 8 and 9. 

A bounce from a single curvature potential Ug""^ is described by the exact 
type II diagonal vacuum solution found by Taub [11]; the supertime t is afhnely 
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related to Taub's time function (by the factor of 12). Consider the case a = 1 for 
comparison with the Lifshitz-Khalatnikov notation; the remaining cases follow 
by cyclic permutation. The Hamiltonian in supertime time gauge is expressed 
in terms of the coordinates {l3f,pi±} by 

H = U-Po' + + PI-') + 6(n(i))^e^(''°-^'5^) , (4.5) 

which is the Hamiltonian of a relativistic particle (with timelike momentum 
due to the super-Hamiltonian constraint) in an exponential potential which 
is moving uniformly in the positive /3]^-dircction with constant /3'^-velocity of 
magnitude v = ^. A boost in this direction with 7-factor (1 — w^)^i = 2/\/3 
transforms the inertial coordinates to the rest frame of the potential [16] 

(P°,P+,Po,Pi+) = 3-^(2/?o -/3+,2/3+ -/?o,2po+P+,2p++po) 

g' = e^^ , fi = 2y3diag(-^+ ,/?"+/?++ ,W +Pt - Pi ) (4.6) 

The problem therefore reduces to free motion in the /3°/3]~ plane and an ordinary 
1-dimensional scattering problem in the /3f -direction with positive energy £0 = 
|(Po^ — Pi-^) = |(Pi+)^ where > is a constant. The super-Hamiltonian 
constraint 

(^tf + 12(n(i))2e-*^'5i = (pr+f (4.7) 

may be integrated to obtain t as a function of f]^; inverting the result yields 
the solution 

(0°{t),f3^{t)) = {-p^t + ^o,Pi-i+ iPiM 

e2^^^(*) = 4\/3|n(i)|(p?^)-^cosh4\/3p?^(f- fo) (4-8) 

poit) = (2po - pr+tanh4V3p^+{t- to)) 

where to is the turnaround time. This solution interpolates between two asymp- 
totic Kasner solutions characterized by the null momenta 

(po(±oo),pi+(±oo),pi_) = 3-5(2po TP?°+,±2p?^ -Po,^/3pl-) (4.9) 

The scattering problem is therefore most naturally interpreted as a map of the 
unit circle in the p+P- plane into itself, which is the result of the Lorentz 
transformation of ordinary reflection from the static potential in its rest frame. 

The Lifshitz-Khalatnikov parametrization of the unit circle in the p+P- plane 
leads to the following representation of the unit momentum coordinates 

(pi+,Pi_) = i(l+w + «2)-i((i+„)2,V3(«2-l)) , (4.10) 

which enables one to define u±oo for the asymptotic values (4.11) of this unit 
vector at i = ±00. The final states with Uoo € [l,oo) correspond to initial 
states with u G (—00, —1]; similarly Uoo S (0, 1] corresponds to U-00 € [—1,0). 
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Figure 10: The Lorentz transforms of simple reflection by the gravitational 
potentials U^^^ and f/^^^ in their rest frames. The asymptotic Kasner states are 
related by a change of sign of the parameter u in the first case and of U(2) in 
the second case. 

No asymptotic states exist with u = or u = ±00 which represent the Lorentz 
transforms of directions parallel to the equipotential lines of the static potential. 
Figure 10 illustrates this situation. 
Since 

Pi-/\Po\ = {pi-/\po\)[2-pi+/\po\]-' (4.11) 
is a constant of the motion, its asymptotic values 

^i_(±oo) = pi_(±oo)[2 -pi+(±oo)]-i = 3-5(1 - u^^) (4.12) 

must coincide, which can only be true if u_oo = —Uao- (Woo = w_oo is ruled out 
by the case pi_ = of normal incidence.) Thus scattering from the curvature 
wall Ug^^ simply leads to a change of sign of the Lifshitz-Khalatnikov parameter 
= u. Similarly scattering from the rotated potentials Ug^^ and Ug^^ leads to 
a change of sign of the rotated Lifshitz-Khalatnikov parameters U(2) = P231 (u) 
and U(3) = P3i2(w). 

Suppose Uoo € [2,00) and hence U-00 € (—00,-2], so that the Kasner ex- 
ponents at t = 00 are ordered: pi{uoo) < P2{u<x) < P3(W(x>)- The transposi- 
tion P12 then orders the Kasner exponents at i = —00 by reflecting the sector 
u e (—00, —2] across the p+-axis into the ordered sector 

Pl2{U-oo) = -(1 + U-00) = Uoo - 1 

(Pl(Woo - 1),P2(u2(Moo - l),P3(Uoo - 1) = (P2(U-oo),Pl(U-oo),P3(U-oo) 
Pl{Uoo - 1) < P2{Uoo - 1) < PaiUoo - 1) • 

(4.13) 

Thus the "ordered" Lifshitz-Khalatnikov parameter decreases by 1 in the neg- 
ative time direction, i.e. as increases. On the other hand if u^o G [1,2), so 
that U-00 € (—2,-1] and the Kasner parameters at ? = 00 are ordered, then 
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the cyclic permutation P312 orders the Kasner exponents at f = 00 by rotating 
the sector u e [—2, —1] into the ordered sector 

-P312(M-oo) = -(1 + U_oo)"^ = (Uoo - 1)"^ 

1 1 , (4-14) 

Piiiuoo - ir^) < P2iiuoo - < P3{{uoo - I'r') ■ 

The "ordered" Lifshitz-Khalatnikov parameter thus decreases by 1, entering the 
interval [0,1) and then inverts. 

The connection between a Lie algebra contraction and an anisotropic sin- 
gularity may be illustrated by introducing the two special turnaround times 
defined by = ±(4\/3)-^ ln(\/3|n(^VPi+)- For to = i^, the future Kasner 
limit P^{t) ^Pi+t is obtained as t ^ 00; the same asymptotic behavior 
arises from the contraction n^^^ ^0 which sends the turnaround time to —00. 
In other words long after the collision the solution approaches the solution for 
the contracted Bianc;hi type corresponding to the anisotropic singular scaling 
/3(f) 0{oo). Similarly for to = , the past Kasner limit /3f (f) — * Pi+t is 
obtained either as i — > —00 or n^^^ 0. 

A collision with a single centrifugal potential J/i"^ is described by an exact 
Bianchi type I vacuum solution with nonzero angular momentum, namely the 
symmetric case A4s{a)- The solution may be found either by transforming the 
Kasner solutions or by solving the equations of motion. For definiteness consider 
the case a = 3. The Hamiltonian in supertime time gauge is 

H = i(-po' + Sa-'^'Ps' (4.15) 

An S'-'^-parametrized family of potentials arc possible, corresponding to the pos- 
sible values of K3 given by (3.4) with </> G [0, 2tt). P3 is a constant of the motion 
and since H only depends on /3~ , both /3° and /3+ obey the free dynamics with 
Po and p+ constants of the motion, reducing the problem to a 1-dimcnsional 
scattering problem for with constant energy £0 = ^{Po^ —P+'^) = {p^}'^ and 

(3) 

static potential Uc ■ The solution may be obtained by integrating the super- 
Hamiltonian constraint for the supertime as a function of (3~ and then inverting 
the result (note that supertime and 0-time are affinely related here). 
The following variables make this calculation simple 

v = 2V^(3-, =2-i(cos</)e''±sin(?!.e-"), 
(m^)^ — (w^)^ = sin20 , du^=u~dv, 

\Q33 = {u-f = {u+f - sin2</,, So = ^(p-f + iPsH^-)-^ , 

cosh~"^|u+/6| e = 1 
sinh-^(u+/6) e = -l 
. .it+l e = 

eb^ = sin 2^ + fPs V^o , 2^3/^0 = 5 In | tan (j)} . 

(4.16) 

The solutions e = and e = — 1 occur only for the type (m) potential when 
the energy equals or is greater than, respectively, the maximum value of the 



{24£o)Hi-to) = Jdu+[{u+f - eb^]"* = < 
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potential Umax = f-Ps^l sin2(/)|~-'^ for this case; the latter solutions pass over the 
potential barrier while the former ones approach or leave the point of unstable 
equilibrium /3~ = 0o ■ The e = 1 bounce solutions are then 

\u+\=bcosh{24£o)Ht-to) p- ^ ±\2£o\^l - iPs/Soiu-)-^] 

+ (2V3)-icosh"V^|sin20|-5| sin2</)>0 
= S + (2V3)-^sinh-^(M+|sin20|-5) sm2(j) < 
.(2\/3)-ncos2^)ln|\/2w+| sin2(;i = 

Note that for a fixed value of P3, once £0 decreases below Umax, only the bounce 
solutions are relevant, i.e. for small enough £0, all three types of potential reverse 

the orthogonal component of the motion. Similarly as £q decreases, the bounce 
occurs far enough from (3^ that the potentials (i) and (ii) essentially reduce to 
pure exponentials. 

These solutions have Kasner asymptotes as t ^ ±00 which are characterized 
by the momenta (po(±oo), p+(±oo), p_(±oo)) = (po, P+, ±p5??) for a bounce 
solution and (po, P+, P^) for an e = —1 solution, but (p_(oo),p_(— 00)) = 
{0,p°^) or (p^,0) for an e = solution. The change in asymptotic momenta 
associated with the bounce solutions corresponds to reflection across the p_- 
axis. 

For this symmetric case one has {9^, 9"^, 9^) = (0,0,9) and W°' = d'^39 so 
from (3.21) one has 

9 = 6P3 J dig-^^^ , z = Ctanh(24f)5(f- l^) , ^ = sin2</.|(P3)"^)^ , 

tan~^ z sin 20 > 

tanh"^^; sin 20 < 

z sin 20 = 

(4._18) 

In the compact case sin 20 > 0, let 9 = |sin20|2^; since 6^"=* = — 1 for ^ = 

tt), 9 E [O.tt) represents a closed path in Ais(3), although twice this interval 
represents a closed path in G. Note that for (:$>!, tan2(^ — ^0) begins at 
a very large negative value a.t t — to ^ —00 and approaches the sign reversed 
value at f — to ~^ 00, so that the total change in 29 during one bounce is just 
slightly less that tt. As ^ decreases the total change in 9 decreases from tt/2 and 
approaches as C ^ 0. This is clearly evident in the numerical calculation for 
the canonical type IX symmetric case A4s{3) [22], where P3 is constant but \pq\ 
and therefore £q decreases as 17 ^ oc. In all cases {A9)bounce ^ as C ^ 0. 
The rescaled DeWitt metric on Ms{:i) (set {9'^,9'^,9^) = (0,0,61)) is 

-dp° d(3° + dp+ (g> d0+ + '^i . '^h = dP" (g) dp- + i^asd^ d9 . (4.19) 
The affinely parametrized null geodesies of this metric describe the solutions of 

(3) 

tlie collision with the single potential Uc in supertime time gauge. The P 



9-9o = 3P3 j rft[(u+)^ -sin20]-i = i|sin20|-5 < 
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scattering problem is then equivalent to finding the geodesies of ^h. This is just 
the induced metric on the future hypcrboloid 

H^2V-s)-^ = {{y°,y\y') e \ {y'f-{y'f-{y'f = (2^3)-^ > o} (4.20) 

in 3-dimensional Minskowski spacetime with metric —dif (g) dy^ + dy^ (g) dy^ + 
dy'^^dy'^. Defining 

maps the coordinate surfaces of constant {l3^,/3~^) onto H^^^^i. Letting x ^ 6 
in {AA) and letting (03(2), 53(2)) be obtained from (03,53) by doubling the 
argument, and using the double angle identities 

cs' = i(c3(2) + 1) .3' = -Urh^'Y'{c3{2) - 1) , (4.22) 

one obtains explicitly for the case m"^ = (n(i)n(2)))5 ^ (set e = sgn(n(i)n(2))) 

yO +2^1 = e^^^o'dfe"-"" - ee-(''-'"')]c3(2) + i[e''-''° + ee-^^""")]) , 

yO - yi = ee^^'^o'd [e"-"" - ee-(''-^°']c3(2) + ie"""" + ee-^"-"")]) , 

y^ = sgnn(i)(i[e''-^'' - ee(''-''°'])s3(2) ^; - t;o = 2\/3(/3- - /Jq") • 

(4.23) 

The special linear subgroup SL{2)s acting on Ms{3) maps onto the action of the 
3-dimensional Lorentz group on H^^^-^ _ 1 ; the coordinates ,6} are comoving 
with respect to the subgroup generated by /C3. This subgroup corresponds to 
null rotations for (/) = 0,7r, rotations in the y^ — y^ plane for </> = 7r/4, 57r/4, 
boosts along for ^ = Z-K/i:,7'K/A and combinations of these for other values 
(e_ generates boosts along y^). 

Geodesies with 9 = 0q have zero angular momentum with respect to this 
subgroup, i.e. P3 = 0, and are diagonalized by the constant active transfor- 
mation e^°'^^. However, these geodesies have nonzero angular momentum with 
respect to the subgroups generated by other values of K3 (modulo sign). In the 
canonical type IX or VIII cases, 6 ^ + O0 corresponds to a rotation or boost 
by the angle or natural boost parameter 2610 because of the double angle func- 
tions. For small angular momentum in the canonical type IX geodesic 
has A{26) tt as if the projection to the y^ — y^ plane were flat, but as the 
angular momentum increases, the curvature of the hyperboloid decreases this 
angle. 

The classic discussion of Lifshitz and Khalatnikov of an "era" of successive 
alternating Kasner regimes describes a sequence of bounces between the grav- 
itational wall J/'^-* in the upper half plane and U^^^ in the lower half plane as 
the universe moves out of the corner where these walls intersect with increasing 
supertime. Reversing the direction of time in order to approach the initial sin- 
gularity, as in f2-time, one must reverse the sign of the momenta; "initial" and 
"flnal" will be correlated with this reversed direction of time in the following 
discussion. An initial Kasner phase with u°^'^'^'^'^'^ = u e [1, 00) is headed for the 
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wall /?+ = and then bounces from this wall into the second Kasner 

phase with = u - 1 so that it is headed for the wall = (^^ )^''". It 

then bounces from this wall into the Kasner phase with u"^'^'^^'^'^ = u — 2, etc. 
When u°^'^''^'^'^ enters the interval [0,1], the final Kasner phase (or "epoch") 
heads away from both walls. Figure 10 indicates the scattering problem for 
each bounce. 

This approximation is valid as long as the universe point does not penetrate 
deep enough into the vertex or corner between the two asymptotic potentials to 
feel the effects of the channel where the straight line contour approximation of 
Figure 5 breaks down. When the effect of the channel is important, the regime 
is called a "long era" or "mixing bounce" and the type VII or VI gravitational 
potential for an open or closed channel respectively may be used as an approxi- 
mation for the channels of types VIII and IX outside their common intersection. 
A qualitative treatment of the exact type VIIq and VIq diagonal vacuum case 
in this context was first given by Lifshitz and Khalatnikov [26] and then in a 
better approximation by Khalatnikov and Pokrovsky [29]. 

The Hamiltonian for this case in supertime time gauge at the canonical 
points of is 



VIIo 
VIo 



(4.24) 

The potential moves with constant velocity dp~^/d^ = — 1 (a null velocity 
in diag{3,R)) so a rest frame does not exist and one must instead use null 
coordinates 

iw,v,p^,Py) = (/3° + /3+,/3°-/3+,i(po+P+),|(Po-p+)) , 

iP°,P'^,Pa,P+) = {■^{w + v), ^{w - v),p^ +Pv,Pw -Pv) , 

g = diag(e2-, e^™, e-'"+3'') , g ^ /"^ 3, (g) : {w, v)^{w + a,v+ ^a) , 



e" 



Vcosh2 2V3/3-/ ^ 
w — wq = —2pyt , pt, = , Pt, < (expansion) , 



V = -2p^ = 8 • 246 



4w 



/sinh22\/3/3- 
Vcosh22V3/3- 



(4.25) 

The coordinate w comoves with the potential function, while v describes the 
motion of the system relative to the potential. 

The Taubliko typo VIo case with = = p_ is an exact solution which is 
useful to consider in detail and represents the path at the center of the closed 
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channel 

v-vo-Vt = ipy-^e^"" = ik^^ ' Q = H = -2p^V -^V = , 

= i^w + Ip.-^e^ , - y.-^e^) , d/3+/df](±oo) = t1 , 
(Po,P+) = (-KI(l + 12p.-V'"),-|p,|(l-12p,-2e4'")) . 

(4.26) 

The turnaround time is defined by Po(fo) = or 12p^~^e*"'^*°^ = 1, while 
setting ipo^ 24e'*(^°+'^iaii) = 24e'*'"™<"" defines the location of the wall for 
this /3~ = case 

PLa = ^=1 ln[bo|(4y3)-i] , dpl^,,dn = 1 - 4(1 + {UyWe^"")-^ , 

(4.27) 

The system has the past Kasner limit a.t t — —00 — w characterized by the 
Lifshitz-Khalatnikov parameter u = 5 (free motion in the positive /3+ -direction 
towards the wall), and the Kasner limit at f = 00 = w characterized by the L-K 
parameter u = 00 (free motion away from the wall) long after the collision when 
it has reversed direction and is again moving along a null curve in diag{3,R). 
The future Kasner limit does not look like free motion in the supertime time 
gauge since the supertime remains sensitive to the small curvature in this limit 
and does not reach a stage where it is affinely related to fl-time and instead 
diverges from f2-time. 

From the point of view of fi-time, both the universe point and the wall ini- 
tially have unit f2-velocity and hence the universe point cannot overtake the 
wall. (This is the origin of the term "long era'' , since this behavior will persist 
for a long period of O-time.) However, since the universe point and potential 
have no relative motion, the very small value of the potential at the location of 
the universe point has a continued effect and over a long time begins to reverse 
the motion and decrease the expansion energy which in turn slows down the wall 
enough so that an actual collision occurs at the turnaround point when both 
the wall and universe point have zero f2-velocity and the same position, after 
which the wall resumes its unit fi-vclocity in the same direction and the uni- 
verse point returns to its unit f2-velocity but in the opposite direction. Modulo 
constants this same solution describes a head on collision with the approximate 
tilt potential U^-H for constant \v'^\. 

By introducing the new time variable ^ = |6/p„|e^^ = |6/p„|e^"'°"'"^l^'"l* and 
the abbreviation q — 4v^/3~ , the equation of motion for (3~ is given by [29] 

d'^q/d^^ + C^dq/di + svahq = . (4.28) 

For (7^1 where sinhg w q, this is just Bcsscl's equation of order zero [16,26], 
while for ^ ^ 1 (at very late supertime), this has an exact solution in terms 
of Jacobi elliptic functions [29] . Khalatikov and Pokrovsky use a WKB-like ap- 
proximation involving this approximate solution to obtain a qualitative solution 
valid for all times and whose limit at t ^ —00 specifies the asymptotic Kasner 



61 



phase which precedes the channel corner run in supertime (but follows it in 
f2-time) . 

Consider the case /3~ <^ 1 where the approximation sinh4\/3/3~ 4\/3/3~ 

is valid and the part of the Haniiltonian which governs the /3~-motion is just 
that of a harmonic oscillator with the supertime frequency u = 246^^" = 4|pt,|^ 

H = -p^p, + 24e''^6^jjj +H_ H_ = \pJ + \uj\l3-f . (4.29) 

In the adiabatic limit which occurs at late supertime (^ ^ 1) when the change 
in frequency per period of ocsillation goes to zero, one has the approximate 
solution 

B~ = Acosi I ojdt) = Acos(£ - £o) , 

J (4.30) 

A = (6e-2-/12)^ = \2p,^/b\-^ = [Stt^/C]-^ , 

where the expression for the amplitude follows from the constancy of the adia- 
batic invariant H-/u! = \A^uj = b, as discussed by Misner [15]. The constant 
C is the constant introduced by Khalatnikov and Pokrovsky [29] . 

The equation of motion for v may be averaged, using ((/3^)^) = ■^A'^ 

V = 192e^'"S^jjj + {A8fe^^{f3-f , (v) = 1926^"" 6^ jjj + 96fee^^ , 

{v) = 2A\p,\-^e*^S^jjj + Ube^^ + V , (4.31) 

{v - vo) + Vt = 3\p.\-\^^5^jjj + Zbe^-^ = \b\p,\£, + ^es^iii ■ 

The averaged super-Hamiltonian constraint requires V = Q. Note that the 
leading approximation to the (ui, v) motion in the closed channel case is just 
the Taublike type VIq vacuum solution. 

At early f2-time, the universe point in its initial asymptotic Kasner phase 
and the potential both have unit 0-velocity and keep pace with each other, so 
w = (3~^ — il changes only slowly in fi-time. (Supertime and fi-time diverge 
at t ^ 00.) As the frequency decreases with increasing f2-time, the amplitude 
increases and the (3'^ motion slows down relative to the potential until the small 
/? approximation breaks down, eventually reversing its direction, and after a 
final bounce goes into an asymptotic Kasner phase which leaves the channel 
potential behind. 

One can also consider a mixing bounce with a centrifugal wall present in 
the channel. This corresponds to the Bianchi type VIIq or VIq symmetric case 
MscA)- Adding the term W^^Q'^^^ to (4.24) only effects the ^~ motion di- 
rectly, while the constant P3 enables one to find the variable defined as in 
(4.18). Matzner and Chitre [65] discuss this case in the adiabatic limit for 
/3" ^ 1. The mixing bounce is confined to the region between the channel 
gravitational wall and one centrifugal wall in the open channel case where the 
(3~ motion corresponds to a 2-dimensional oscillator with nonzero angular mo- 
mentum. In the closed channel case the situation is much different. 

When the curvature of the channel is not important, one simply has a se- 
quence of bounces between one gravitational wall and a centrifugal wall. During 
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the bounces from the latter waU the ordered Lifshitz-Khalatnikov parameter 
^ordered ^jggg change. Such a sequence of bounces with an exponential cen- 
trifugal wall is described in terms of transformations of the first and second type 
in rcf.(24). 

The importance of the scalar curvature term on the dynamics also can 
be described by a wall 

^classB ^ 72a2e4(/3°+/3+) , 

» ' (4.32) 

PLu = + I In bo/12a| d^+^Jdn = 1 + irfln \po\/dn . 

This potential becomes important only during collisions with its associated wall, 
but since this wall has unit fi-velocity away from the location of the universe 

point when it is in free motion, it becomes less and less important as f2 ^ oo. 
The exact diagonal Bianchi type V solution describes a collision with this wall 

H = i(-po' + + P-^) + i(12a)2e^('5°+'5+) . (4.33) 

obey the free equations of motion (since V* = 0) with constant p±, leaving 
the super-Hamiltonian constraint to determine /J". This is easily integrated if 
one imposes the supermomentum constraint p+ — 0, thus setting (3^ = f]^ . 
(The equations of motion for 7^ are integrable in null coordinates, but the 
super-Hamiltonian constraint then forces p+ = 0.) 

r = p.{i- to) + 00 , - (12a)^e^('^°+'^o+) = pJ , 

e2(/3°+/3o+) = -|p_/12a|csch|p_|(f - to) , -po = -b-|coth2|p_|(i - to) , 

PLii =P^ + 1 In cosh 2\p_ I {t - to) , dp+^Jdn = tanh22|p_ | {t -to). 

(4.34) 

The evolution of corresponds to motion in the 1-dimensional potential 
— i(12a)2e''^'' ^ with constant energy £0 = \p-'^- The expanding solution 
has t — to e (—00,0); at i — > —00 the solution has a Kasner limit but it runs 
out of supertime at f = fo when a collision takes place. The wall moves in from 
/3+ = 00 at t = — oc' when it has unit O- velocity to = (i^ where its velocity 
vanishes and the wall makes contact with the universe point at a finite value of 
However, the proper time goes to 00 at F = fo since near t'^io 

t-to = j e^^''rfi~ j \i-io\~^dt<^\i-io\~^ ■ (4.35) 

Thus even though (3^ obeys the free dynamics in supertime time gauge, it runs 
out of supertime at a finite value 0q , i.e. its value becomes frozen in proper 
time or il-time. 

It is basically an effect of this kind which seems to result in isotropization of 
type VII/i models with "frozen in" values of /3~ as /?" 00 [40] . A similar effect 
occurs for all Bianchi types when the matter contributions to the Hamiltonian 
dominate the dynamics as /3° — > 00, as noted by Doroshkevich, Lukash and 
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Novikov [40]. For example, consider the exact diagonal type I case where Va = 

= Pa and the Hamiltonian is 

H = li-po' + p+2 + + 12e^0''HM , 12e^^°HM = 2AW e^^" ^^-'^ . 

(4.36) 

Here p± are constants and /J^ undergo free motion, reducing problem of the 

evolution of to motion in the potential —24,kPe^^^ (2-7) with constant energy 
£q — One may directly integrate the super-Hamiltonian constraint 



r|(24fcr/£:o)*sinh(2£:o)^(io-i)|-'('-^^"/' 7 7^2 (4.37) 

|^g[2(£o+24feP)li(t-to) 7 = 2 ' 

The growth of the matter Hamiltonian as /3° ^ 00 again leads to a singularity 
in supertime and hence a "freezing in" of the values of occurs in proper time 
or 0-time when 7^2. However, these effects are important for the direction 
/3° — ^ 00 but diminish in importance as O — > 00, since the solution quickly 
approaches its Kasner asymptote. 

On the other hand when 7 = 2, the system again experiences free motion 
but along a timelike direction in diag{3, R) 

V'^^PAPB = -48kl^ , d/3/dn = [(d(3+ /dnf + {df3- /dnf]^ = [l-48A;;V-2]^ . 

(4.38) 

A stiff perfect fluid with Va = only changes the supertime time gauge Hamil- 
tonian by the addition of a positive constant; the only change one must make 
in the description of a collision with a gravitational or centrifugal wall is £0 
£0 + 24kP, where £0 is the constant energy defined in each of these problems. 
(It is exactly this change which allows one to easily insert a stiff perfect fluid 
in the Taub-Nut spacetime [41].) Collisions with gravitational walls as ^ 00 
decrease |po | and therefore reduce the f2- velocity of the universe point between 
collisions. If the universe point does not escape from the gravitational walls 
first, then eventually one has d(3/dVl < ^ . after which the universe point can no 
longer catch the faster moving gravitational walls and remains in this free mo- 
tion phase as Q — » 00. Thus such a stiff perfect fluid (or scalar field [66]) leads 
to an end of an "oscillatory approach to the initial singularity". (See below.) 
Similar calculations may be used to study the effect of a nonzero cosmological 
constant which only directly affects the equation of motion of (3^ and then in- 
directly affects the remaining variables through the appearance of (3'^ in their 
equations of motion [67]. 

What has been described so far are exact or qualitative solutions for the case 
in which only one or two potentials in the supertime time gauge Hamiltonian 
play an important role at any given time. The importance of a given potential is 
described by a collision of the universe point with the wall associated with that 
potential or with a corner where two walls meet and a more careful description 
may be required. Each such phase of the motion is preceded and followed 
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by asymptotic Kasncr phases which arc related to each other by the exact or 
quahtative solution for that collision. Each collision is therefore interpreted as 
a scattering problem. As f2 — > oo one quickly reaches a stage in which the 
description in terms of successive collisions with different walls is valid and one 
may approximate the evolution by matching together the Kasner asymptotes 
between collisions. 

The class A diagonal case is easiest to describe in this manner. As (4.36) 
shows, the matter potential is exponentially cutoff as ^ oo and even when it 
dominates the spatial curvature, (4.37) shows that it has a very shortlived effect 
as n ^ oo, i.e. "the matter is dynamically unimportant" as ^ oo. In this 
limit the system obeys the vacuum dynamics, leaving only the gravitational walls 
to determine the qualitative solution. Since these walls move with f2-velocity 
i between collisions while the universe point moves with unit fJ-vclocity, in 
the semisimple case of Bianchi types VIII and IX where these walls form closed 
trapping regions, the collisions never stop, leading to an "oscillatory approach to 
the initial singularity" in the words of Lifshitz, Khalatnikov and Bclinsky. The 
oscillations occur in the values of the variables , and e^' , a half period of 
an oscillation corresponding to two Kasner phases ("Kasner epochs" [24]) joined 
by a collision at the maximum or minimum point. A single excursion of the 
universe point into the vertex between two walls is called a Kasner era. Except 
for collisions with corners between two walls, these two types experience the 
same motion. However, as ^ oo the size of the corner region where channel 
effects are important becomes increasingly smaller in relation to the triangular 
trapping region. For types VIq and VIIq the missing third wall provides an 
escape route which leads to only a single Kasner era before the system enters a 
final Kasner phase as f2 ^ cx), while for type II only one collision with the single 
gravitational wall occurs between two asymptotic Kasner epochs. The Taublike 
cases for types VIq, VIII and IX are special cases characterized by channel 
trajectories which undergo a single collision with one gravitational wall followed 
by an escape down an open channel in the semisimple case, while the isotropic 
cases have no motion at all. The type IX isotropic case is exceptional in that 
the matter is always dynamically important since the positive spatial curvature 
makes the vacuum super-Hamiltonian negative-definite. The Kasner axes for 
the diagonal case coincide with the body-fixed and space-fixed axes and are 
eigendirections of the spatial Ricci curvature tensor. 

The symmetric cases ( "nontumbling models" [68]] are much easier to de- 
scribe qualitatively than the general case since the G-angular momentum Pa 
is aligned with one of the diagonal axes which is in fact an eigendirection of 
the spatial Ricci curvature in all instances. This axis is also an eigendirection 
of both g and K, so it is a Kasner axis. Only the 1-parameter subgroup of G 
which leaves this preferred axis invariant is relevant to the dynamics; for this 
subgroup Pa and Pa coincide and have a single nonvanishing component. The 
same is true of v'^ and Va and of the components of the supermomentum P' {8a) 
and P{Sa) with respect to the primed and unprimed frames, all of which have 
a single (possibly) nonvanishing component along the preferred axis. For the 
symmetric case A4s{a)j the index a is associated with this preferred axis; as 



65 



usual let (a, 6, c) be the associated cyclic permutation of (1,2,3). (Because of 
the block diagonal form of g and K, the remaining two Kasner axes are linear 
combinations of and Bc-) Then Va = is a constant. In the class A case 
excluding the degenerate types I and II, a may assume any value and I and 
Pa = Pa are constants such that e" Pa = —2klv'g^, leading to a single nonzero 
centrifugal potential which is static in supertime time gauge. In the class B 
case, only a = 3 is allowed and e" P3 + ap+ = —2klv'^, but neither P3 nor I is 
a constant and the centrifugal potential is not static. 

For the case n^"^ ^ and rank n > 1, the walls associated with the grav- 
itational potentials and the single centrifugal potential define closed trapping 
regions as is easily seen from Figure 9. For the type {ii) centrifugal potential, 
the pair of walls do not come into existence until C <1 (see (4.3)) but since |po| 
decreases towards zero as $1 ^ 00, this condition is eventually satisfied (unless 
|po I has a limit away from zero as may occur in the case of a stiff perfect fluid) . 
In the class A case one may describe the Q — > 00 evolution in terms of collisions 
with the centrifugal wall and the gravitational walls between regimes of free 
motion. Since these walls move with fi-velocity and ^ respectively between 
collisions, while the universe point moves with unit f2-velocity, in the case where 
closed trapping regions exist the collisions never stop, leading to an "oscillatory 
approach to the initial singularity" . In this case one must make a distinction 
between a Kasner epoch and an ordinary epoch. An ordinary epoch refers to a 
period of free motion in the plane, while a Kasner epoch refers to a period 

of free motion with respect to the full rescaled DeWitt metric which includes 
collisions with centrifugal walls. In the latter case one can always transform 
away the constant G-angular momentum Pa by a constant linear transforma- 
tion from the space-fixed axes and Cc to the Kasner axes. Since the Kasner 
axes are time independent in a period of free motion, they can change only 
during collisions with gravitational walls. 

As |po| decreases due to collisions with gravitational walls, the single cen- 
trifugal wall moves outward toward decreasing values of the centrifugal potential 
where the type (i) and (ii) potentials rapidly approach the asymptotic expo- 
nential potentials of type (Hi) and the distinction between the three types of 
centrifugal potentials becomes unimportant. Note also that the value of Q-^°-'^ 
at the universe point experiences local maxima at each collision with the cor- 
responding centrifugal wall and these local maxima approach zero as the wall 
moves out. Since 0-1"" goes to zero and since Pa is a constant in the class A 
symmetric case A^d, the offdiagonal velocity W^" = QQ~^°'°'Pa also goes to zero. 
The single nonvanishing automorphism variable 9 is given by the integral of 
W'^. Because 0-i«" is exponentially cutoff, 9 essentially changes only during a 
collision with the centrifugal wall where the value of g-^°-°- at the universe point 
assumes a local maximum which decreases with each collision. The discussion 
of (4.18) in fact shows that the total change {l^9)hounce of this variable during 
a single collision approaches zero as |po| — 0. An argument due to Lifshitz, 
Khalatnikov and Belinsky [26] shows that the sum of these changes from any 
given time as ^ 00 is finite, so 6 approaches a limiting value as il — » 00. 

If rS"-^ = 0, then the universe point bounces around between the centrifugal 
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wall and a gravitational wall until it eventually enters a final Kasner phase which 
leaves both walls behind forever as ^ oo. In the type VIII case with a type 
(ii) centrifugal potential, the universe point will collide with the gravitational 
walls until |po| is lowered enough for the pair of centrifugal walls to come into 
existence, after which the motion will be confined to a smaller trapping region. 
Similarly in the type VIq symmetric case with a type (ii) centrifugal potential 
(associated with A4s(3) in the canonical case), it is possible that a final Kasner 
phase is reached before the centrifugal walls appear. 

Consider the analogy with the central force problem. If one imagines a 
nonrelativistic particle in free motion on whose radial motion is continu- 
ally reversed at large distances from the origin by collisions with an expanding 
spherical shell centered at the origin, the angular momentum will be conserved 
but the impact vector connecting the origin to the point of closest approach of 
the particle to the origin will change at each collision. A different translation 
of the origin will then be required in each free motion phase between collisions 
to transform away the angular momentum. The Kasner axes are changed by 
collisions with the gravitational walls for a similar reason. Thus even after the 
motion in the offdiagonal variable 6 essentially stops and one transforms to 
a diagonal gauge frame where its limiting value 0q is zero (so that = ejj), 
two of the Kasner axes do not coincide with the diagonal axes (otherwise the 
G-angular momentum would be zero) and continue to change at each collision 
with a gravitational wall. 

As long as the a^-curvature potential is important in the class B case, this 
description is not valid since the system will not experience free motion between 
collisions with the walls and the analytic description of an individual collision 
given above will not be valid. However, eventually this curvature potential 
becomes unimportant as O ^ oo as indicated by the above wall description 
and the evolution reduces to that of the corresponding class A model with one 
important exception. Only the symmetric case A^s(3) is allowed and it is not 
the G-angular momentum P3 = P3 which is a constant for the vacuum dynamics 
but rather the supermomentum e" P3 + ap+. The equation of motion for P3 
in this case, namely (P3)' = NQ^ = 48ae^'^^ ^e" ^33, is exactly — ae~" 
times the equation of motion for p+ (excluding the degenerate type V case), 
so even in the limit that the curvature potential is unimportant, P3 still 
changes during collisions with gravitational walls to compensate for the change 
in p+ which occurs during these collisions, causing additional motion of the 
centrifugal wall. This additional change in the G-angular momentum compared 
to the corresponding class A case then leads to an additional change in the 
Kasner axes relative to that case. Also since n^^^ = for the class B case, 
closed trapping zones do not occur and a final Kasner phase exists as O — > 00 
when the universe eventually leaves the walls behind. 

This discussion has assumed that the matter is not dynamically important 
as in the diagonal case (when 7 7^ 2). The fluid variable ri may be expressed 
in the form n = e^^l/u^ where > 1. In the class A case I is a constant; for 
given values of I and O, n can only be less than its corresponding value in the 
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diagonal case where = 1 since ?;„ 7^ implies > 1. A similar statement 
holds for ix and p. This means that the only difference the nonzero value of Va can 
make that might lead to the dynamical importance of matter in the symmetric 
case M-sCi) would be if m"*- ^ 1. In this case the matter super-Hamiltonian is 

approximated by the exponential tilt potential U^^i^, but since is a constant, 
the associated tilt wall moves out with unit fi-velocity between collisions of 
the universe point with the gravitational walls and becomes increasingly less 
important as ^1 ^ 00. In this limit one is justified in ignoring the matter 
potentials. However, the constant nonzero value of the fluid supermomentum 
is required to allow a nonvanishing centrifugal potential (when the symmetric 
case is nontrivial, i.e. not equivalent to a diagonal case). 

In the class B case, I is not constant but satisfies the equation of motion 
(InZ)' = 2Aae^ '^^^ vs/v-^, so this argument does not go through. Because 
of the supermomentum constraint, I must approach a constant as ^1 ^ 00 if 
the gravitational supermomentum is to approach a constant as assumed above. 
In fact since the equation of motion for the supermomentum constraint is 
{2klv^)' = (e" + ap+)' = a{p+,H^*}, / is negligible only when the tilt po- 
tential has a negligible effect on the jS'^ motion. Apparently such a limiting 
stage is reached according to the work of Percsetsky [35] . 

In the general case ("tumbling models" [68]) the angular momentum is not 
aligned with a diagonal axis and the time dependence of Pa complicates matters 
considerably by making the three centrifugal potentials implicitly time depen- 
dent. Although one or two components of Pa may vanish at one instant of time, 
this can only occur at isolated times without reducing the problem to a symmet- 
ric case, so except for these exceptional moments of time, all three centrifugal 
potentials are present and the varying values of the angular momenta induce 
time dependent translations of each of these potentials. For the class A case, 
e°'"+°'''n"-''PaPb = e"-''+°'''n°'^PaPb is a constant of the motion [43], so only two 
degrees of freedom are involved in translating the centrifugal potentials. Closed 
trapping regions exist for all Bianchi types but I and V in the general case, 
leading to a never ending sequence of bounces from the various potentials, or 
an "oscillatory approach to the initial singularity" . 

Examining Figure 9, one sees that the closed trapping regions are all bor- 
dered by one gravitational wall and cither two or three centrifugal potentials, 
depending on the relative magnitudes of the components of the offdiagonal mo- 
menta. Because the walls associated with ui^^ and ul?^ have been drawn in 
symmetrical positions for the canonical types VI, VII and II, only two bounding 
centrifugal walls are shown for these types, but in general the trapping region 
situation is described by the type IV figure with three bounding centrifugal 
walls. This is true even for the types VIII and IX where a more detailed figure 
showing pairs of walls at varying distances from the /?+-axcs would make this 
possibility obvious. However, since the gravitational walls move with il-velocity 
\ between collisions, while the motion of the centrifugal walls is only due to the 
time dependence of the offdiagonal momenta, the percentage of the trapping 
region that might be cut off the corner opposite the gravitational wall by the 
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third centrifugal wall quickly goes to zero and the situation effectively looks like 
the one described in the upper diagrams of Figure 9, unless something unusual 
happens in the evolution of the offdiagonal momenta. 

This is analogous to the argument which shows that collisions with the tilt 
potential should be unimportant as Q ^ oo. Neglecting the velocity due to the 
time dependence of w^, these walls move with unit f2- velocity between collisions 
so the percentage of the trapping region cut off its corner by the tilt wall rapidly 
goes to zero. Put slightly differently, the truncated region quickly disappears on 
a scale which expands to keep up with the trapping region borders. In the type 
IX case, compactness enables one to eliminate the qualifications regarding the 
time dependence of Pa and v'^ in these arguments; the existence of the positive- 
definite quadratic constant of the motion \{2kl)'^n°'''v'^v'i,\ = |n"^e"''+" PaPb\ 
[43] provides an upper bound for the absolute value of each component. One 
may compare the locations of the centrifugal and tilt walls with the locations 
they would have for each symmetric case characterized by the constant upper 
bound of the corresponding component. Due to the smaller value of Pa^ or 
v'g^ compared to its upper bound, the actual wall locations must be outside 
the comparison walls relative to the trapping region. The tilt potentials can 
therefore have less effect than in the comparison symmetric cases where it has 
already been established that they are unimportant in the limit Q ^ oo. Sim- 
ilarly the motion of the centrifugal walls outward towards decreasing values of 
the associated potential is limited by the motion of the comparison wall which 
moves only due to changes in \po \ . Thus the shape of a trapping region in the 
type IX case is determined by one gravitational wall and two centrifugal walls 
in the limit Q ^ oo and the arguments given for the symmetric case show that 
matter is not dynamically important in this limit. This reasoning does not ex- 
tend to the remaining types but the work of Peresetsky [35] seems to indicate 
that the results remain valid. This question requires further study in the present 
approach. 

As SI — > 00, ^-1"" assumes it greatest values during collisions of the uni- 
verse point with the corresponding centrifugal wall, leading to significant values 
of W'^ only during such collisions according to (4.32). This means that the auto- 
morphism matrix S essentially changes only during collisions with a centrifugal 
potential, during which it undergoes a right translation by the matrix e^"" for 
a collision with the potential Uc"'\ where 9 is given by the symmetric case de- 
scription of the collision already discussed. The corresponding component Pa 
does not change much during this collision, justifying the approximation. If 
the outward motion of the centrifugal walls towards decreasing values of their 
potentials due to decreasing \po \ is not compensated for by the time dependence 
of the offdiagonal momenta, then approaches zero. In the type IX case 

the bounds on the offdiagonal momenta lead to the vanishing of the offdiagonal 
velocities and an end to the offdiagonal motion. However, this seems to be an 
effect which docs not depend on the Bianchi type. 

Lifshitz, Khalatnikov and Belinsky [64] claim that the unprimed and primed 
gravitational and fluid supermomenta and the matrix S approach constants 
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as $1 ^ oo. In the nonabclian class A case, all of the centrifugal potentials 
would then become static since the G-angular momenta Pa would approach 
constants. The evolution then would depend only on the arrangement of the 
walls which enclose the trapping region. All of these walls would soon find 
themselves in the regime in which their potentials are all exponential and the 
differences between Bianchi type are washed out. For corresponding trapping 
regions in these nonabclian class A Bianchi type cases, the evolution becomes 
identical in this limit. In particular the change in Kasner axes which occurs 
during a collision with the single gravitational wall as computed by Lifshitz, 
Khalatnikov and Bclinsky [64] docs not depend on Bianchi type. In the class B 
case, excluding type V, Pi and P2 would approach constants but P3 would still 
change during gravitational wall collisions as in the symmetric case, leading to 
an additional change in the Kasner axes relative to the class A models. 

The BLK discussion involves a trapping region formed by the potentials 
Ug^\ Uc^^ and f/i^^; the outward motion of the two centrifugal walls as O — > 00 
leads to their condition ':$> » being satisfied by points in this trap- 
ping region, the inequality increasing with fi. The BLK limit for this trapping 
region (n^^-* ^ 0) is governed by the following supertime time gauge Hamilto- 
nian for the diagonal variables, with Pi, P2, I and v'^ constants related by the 
supermomentum constraints, while P3 is determined by those constraints 

e^^'Pa = -ap+- 2klv'3 . 

(4.39) 

Furthermore, as \po\ and therefore decrease towards zero with increasing 
f2, the term ap+ becomes negligible in the supermomentum constraint so that 
P3 also approaches a constant. In other words the BLK limit for the generic 
dynamics of all spatially homogeneous perfect fluid spacetimes (7 7^ 2) allowing 
anisotropic spatial curvature (all but I and V) is eventually described by a 
Bianchi type II model with frozen in values of the offdiagonal variables, class 
B spacetimes passing through a type IV stage before reaching this limit. The 
relevant type II value of g for this limit in the case n*^^^ 7^ is the matrix Lie 
algebra of superdiagonal matrices, a 3-dimensional Lie algebra of Bianchi type 
II, for which the centrifugal potentials are all of type (iii). 

In order to understand the existence of this BLK limit using the present 
approach, one must examine in detail the equations of motion for the offdiagonal 
momenta and fluid variables. Although the BLK analysis may be correct, it 
would certainly help to have a clearer discussion of these questions. However, 
this requires a more careful treatment than the limitations of time and space 
permit in this article. 
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5 Concluding Remarks 



At one time or another spatially homogeneous cosmological models have cap- 
tured the imaginations of a significant fraction of the relativity commimity. 
Independent of the numerous reasons which have motivated their study, this 
finite dimensional class of spacetimes within the context of a particular rela- 
tivistic theory of gravitation represents an extremely rich mathematical system 
worthy of a careful and systematic treatment. Its description brings together 
the fields of differential geometry, Lie group theory and classical mechanics in 
an elegant example revealing many facets of each individual field and provid- 
ing a finite dimensional setting within which certain questions about relativistic 
theories of gravitation may be more easily explored. 

The present article has endeavored to tie together numerous individual re- 
sults concerning solutions of the Einstein equations of general relativity for this 
class of cosmological models by constructing a single unifying framework based 
on a few simple ideas. Developing the consequences of these ideas unavoidably 
leads to an unusual amount of detail. These details have not been sacrificed in 
the hope of establishing a clearer perspective of this area of cosmology. 

The key to the present discussion of spatially homogeneous dynamics is 
of course diagonal gauge. In the semisimple case of Bianchi types VIII and 
IX, this gauge coincides with both the minimal strain and minimal distortion 
gauges introduced by Smarr and York for spatially compact and asymptotically 
flat spacetimes [69]. For the remaining types where the correspondence with 
the general theory breaks down [2], diagonal gauge generalizes the attractive 
properties of the minimal strain minimal distortion gauge in the semisimple 
case. The decomposition of the usual synchronous gauge gravitational variables 
into diagonal and ofFdiagonal variables which accompanies the transformation 
to diagonal gauge in the semisimple case is merely an application of the ideas 
put forth by Fischer, Marsden and York, among others (see ref.(62) for a sum- 
mary and further references), in describing the true degrees of freedom of the 
gravitational field and related questions. It is therefore clear that the power 
of the present approach has come simply from taking seriously the ideas that 
have been developed in the context of the three-plus-one approach to general 
relativity and modifying them when necessary to fit the spatially homogeneous 
symmetry. 

A The Adjoint Representation Trick 

If e = {ca} is a basis of the Lie algebra g oi a. 3-dimensional Lie group G, 
with structure constant tensor components C"^bc, one may introduce canonical 
coordinates {x"'} of the second kind on the group [52] in a local patch centered 
at the identity (a;° = 0) by the parametrization 

a; = exp(a;^ei) exp(a:^e2) exp(a;'^e3) e G. (^-1) 
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The left invariant and right invariant frames e and e on G are related by the 
matrix adjoint representation of G 

where the expression for R(.t) holds in the given coordinate patch and = 
ade(ea) = {C^ac) are the adjoint matrices for the given basis e, spanning the 
adjoint matrix Lie algebra ade(g') and satisfying [ka,kf,] = C^ab^c which is the 
matrix form of the Jacobi identity. The matrix R belongs to the linear adjoint 
matrix group Ade{G) of G, which for a connected Lie group equals the inner 
automorphism matrix group IAute{g) of the Lie algebra g. When {ka} are 
linearly dependent matrices (i.e. the center of g is trivial) and thus a basis of 
&de{g) having the same structure constant tensor components as e, then a,de{g) 
and g are isomorphic and Ade{G) is locally isomorphic to G, with {a;"} serving 
as the corresponding coordinates on the matrix group through (A. 2). This is 
true for all Bianchi types except I, II and III. 

When C"(,c belongs to Cjj, the adjoint matrices are 



ki = 






- ae^i, 




k2 = 


-n(3)ei3 




- ae^2, 


{A.'S) 


k3 = 






+ a{e\+e\), 





where {e'^i,} is the standard basis of 5^(3, i?), enabling a matrix to be expressed 
as A = A'^h^^a- The matrix R(x) is then easily evaluated 

/l -ax^ \ / C2 n(i)s2\ / 036"^' -71(1)536''^' 0^ 
^{x) = Cl -n(2)si 1 -ax^ n(2)s3e«^' 036"^' 

\0 n(3)si Cl / \-n^^h2 C2 / V ly 

(AA) 

The following abbreviations and identities are useful, where (a, b, c) is any 
cyclic permutation of (1, 2, 3) and no indices are summed 

TO^"^ = (-nWn('=))5 , m(") = (n(*)n(^))5 , c^ = cosh m^^^a;" = cosm^^^a;'' 
s„ = (m(''')-^sinh m^^^x" = (m*"))"^ sinm^^'a;" , lim (c„,s„) = (l,a;") , 

m(»)->0 

{caf - (m(«)s„)2 = {c,f + (m(«)s„)2 = 1 , 

dCa = {vrS"'^)'^ Sadx*^ , dSa = Cadx'^ . 

The canonical adjoint matrix groups (evaluating (A. 4) at the canonical 
points of Cd) are seen to be IX: 5*0(3, R) with k^ generating an active ro- 
tation about the a*'* axis, VIII: 5*0(2. 1) with ka generating an active rotation 
about the 3^'' axis and ki and k2 passive and active boosts along the f** and 
2"'^ axes, respectively, VIIq: IS0{2, R) ~ Euclidean group of the plane and VIq: 
750(1, 1) ~ Poincare group in 2 dimensions. The latter two matrix groups are 



72 



better interpreted in terms of their inhomogeneous action on by letting them 
act on the column vector (jji,y2,l), in which case ka generates an active rotation 
and a passive boost respectively of while {ki, generate the translations. 
The canonical groups for Bianchi types VII^, Ylh^-i, V and IV differ in inter- 
pretation from types VIIq and VIq only regarding the corresponding action of 
k3 on i?^. In the first two cases the rotation or boost e^'*'*^ of Bianchi types 
VIIq and VIo Ik acc;onipanied by a dilation of R^ by the factor e"^ in types 
VII;j and Vlh, while only the dilation is present for type V, and for type IV 
this dilation is instead accompanied by a null rotation of R^ considered as a 
null 2-plane in Minkowski spacetime. In other words all of these matrix groups 
are isometric to 3-dimensional subgroups of the conformal group of Minkowski 
spacetime. The canonical adjoint matrix groups for the remaining types I, II 
and III are respectively trivial(I) or abelian(II) or nonabelian(III) 2-dimensional 
groups of GL{3, R). 

The Bianchi types fall naturally into three categories in this context: (i) 
the types I, II and III, (ii) the remaining nonsemisimplc types and (iii) the 
semisimple types VIII and IX. Consider category (ii) and assume n'^^^ = for 
simplicity. The adjoint matrix group T2 Xg Adz which through its action on 
the column vector {y\_.y2; 1) is identified with the group of translations of R^ 
together with a 1-dimensional subgroup of linear transformations of R? taken as 
a vector space. Here T2 = exp(span{e^i, 6^2}) corresponds to the translations 
and Adz = cxp(span{k3}) to the linear transformations. Carrying out the 
matrix multiplication indicated in (A. 4) 

R{x) = „(2)53e-' cze'^^' -n(^)x^ - ax^ (^.6) 
V 1 / 

and then considering the product R(a;3) = R(a;i)R(a;2) of two such matrices, 

one may explicitly evaluate the the multiplication function xz = (p{xi,X2) in 
these coordinates. Introducing the obvious notation cz,a = cosh m'^^^x^, etc., 
one finds 



X 



„2 



^(03,1X2 - 71*^-^^53,12-2) + x\ 



e''-Hn^''>sz,ixl+cz,ixl)+xl {A.7) 

^3 — ™3 _|_ ™3 

Similarly by permutation of these formulas one can obtain the multiplication 
law for all points of Cd falling in category (ii). With a little eff'ort one can write 

down a complicated generalization valid for case (iii) as well, thus obtaining a 
local expression for the multiplication law on G which is in fact valid for all 
points of Cd. 

For any Lie group the diflFerential of the adjoint matrix satisfies the relations 

R-^rfR = kaw", dH R-^ = kaw" , (^4.8) 

where {w"} and {u)"'} are the invariant 1-form bases dual to e and e respectively. 
When {ka} are linearly independent matrices, these relations may be used to 
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evaluate these l-forms in the local coordinates under consideration and then 
the corresponding bases e and e may be obtained from these expressions using 
duality. The result is 



e {c2C3dx^ + n^^hsdx'^) ei = e"-^' {n^^h3d2 + C3C2 ^ {di — n^^^ 3283 



u)^ = n^^h2dx^ + dx^ 63 = ^3 

u)^ = dx^ + (n^^'s2 — ax^)dx^ e\ = d\ 

= cxdx^ - {vP-UxC2 + ax^)dx^ 62 = 01^2 - rP^ s^c^ "^(^s - -nl^^^ S2dx)) 
(2,3 ^ n'^'sidx^ + c\C2dx^ 63 = n^'^^ S\d2 + C1C2 "^(^3 - n(^^S25i) 

+ aix^dx + x'^di) . 

{A.9) 

These formulas also hold for all points of Cd, providing C^-parametrized ex- 
pressions for the collection of invariant fields on the Cc-parametrized family 
of Lie groups G, with the C/j-parametrizcd local expression for the multipli- 
cation law given by a suitable generalization of (A. 7). If G is assumed to be 
simply connected, then {x"} are global coordinates on G ~ i?^ for all but 
the type IX orbit where G ^ ^ SU{2), and the group G may simply 
be defined by these formulas, with restrictions on the ranges of the coordi- 
nates for points in the type IX orbit obtained by consideration of the basis 
e = — ii(n(^)n(^)) 2(7'2, — iz(n(^'n(^)) 2(73} of the Lie alge- 

bra of 5i7(2),where {Ca} are the standard Pauh matrices. In short, the Cd- 
parametrized simply connected 3-dimensional Lie group Gco has been con- 
structed. 

By slightly modifying all of the above formulas, one may obtain expressions 
for the quantities associated with the multivalued matrix group function G on 
Cd introduced in the third section. Rather than giving formulas valid for all 
points of Cd, attention is restricted to those points for which n^^^ = for clarity. 
Referring to (3.4), make the following definitions 

{Ka} = {-6(2) e32, 5(^^631, -a^'^^e21 + a^^^eU} 

{K} = {-n(2)e32, n^^^eSl, -nWe21 + n(2)el2} (AlO) 
G3 = cos (a(i)a(2))^^3 ^ (a(i)a(2))-^ sin (a^i^af^))^^^ ^ 

and let {£3, §3) be obtained from (C3, S3) by the replacement (n^^^ , x") ^ (n*^") ,0°')- 
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Then one can write down the following formulas by inspection 

\o 1/ Vo 1^ 



w-1 


= C3^1 +n(i).S3^2 


Pi 


= n(2)s3pi + C3P2 




= -n(2)s3^1 + C3^2 


P2 


= C3P1 - n(^^S3P2 






P3 


= P3 






Pi 


= Pl 




= ^2 _ ^(2)^1^3 


P2 


= P2 






P3 





(All) 

It is worth noting that the only values of G which arc not adjoint matrix groups 
of some group are those for which rank n = 1. These are matrix groups of 
Bianchi type II. The centrifugal potentials are all exponential (type (iii)) for 
such values of G. For example, in the case n'^''^ = 0, setting n^^^ = makes 
a(^) = and gives such a value of G; its Lie algebra consists of superdiagonal 
matrices. 

The invariance of the structure constant tensor under the action of the au- 
tomorphism group 

jA(^^c) = C\c , A e Auteig) {A.12) 

may be written in matrix form using the adjoint matrices 

Ak„A-i = kiA^ . (A.13) 

Similarly the tensor 6j'^ = C^ac — ^S^aCic which is related to the matrices Sa 
of (2.49) in the same way C^ac is related to ka is also invariant under the 
automorphism group so the matrices Sa also satisfy (A. 13). This fact was used 
in (3.37). 



B Automorphism Matrix Groups 

Discussion of the matrix group of the Cc-parametrized Lie algebra gc^ leads 
one to consider the four categories of table II labeled by the integers 0, 1, 2, 
3 which represent the orbit dimensions within Cd- Discrete automorphisms 
will be ignored here, being discussed at length in ref. [43], so only the identity 
component Aute{g)~^ C Aute{g) of the matrix automorphism group and the 
identity component SAute{g) = SL{3,R) n Aute{g)~^ of the special automor- 
phism matrix group will be considered. Deqignate their matrix Lie algebras by 
auteig) and saute{g). The adjoint-matrix Lie algebra ade(s') is a Lie subalgebra 
of aute{g) = dere{g) which consists of the matrices of derivations of the Lie 
algebra g. 
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For the first (Abclian) category, the fuU matrix groups are GL{3, R) and 
SL{3, R) respectively, while for the last (semi-simple category) one has Aute{g) = 
SAute{g)~^ = Ade{G) which has been discussed in appendix A. For the third 
category, assume n*^^^ = for uniformity of discussion. Here Auti,{g) is 4- 
dimensional, the extra dimension relative to the adjoint subgroup Ade{G) aris- 
ing from the addition of the matrix I3 = diag(l,l,0) to the basis of ade(g'), 
generating the scaling matrix e^^^ = diag(e^, e^, 1) G Diag{3, For the class 
A types of this category SAute{g) = Ade{G), while, in the class B case, replac- 
ing the generator ka of ade(5f) by kg leads to a basis of aute{g), where {k°} are 
the matrices obtained from (A. 3) by setting the structure constant a to zero. 
Thus for this category SAute{g)~^ = T2 Xs Ad^, where Adg = exp[span{k3}] 
and the group T2 = exp(span{e^i, 6^2}) is generated by span{ki,k2} (except 
for Bianchi type VI_i = III, where ki and k2 arc linearly dependent, satisfying 
aki -|- n''^''k2 = 0, so only a 1-dimensional subgroup of T2 corresponds to an 
adjoint transformation). 

For the remaining category containing Bianchi types II and V, some notation 
is required. Let GL2,3 be the subgroup of GL(3, R) isomorphic to GL2,r which 
leaves the 3rd axis of R^ fixed, let SL2,3 = GL2,3 U SL{3, R) and let be the 
transpose of the matrix group T2: let and be the corresponding Lie algebras. 
Again for uniformity of discussion assume n = diag(0, 0, n'^-') for type 11. For 
Bianchi type V, one has Aute{g) = T2 Xg GL2,3 and S{g) = T2 XsSL2 while for 
type II SAut,,{g) = Aute{g)^ is obtained by adding the diagonal 

automorphism generator diag(l, 1, 2) to a basis of the Lie algebra of SAute{g), 
generating the scaling diag(e^, e^, e^^) e Diag{3, J?)+. 

Introduce the map (homomorphism) ad : aut{G) o.ut{g) by ad(^)X = 
[^,X] for ^ e aut{G) and X G g, and let ade(^) be the matrix of ad(^) with 
respect to the basis e of g; restricting the domain of ad to g gives the adjoint 
representation of g. For a simply connected Lie group G this map ad is a 
Lie algebra isomorphism; one can therefore, consider the inverse map ad~^ : 
a.d{g) — > aut{G) which associates a generating vector field of the automorphism 
group of G with each matrix of the matrix Lie algebra of the automorphism 
group of g. For example, the adjoint matrices map onto elements of ad(G') with 
ad~^ (ka) =60 — 60, the coordinate representation of which may be read off from 
(A. 9). For the nonsemisimplc Bianchi types where aut(g) is larger than ad(G), 
the coordinate representation of the remaining automorphism generators is also 
easily obtained. For the third category only one additional linearly independent 
such generator exists given by ad~^{l3) = —{x^di +x'^d2) for those points of Cd 
for which n*^"^) = 0. For the first (Abehan) category where aute{g) = gl{3, R) 9 
B, one has ade(B) = ~B°-i,x^da- The same formula holds for Bianchi type V if B 
belongs to the Lie algebra 5/2,3, leaving only the remaining Bianchi type II. Here 
only the Lie subalgebra 5/2,3® span{diag(l, 1, 2)} C autf,{g) need be considered, 
assuming n = diag(0, 0, n^^^) for simplicity. Again the Abelian formula holds for 
all diagonal elements of this Lie subalgebra, leaving only off-diagonal elements 
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to be considered. Here the formula 

ad;\-q^'^e\+q^^^e\) = q^'^x^d^ - q^^h'd2 + ^n^^'^iq^^Kx'f - q^'\x^)^]d3 

{B.l) 

is essentially due to Bianchi [10]. An important thing to understand about the 
automorphism group is how it acts on the dual space Qq^ of left invariant 1- 
forms on Cd ■ This is important since both the supcrmomentum 'H^co"' and fluid 
current 1-form VaUj"^ transform under this action, namely 

Sauj" Gg*^ ZbA-^\uj^ , A e Aut^ig) . {B.2) 

The orbit space then gives one information about incquivalcnt initial data, while 
the orbits themselves give information about constants of the motion if the 
action is intransitive [43] . This automorphism group action is also important in 
establishing a canonical form for the Kasncr axes in the BLK limit [35, 64]. 

Define z = \n°'^ZaZ},\^/'^ and 7"'' = sgn(detn)n''''. For the semi-simple types 
IX and VIII, 7"'' are the components of an inner product of signature (+ + +) 
and ( — h +), respectively, and the norm z associated with this inner product is 
an invariant of the automorphism action. The orbits are just surfaces of constant 
z. In the canonical case these are just the origin-centered spheres of radius r 
for type IX and pscudosphcrcs of radius z for type VIII. For the latter type, 
£ = sgn{'y°'^ ZaZb) may assume the invariant values 1, and -1 for the two simply 
connected timelike hyperboloids, the null cone and the nonsimply connected 
spacclike hyperboloids, respectively. (The timclikc directions arc associated with 
the rotations and the spacelike directions with the boosts of S02,i, while the 
null directions are associated with the null rotations.) For the nonsemi-simple 
types; n°'^ are the components of a degenerate quadratic form. The norm z is a 
constant only for the special automorphism action. For the canonical type- VII 
point, the orbits of SAutf.{g) are cylinders of radius z about the 2;3-axis which 
consists of fixed points, while for the canonical type- VI point each value of z 7^ 
consists of 4 disconnected orbits ('hyperbolic cylinders') equivalent under the 
action of discrete automorphisms, with z = consisting of 4 disconnected half- 
plane orbits meeting at the z^ axis of fixed points. Under the action of Aut^^g) 
all points with z ^ are equivalent. For type-IV points with n*^'^ ^ 0, the orbits 
of SAute{g) are the planes of constant \zi \ for 2; ^ 0, but the lines parallel to the 
23-axis in the plane zi = for z = 0. For typc-II points with rS^^ ^ 0, SAute{g) 
acts on each plane of constant [zal as the inhomogeneous special linear group of 
the plane, with all points Of Z3 ^ being equivalent under Autf.{g), but only 
as the special linear group on the plane Z3 = with all points of this plane 
except the origin equivalent under Autf.{g). For type V the Z3-axis consists of 
fixed points while all other points belong to the same orbit for both Aute{g) 
and SAute{g). For type I only two orbits exist for both groups, the origin and 
all other points. 
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C A Spatially Homogeneous Perfect Fluid 

Following the notation of chapter 23 of [18], consider a spatially homogeneous 
perfect fluid with energy density p, pressure p, 4- velocity field u; baryon number 
density n, chemical potential /i = {p + p)/n and synchronous-gauge energy- 
momentum tensor 

T''0^{p + p)u"u0+pS''(3 . (C.l) 

For the equation of state p = (7 — I)/?, n may be taken as the single independent 
thermodynamic variable, in terms of which the others may be expressed as 
follows 

p = n^ , p=i"/- l)n^ , n = 771^-1 , (C.2) 

where a constant of integration has been eliminated by a redefinition of n. Fol- 
lowing Taub [23,70], introduce the spatial circulation 1-form v with components 
Va = fJ-Ua and the spatial scalar density 

e = ng^/\-^ = niJ,-^g^/\-^ , where u-^ = (1 + Waw")^/^ . 

Under the change of basis Ca = A~^^a&b, the fluid variables {£, Va) transform in 
the following way 

{i,Va) = .fA{e,Va) = {\det{A)\-'£,vt,A-'i) . (C.3) 

The matter super-Hamiltonian and supermomentum may be expressed in terms 
of these variables as follows 

= 2Mv^ - 2fcpgi/2 ^ j^M ^ 2k£va . (C.4) 

The matter super-Hamiltonian, considered as an independent function of the 
spatial metric and the matter variables {n,i,Va), satisfies (2.50) and, therefore, 
acts as a potential for the matter driving force T* . 

The equations of motion for the fluid variables {£,Va) in almost synchronous 
gauge are [43] 

£ = N£(v^)-'2a%, , Va = N{v^)-'C\av' , (C.5) 

while the defining relation for £ may be used as an equation implicitly defining 
n in terms of (gab, £,Va)- The fluid constants of motion in almost synchronous 
gauge are described elsewhere [43]. In the class A case, for example, I and 
V = \n"'''vaVi)\^^^ are constants of the motion, while £^-^/^V is a class B constant 
of the motion; others exist, however, with at least two nontrivial constants of 
the motion of the motion in all cases. 

D The Kantowski- Sachs models. 

The only spatially homogeneous space-times not described by a Bianchi type 
model are the Kantowski-Sachs space-times [71,72] which have a 4-dimensional 
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isomctry group acting transitively on the hypersurfaces of homogeneity (im- 
plying local rotational symmetry) but no 3-dimensional subgroup which acts 
simply transitively on these hypersurfaces. The spatial metrics of such space- 
times were also classified by Bianchi [10] in his categorization of all Riemannian 
3-manifolds which admit a Lie group of isometries. Their isometry group is the 
direct product group R x SO {3, R) acting on the 3-manifold R x S2 (in the 
simply connected case) , where the additive group of real numbers acts on R by 
translation and the special orthogonal group acts isometrically on 5*2 with its 
standard metric. 

To appreciate the relationship of these models to the Bianchi type models, 
consider the following choice of Euler angle coordinates for the class A subman- 
ifold of Cd for which n^^) 7^ 0: 

X = cxp(t/^e3) cxp(y'e2) cxp{j/e;i) G G . (-C-1) 

Using the trick of appendix A, one finds the coordinate expressions for the 
invariant fields to be 

= n(i)(s3rfy' - C3S2,idy') , ei = {n^'yy'Um^^yfssdi 

- (S2,l)"^(92 - C2,l93)] , 

= c^dy'^ + (m'3))^S3S2,irfj/^ , 62 = 0381 + S3(s2,i)~^(92 - 02,193) , 

= C2,idy^ + dy^ , es = 83 , 

0)1 = n(i)(-S3,2rfy' + C3,2S2,idy^) , ei = in^^^-^[-{m^^^f 83,281 

+ C3,2(S2,l)"^(93 - C2,l92)] , 
= C3,2dy^ + (m^^^ fs3,2S2,\dy^ , 62 = 03,281 + S3,2(S2,l)"^(93 - 02,1^2) , 



{D.2) 

where the notation of (A. 5) and (A. 7) is used with the replacement of x by y. 
Note that these coordinates are singular at the identity (y' = 0). 

The left coset space X = G/ exp(span{e3}) is obtained by identifying points 
of G along integral curves of 63 = 83 (the orbits of right translation by the sub- 
group exp(span{e3}), namely the y'^-coordinate lines of these local coordinates 
which are comoving with respect to 63). are local coordinates on X 

which reduce to standard spherical coordinates {0, 0} on S2 at the canonical 
type- IX point of C^j. The right invariant vector fields Ca, since they are invari- 
ant along 63, project to fields = ^{oa) on the quotient space X obtained by 
ignoring their third components in these local coordinates, {^a} is the image 
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basis of generators of the natural left translation action of G on X 

6 = (n(^))-M-(m(3^)'s3,29l - S3,2C2,l(S3,l)"'a2 , 
6 = C3,29i - S3,2C2,l92 , 
^3 = 1^2 , 



(^•3) 



The isotropy group at the identity coset exp(span{e3}) € X of this left 
translation action of G is just exp(span{e3}). Now consider the following left 
invariant second-rank symmetric covariant tensor ^g^ on G, with k = m^^^ : 

= (m(2))2[rfj/i dy^ + (m(3)s2,i)^ V ® ^2/^] • ^^'^"^ 

This is also invariant under right translation by the subgroup exp(span{e3}) and 
so projects to a left invariant tensor on X . For (to^''^)^ > this is a Riemannian 
metric on X of constant Gaussian curvature -R^^i2 = —k^. (For (m^^^)'^ < 
it is a pseudo-Riemannian metric of constant curvature, but at m^^^ = it is 
degenerate. 

If Riemannian metrics are of interest, one might as well consider only the 
line segment 

T = {{a,n}€CD\a = 0; n = diag(l, 1, n^^^), n^^^ehl,!]} {D.5) 

connecting the canonical type- VIII,- VIIq, and -IX points of Cij. The most gen- 
eral Riemannian metric on the product manifold R x X invariant under the 
natural left action of the direct-product group R, x G (where R is the additive 
group of real numbers, with coordinate u, and only the groups G parametrized 
by the line segment T are considered) is 

^Ek = 522 du0du + 522 ^Sk , (^-6) 

where ^22 ^.nd ^33 are constants and now = —n^^\ This is in fact the 
class of metrics studied by Kantowski and Sachs. However, when n^^^ < 0, 
these metrics are locally rotationally symmetric Bianchi type metrics, as shown 
by Bianchi [10]. The case n^^^ = is obviously just the locally rotationally 
symmetric type-I or type-VIIo metric expressed in cylindrical coordinates. For 
n^^-* < 0, three inequivalent classes of standard coordinate systems exist for the 
constant-negative-curvature 2-manifold X [73,74], namely coordinates chosen 
to be comoving with respect to ^1, ^3 or ^1 -|- ^3 which correspond, respectively, 
to boosts, rotations and null rotations in the simply-connected covering group 
502,1 of 502,1 which is the value G assumes at n^^^ = —1. The coordinates 
{y'l y^} = {^2, fl2} are comoving with respect to ^3. For the canonical type- VIII 
case n*^^-* = —1, choosing new coordinates {x\x'^} = {vz^rj^} as in (2.34) of 
ref. [73] which are comoving with respect to ^1 -I- ^3 and defining = u, one 
has 

^Si = 922<^^ fSidx"^ + 533 (c&^ dx^ -h e~^^^dx^ ® dx^) , {D.7) 
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which has the component matrix g = diag(g33,p22i533) ^ A^t(2) with respect 
to the type-Ill = VIi frame e of (3.52) evaluated at g = a = 1/2, i.e. this is 
just the locally rotationally symmetric Bianchi type-III metric, whose spatial 
curvature matrix (in the same frame) is 

533R = -K2diag(l, 0, 1) , K = l. {D.8) 

The Einstein equations for a spatially homogeneous perfect-fluid spacetime 
with metric 

= N^dt ® dt , e [-1, 1] , {D.9) 

differ from this case only in that one retains the factor in the spatial curva- 
ture. The value k = corresponds to the locally rotationally symmetric typc-I 
or typc-VIIo, case, while the complex rotation k = I ^ k = i takes one to the 
Kantowski-Sachs case. The equations for (-/V, 322) 533) ^■'^d the fluid variables 
in the Kantowski-Sachs case k = i are identical with the k = 1 case with the 
exception that the sign of the spatial curvature changes. In other words, if 
one solves the locally rotationally symmetric type-III equations with k left as 
an arbitrary parameter, one may obtain the Kantowski-Sachs solutions by the 
analytic continuation k = 1 ^ k = i [74]. The Kantowski-Sachs case could, 
therefore, be included in table II connected by horizontal dots to Bianchi type 
III. The Kantowski-Sachs metric is also related to the locally rotationally sym- 
metric Bianchi type-IX metric by a contraction of the group action in which 
the length of the subgroup exp(span{e3}) ~ becomes infinite; a similar con- 
traction leads from the locally rotationally symmetric Bianchi type- VIII metric 
to the locally rotationally symmetric type-III metric. The singular transforma- 
tion of the vacuum solutions induced by these group contractions is described 
explicitly in ref. [74]. 

References 

[1] R.T. Jantzen, Commun. Math. Phys. 64, 211 (1979). 
[2] R.T. Jantzen, Ann. Inst. H. Pomcare A33, 121 (1980). 
[3] R.T. Jantzen, Nuovo Cim. B55, 161 (1980). 

[4] M.P. Ryan, Jr., Hamiltonian Cosmology, Lecture Notes in Physics, Vol. 13 
(Springer- Verlag, Berlin, 1979). 

[5] M.P. Ryan, Jr. and L.C. Shepley, Homogeneous Relativistic Cosmologies 
(Princeton University Press, Princeton, 1975). 

[6] M.A.H. MacCallum, Cosmological models from a geometric point of view, 
in Cargese Lectures in Physics Vol. 16, edited by E. Shatzman (Gordon 
and Breach, New York, 1973). 



81 



[7] M.A.H. MacCallum, Anisotropic and inhomogeneous relativistic cosmolo- 
gies, in General Relativity: an Einstein Centennial, edited by S.H. Hawking 
and W. Israel (Cambridge University Press, Cambridge, 1979). 

[8] M.A.H. MacCallum, The mathematics of anisotropic spatially homoge- 
neous cosmologies, in Physics of the Expanding Universe, edited by M. 
Demianski (Springer- Verlag, Berlin, 1979). 

[9] M.A.H. MacCallum, Rc^lativistic cosmology for astrophysicists, lectures 
given at the VII International School of Cosmology and Gravitation (Erice, 
Sicily, 1981): in The Origin and Evolution of Galaxies, edited by de Sab- 
bata (World Scientific, Singapore, 1982). 

[10] L. Bianchi, Sugli spazi a tre dimensions the ammettono un gruppo continuo 
di movimenti (1897), in Opere, Vol. 9 (Edizione Cremonese, Roma, 1958); 
English translation with editorial notes by R.T. Jantzen in Gen. Rel. Grav. 
33, 2001. 

[11] A.H. Taub, Ann. Math. 53, 472 (1951); reprinted with editorial notes by 
M.A.H. MacCallum in Gen. Rel. Grav. 34, 2002. 

[12] 0. Heckman and E. Schucking, Relativistic cosmology, in Gravitation: an 
Introduction to Gurrent Research, edited by L. Witten (J. Wiley & Sons, 
New York, 1962). 

[13] F. B. Estabrook, H. D. Wahlquist and C. G. Behr, J. Math. Phys. 9, 497 

(1968) ; unpublished notes by Behr to appear in Gen. Rel. Grav. 34, 2001. 

[14] G.F.R. Ellis and M.A.H. MacCallum, Gommun. Math. Phys. 12, 108 

(1969) . 

[15] C.W. Misner, Phys. Rev. Lett. 22 1074 (1969); Phys. Rev. 186, 1319 (1969). 

[16] C.W. Misner, Classical and quantum dynamics of a closed Universe, in 
Relativity, edited by M. Carmeli, S.I. Fickler and L. Witten (Plenum Press, 
New York, 1970). 

[17] C.W. Misner, Minisuperspace, in Magic without Magic, edited by J. R. 
Klauder (Freeman, San Francisco, 1972). 

[18] C.W. Misner, K. S. Thorne and J. A. Wheeler: Gravitation (Freeman, San 
Francisco, 1973), Chapter 30. 

[19] M.P. Ryan, Jr., J. Math. Phys. 10, 1724 (1969). 

[20] M.P. Ryan, Jr., Ann. Phys. (N.Y.) 65, 506 (1971). 

[21] M.P. Ryan, Jr., Ann. Phys. (N.Y.) 65, 68, 541 (1971). 

[22] M.P. Ryan, Jr., A.R. Mosier and R.A. Matzner, Ann. Phys. (N.Y.) 79, 558 
(1973). 



82 



[23] M.P. Ryan, Jr., J. Math. Phys. 15, 812 (1974); A.H. Taub and M.A.H. 
MacCallum, Commun. Math. Phys. 25, 173 (1972); G.A. Sneddon, J. Phys. 
A9, 229 (1976). 

[24] M.P. Ryan, Jr., Ann. Phys. (N.Y.) 70, 301 (1972). 

[25] Y.A. Bclinsky and I.M. Khalatnikov, Sov. Phys. J.E.T.P. 29, 911 (1969). 

[26] Y.A. Belinsky, I.M. Khalatnikov and E.M. Lifshitz, Ann. Phys. (N. Y.) 19, 
225 (1970). 

[27] Y.A. Belinsky, I.M. Khalatnikov and E.M. Lifshitz, Sov. Phys. J.E. T.P. 

32, 173 (1971). 

[28] Y.A. Bclinsky, I.M. Khalatnikov and E.M. Lifshitz, Sov. Phys. J.E. T.P. 

33, 1061 (1971). 

[29] I.M. Khalatnikov and V.L. Pokrovski, A contribution to the theory of 
homogeneous Einstein spaces, in Magic without Magic, edited by J. R. 
Klauder (Freeman, San Prancisco, 1972). 

[30] S.P. Novikov, Sov. Phys. J.E.T.P. 35, 1031 (1972). 

[31] 0. Bogoyavlensky and S.P. Novikov, Sov. Phys. J.E.T.P. 37, 747 (1973). 

[32] 0. Bogoyavlensky, Sov. Phys. J.E.T.P. 43, 187 (1976). 

[33] 0. Bogoyavlensky and S.P. Novikov, Buss. Math. Surv., 31, 31 (1976). 

[34] 0. Bogoyavlensky and S.P. Novikov, Qualitative Theory of Homogeneous 
Cosmological Models, in Proceedings of the I. G. Petrovsky Seminar, Vol. 7 
(1973), English translation in Sel. Math. Sov. 2 (1982). 

[35] A. A. Peresetsky, Russ. Math. Notes 21, 39 (1977). 

[36] C.B. Collins, Commun. Math. Phys. 23, 137 (1971). 

[37] C.B. Collins, Commun. Math. Phys. 39, 131 (1974); C.B. Collins and 
G.F.R. ELLIS, Phys. Rep. 56, 65 (1979). 

[38] C.B. Collins and S.W. Hawking, Astron. J. 180, 317 (1973). 

[39] M.A.H. MacCallum, Commun. Math. Phys. 20, 51 (1971). 

[40] A.G. Doroshkevich, Y.N. Lukash and I.D. Novikov, Sov. Phys. J.E. T.P. 
37, 739 (1974); L.P. Grishchuck, A.G. Doroshkevich and Y.N. Lukash, 
Sov. Phys. J.E.T.P. 34, 1 (1972). 

[41] R.T. Jantzen, Ann. Phys. (N.Y.) 127, 302 (1980). 

[42] R.T. Jantzen, J. Math. Phys. 28, 1137 (1982). 

[43] R.T. Jantzen, Ann. Phys. (N.Y.) 145, 378 (1983). 



83 



[44] B.S. DeWitt, Phys. Rev. 160, 1113 (1967). 

[45] H. Goldstein, Classical Mechanics (Addison- Wesley, Reading, MA, 1950). 

[46] R. Abraham and J.E. Marsdcn, Foundations of Mechanics, revised second 
edition (Benjamin, New York, 1980), Chapter 4. 

[47] A.E. Fischer and J.E. Marsden, J. Math. Phys. 13, 546 (1972). 

[48] L. D. Landau and E.M. Lifshitz, The Classical Theory of Fields (Addison 
Wesley, Reading, MA, 1971). 

[49] F. W. Warner, Foundations of Differenti,ahle Manifolds and Lie Groups 
(Scott, Foresman and Co., Clcnvicw, IL, 1971). 

[50] R. Gilmore, Lie Groups, Lie Algebras and some of their Applications (Wi- 
ley, New York, 1974). 

[51] B. G. Wybourne, Classical Groups for Physicists (Wiley, New York, 1974). 

[52] A. A. Sagle and R.E. Walde, Introduction to Lie Groups and Lie Algebras 
(Academic Press, New York, 1973). 

[53] A.E. Fischer, Superspace, in Relativity, edited by M. Carmeli, S.I. Fickler 
and L. Witten (Plenum Press, New York, 1970). 

[54] E.M. Lifshitz and LM. Khalatnikov, Sov. Phys. J.E.T.P. 12, 558 (1961). 

[55] V. Belinsky and M. Francaviglia, Gen. Rel. Grav. 16, 1189 (1984). 

[56] M. Henneaux, Phys. Rev. D21, 857 (1980); Ann. Inst. Henri Poincare 
A34, 329 (1981). 

[57] 0. Obregon and M.P. Ryan, Jr., J. Math. Phys. 22, 623 (1981). 

[58] S.T.C. Siklos, Commun. Math. Phys. 58, 255 (1980); Phys. Lett. A76, 19 
(1980); — J. Phys. A14, 395 (1981); Gen. Rel. Grav. 13, 433 (1981). 

[59] E.P. Belasco and H.C. Ohanian, J. Math. Phys. 10, 1503 (1969). 

[60] D. Christodoulou and M. Francaviglia, The geometry of the thin-sandwich 
problem, Proc. S.I.F., Course LXVII, edited by J. Ehlers (North Holland, 
Amsterdam, 1979), p. 480. 

[61] G.F.R. Ellis and A.R. King, Commun. Math. Phys. 38, 119 (1974). 

[62] J. Isenberg and J. Nester, in General Relativity and Gravitation, Vol. 1, 
edited by A. Held (Plenum Press, New York, 1980). 

[63] S.W. Hawking, Mon. Not. R. Ast. Soc. 142, 129 (1969); C.B. Collins and 
S.W. Hawking, Mon. Not. R. Ast. Soc. 162, 307 (1973). 



84 



[64] V.A. Belinsky, I.M. Khalatnikov and E.M. Lifshitz, Adv. Phys. 31, 639 
(1983). 

[65] R.A. Matzner and D.M. Chitre, Commun. Math. Phys. 22, 173 (1971). 
[66] V.A. Belinsky and I.M. Khalatnikov, Sov. Phys. J.E.T.P. 36, 591 (1973). 
[67] H. Sirousse-zia, Gen. Rel. Grav. 14, 751 (1982). 

[68] R.A. Matzner, L.C. Shepley and J.B. Warren, Ann. Phys. (N. Y.) 57, 401 
(1973). 

[69] L. Smarr and J.W. York, Jr., Phys. Rev. D17, 2529 (1978). 

[70] A. H. Taub, Proceedings of the 1967 Colloque on "Fluids et champ gravi- 
tationel en relativite generate," No. 170 (Centre National de la Recherche 
Scientifique, Paris, 1969), p. 57. 

[71] R. Kantowski and R.K. Sachs, J. Math. Phys. 7, 443 (1966); see also the 
reprint with editorial notes: R. Kantowski, "Some Relativistic Cosmological 
Models," Gen. Rel. Grav. 30, 1665 (1998). 

[72] C.B. Collins, J. Math. Phys. 18, 2116 (1977). 

[73] E. G. Kalnins and W. Miller, Jr., J. Math. Phys. 18, 1 (1977). 

[74] R.T. Jantzen, Nuovo Gim. B59, 287 (1980). 

[75] [some later work in this approach: 

R.T. Jantzen, Phys. Rev. D33, 2121 (1986); Phys. Rev. D34, 424 (1986); 
Phys. Rev. 35, 2034-2035 (1987); K. Rosquist and R.T. Jantzen, Phys. 
Rep. 166, 89 (1988); R.T. Jantzen Phys. Rev. D37, 3472 (1988); C. Uggla, 
K. Rosquist and R.T. Jantzen, Phys. Rev. D42, 404 (1990); R.T. Jantzen 
and C. Uggla, Ge n. Rel. Grav. 24, 59 (1992); J. Math. Phys. 40, 353 (1999). 



on-line abstracts: b.ttp : //www . homepage . villanova . edu/ robert . j antzen 



85 



